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i. 

This  report  documents  the  work  that  was  completed  to  develop  the  software  tools  and 
geometric  data  base  that  will  help  Air  Force  Design  Engineers  take  aircraft  and  mis¬ 
siles  from  concept  to  test  specimen.  Navgraph  is  a  geometric  modeling  program  that 
allows  the  user  to  define  complex,  convex  and  non-convex  geometric  solid  objects  for 
use  in  the  various  modules  developed.  A  summary  of  the  capabilities  are  given  here: 

NAVGRAPH  can  be  used  to  create  geometric  models  of  aircraft,  weapons, submunitions, 
pylons,  racks,  etc.  to  whatever  detail  is  necessary  for  the  problem  being  studied. 

NAVGRAPH  can  be  used  to  recall  any  geometric  model(s)  from  the  existing  data  base 
(via  DBMERGE)  to  define  the  desired  complex  model  to  whatever  detail  is  necessary. 

The  existing  data  base  consists  of  the  models  developed  under  this  contract,  mainly, 
an  F— 15  aircraft,  fuel  tank,  racks,  rails,  bombs,  missiles,  and  submuntlons. 

CALIPER  can  be  used  to  study  physical  fit  compatibility.  For  example,  CALIPER 
can  calculate  the  separation  distance  between  a  missile  fin  and  an  aircraft  wing  or 
detect  and  calculate  the  interference  distance  and  direction  of  internal  components 
with  respect  to  the  interior  of  a  missile  body.  The  information  derived  here  can 
be  used  to  translate  the  interfering  object  and  remove  the  Interference. ) 
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13.  ABSTRACT  (CONCLUDED) 


MASSPROP  will  calculate  the  mass  properties  (including  surface  area,  volume, 
mass,  center  of  gravity  and  moments  of  inertia)  of  a  NAVGRAPH  model. 

PACKER  (  along  with  OPTDES)  can  be  used  to  package  submunitions  and/or 
internal  components  in  a  weapon  with  constraints  on  certain  mass  properties. 

For  example,  find  the  optimal  placement  and  orientation  of  internal  components 
and/or  submunitions  that  satisfy  constraints  on  the  roll,  pitch  and  yaw  moments 
of  inertia  as  well  as  constraints  on  the  location  of  the  center  of  gravity  with 
respect  to  the  center  of  pressure. 

NAVGRAPH  can  be  used  to  generate  finite  elements  (polygonal  facets)  on  geo¬ 
metric  models  for  use  in  aerodynamic,  radar  cross  section,  stress,  etc.  prediction 
software  codes. 

NAVGRAPH  can  accept  finite  element  data  from  other  Air  Force  data  Dase  formats 
such  as  ASE  and  VUMADM. 

NAVGRAPH  can  be  used  to  generate  computer  color  graphic  displays  (line  drawings 
or  smooth  rendering)  of  all  the  above  upon  demand  including  interference  highlighting. 

In  addition,  formating  software  can  be  used  to  translate  (via  IGES)  the  data 
base  to  and  from  formats  compatible  with  other  software  codes  for  generating  computer 
controlled  milling  machine  instructions. 


Table  of  Contents 

Report  Of  Documentation  Page  (DD  1473)  i 

Acknowledgements  ii 

Table  of  Contents  iii 

Report  Overview  1 

Part  1:  The  Technical  Description  4 

Graphics,  Mass  Properties  and  Hardware  Compatibility  5 

Hardware  compatibility  5 

Mass  properties  calculator  5 

Highlight  interferences  7 

Geometric  data  base  methodologies  8 

CALIPER:  Separation  and  Interference  Distances  10 

The  Solid  Modeling  Convention  1 1 

Separation  Distance  Calculation  and  Interference  Detection  14 

Interference  Distance  Calculation  15 

CALIPER  Examples  16 

PACKER  with  rotations  and  translation  18 

Geometric  modeling  convention  19 

Rotational  degrees  of  freedom  and  the  quaternion  19 

The  optimal  packaging  problem  20 

PACKER  with  spatial  relationships  and  mass  properties  24 

Geometric  modeling  convention  24 

Rotational  degrees  of  freedom  and  the  quaternion  25 

Mass  properties  for  the  packaging  problem  25 

The  optimal  packaging  problem  26 

Examples  of  packaging  problems  solved  with  PACKER  29 

Capabilities  added  to  increase  speed,  robustness  and  flexibility  32 

Explicit  gradients  32 

NAVGRAPH  group  commands  34 

NAVGRAPH  ployfll  text  font  •  36 

Methodologies  to  create  AEROHEAT  shapes  with  NAVGRAPH  36 
Fast  surface  meshing  capability  36 

MOVIE. BYU  software  supplied  36 

GRG  algorithm  modifications  38 

Solid  groupings  into  objects  38 

Packaging  by  layers  3S 

Other  heuristic  methods  and  trouble  shooting  39 

iii 


Translator  Software  Development  40 

TEKNICAD  To  NAVGRAPH  Translator  40 

Translators  Between  NAVGRAPH  and  VUMADM  40 

IGES  Translators  Between  NAVGRAPH  and  ANVIL-5000  41 
Translator  for  AFATL/ASE  43 

Data  Base  Generation  with  NAVGRAPH  44 

Aircraft  Data  Base  Generation  44 

Weapon  Data  Base  Generation  47 

Submunition  Data  Base  Generation  50 

Summary  51 

Part  2:  The  Users  Manual  and  Tutorial  52 

Users  Manual  53 

New  NAVGRAPH  Commands  54 

Eglin  AFB  Data  Base  Translators  76 

Tutorial  77 

AEROHEAT  shape  .RUN  files  78 

MAVERICK  .RUN  file  87 

CALIPER:  Physical  Fit  Compatibility  Study  Example  91 

Optimal  Packaging  Example  with  PACKER  94 

Trouble  Shooting  Guide  For  Model  Generation  129 


Report  Overview 

This  report  discusses  the  work  completed  under  Air  Force  SBIR  contract  F08635- 
S8-C-0063  titled,  ”An  Interactive,  Optimization-Based,  Computer  Graphics  Software 
Package  -  Phase  II”  by  APTEK,  Inc.  Part  1  contains  a  technical  description  of  the 
work  done.  Part  2  contains  the  users  manual  and  tutorial  with  examples. 

The  overall  objective  for  the  Phase  II  work  was  to  develop  and  deliver  the  software 
tools  that  will  help  the  Air  Force  aircraft  and  missile  designer  take  a  weapon  from 
concept  to  test  specimen  quickly  and  accurately  and  to  provide  an  initial  data  base  of 
important  existing  designs. 

The  software  package  consists  of  the  following  modules  described  here.  Figure  1 
shows  the  inter-relationships  of  the  modules. 

•  NAVGRAPH  is  the  geometric  modeling  module  that  allows  the  user  to  define 
points,  lines,  surfaces  and  solids  that  are  used  to  define  the  solid  objects  that  make 
up  a  physical  fit  compatibility  study  or  optimal  packaging  problem.  NAVGRAPH 
also  has,  among  other  things,  a  friendly  user  interface  and  finite  element  pre-  and 
post-processing  capabilities.  Examples  of  the  modeling  capability  are  contained 
throughout  this  report. 

•  CALIPER  is  the  separation  and  interference  calculator.  It  accesses  a  previously 
defined  NAVGRAPH  data  base  (that  contains  a  model  developed  for  a  physical  fit 
compatibility  study)  and  then  calculates  the  needed  separation  and  interference 
distances  between  user  defined  objects.  CALIPER  can  be  run  in  NAVGRAPH  so 
one  does  not  have  to  exit  a  NAVGRAPH  session  to  detect  and  display  interfer¬ 
ences. 

•  PACKER  is  the  optimal  packaging  software.  It  1)  accesses  a  previously  defined 
model  in  a  NAVGRAPH  data  base,  2)  calculates  the  initial  separations  and  in¬ 
terferences  between  all  objects  and  all  objects  and  container  walls,  3)  computes 
the  composite  mass  properties  for  the  whole  model,  4)  allows  the  user  to  "setup” 
the  desired  packaging  problem  by  defining  variables  and  functions  for  the  opti¬ 
mization,  and  5)  allows  the  user  to  run  OPTDES  to  find  the  optimal  packaging 
solution  for  the  defined  problem. 

•  OPTDES,  a  commercial  optimization  software  package,  is  the  optimization  soft¬ 
ware  used  to  solve  the  packaging  problem.  It  has  several  robust  and  efficient 
mathematical  optimization  search  routines  that  help  find  the  values  of  variables 
defining  optimal  design  solutions. 

•  MASSPROP  calculates  the  mass  properties  of  NAVGRAPH  models. 

•  DBMERGE  allows  the  user  to  merge  models  from  several  NAVGRAPH  data  base 
files  together  into  one  file  for  a  physical  fit  compatibility  study  or  optimal  pack¬ 
aging  problem. 
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Figure  1:  Software  Overview 


The  following  capabilities  are  available  to  the  user  of  the  software  package  outlined 
above: 


1.  NAVGRAPH  can  be  used  to  create  geometric  models  of  aircraft,  weapons,  submu¬ 
nitions,  pylons,  racks,  etc.  to  what  ever  detail  is  necessary  for  the  problem  being 
studied. 

2.  NAVGRAPH  can  be  used  to  recall  any  geometric  model(s)  from  the  existing  data 
base  (via  DBMERGE)  to  define  the  desired  complex  model  to  what  ever  detail 
is  necessary.  The  existing  data  base  consists  of  the  models  developed  under  this 
contract,  namely,  an  F-15  aircraft,  fuel  tank,  racks,  rails,  bombs,  missiles,  and  sub¬ 
munitions.  One  can  use  the  existing  data  base  and  DBMERGE  to  build  complex 
models  by  assembling  simpler  models  without  re-entering  their  data  by  adding 
new  subassemblies,  or  by  modifying  old  data. 

3.  CALIPER  can  be  used  to  study  physical  fit  compatibility.  For  example,  CALIPER 
can  calculate  the  separation  distance  between  a  missile  fin  and  an  aircraft  wing  or 
detect  and  calculate  the  interference  distance  and  direction  of  internal  components 
with  respect  to  the  interior  of  a  missile  body.  The  information  derived  here  can 
be  used  to  translate  the  interfering  object  and  remove  the  interference. 

4.  MASSPROP  calculates  the  mass  properties  (including  surface  area,  volume,  mass. 
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center  of  gravity  and  moments  of  inertia)  of  a  NAVGRAPH  model.  These  func¬ 
tions  can  be  used  as  design  functions  to  match  or  limit  user  defined  design  criteria. 

5.  PACKER  (along  with  OPTDES)  can  be  used  to  package  submunitions  and/or 
internal  components  in  a  dispenser  missile  cargo  bay  with  constraints  on  certain 
mass  properties.  For  example,  it  will  find  the  optimal  placement  and  orientation, 
of  internal  components  and/or  submunitions,  that  satisfies  constraints  on  the  roll, 
pitch  and  yaw  moments  of  inertia  as  well  as  constraints  on  the  location  of  the 
center  of  gravity  with  respect  to  the  center  of  pressure. 

6.  NAVGRAPH  can  be  used  to  generate  finite  elements  (polygonal  facets)  on  geo¬ 
metric  models  for  use  in  several  prediction  software  codes  such  as  aerodynamic, 
radar  cross  section,  stress,  etc.. 

7.  NAVGRAPH  can  provide  finite  element  and  geometric  data  for  other  Air  Force 
data  base  formats  such  as  AFATL/ASE  and  VUMADM.  NAVGRAPH  can  also 
provide  geometric  model  data  (via  IGES)  to  ANVIL  for  numerically  controlled 
milling  machine  instructions.  NAVGRAPH  can  accept  finite  element  and  geome¬ 
try  data  from  other  Air  Force  data  bases  for  further  model  development  or  graphic 
display. 

8.  NAVGRAPH  can  be  used  to  generate  computer  color  graphic  displays  (line  draw¬ 
ings  or  smooth  rendering)  of  all  the  above  upon  demand  including  interference 
highlighting. 
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Part  1 

The  Technical  Description 


The  Technical  Description 

The  Phase  II  work  consisted  of  the  following  main  areas,  namely:  1)  Graphics, 
Mass  Properties  and  Hardware  Compatibility,  2)  Interference  Calculation  and  Optimal 
Packaging,  3)  Translator  Software  Development,  and  4)  Data  Base  Generation,  5)  Other 
Capabilities  and  Enhancements.  The  results  accomplished  in  each  area  are  discussed 
below. 

The  discussion  contained  herein  assumes  an  understanding  of  solid  geometric  mod¬ 
eling  as  provided  by  NAVGRAPH  and  an  understanding  of  the  optimization  techniques 
as  provided  by  OPTDES.  A  complete  discussion  of  these  two  packages  are  contained 
in  their  respective  users  manuals  that  accompany  this  report  or  can  be  obtained  from 
Brigham  Young  University  or  APTEK. 

Graphics,  Mass  Properties  and  Hardware  Compatibility 

Hardware  compatibility. 

The  complete  software  package  was  developed  to  be  compatible  to  run  on  a  Mi- 
croVajc  computer  (and  all  DEC-VMS).  All  graphics  output  is  compatible  with  Tektronix 
4010,  4109,  and  4129  graphics  terminals.  The  Tektronix  4129  terminal  driver  allows  the 
use  of  the  local  hardware  (rotate,  zoom,  pan,  etc.)  commands.  NAVGRAPH  can  ac¬ 
cept  data  input  with  the  Tektronix  4957  and  4958  digitizers  as  well  as  with  the  use  of 
thumbwheels  and  traditional  keyboard  input.  Tektronix  4692  and  4693  copiers  can  be 
used  to  get  hardcopies  of  any  display  created  by  NAVGRAPH. 

All  of  the  software  can  also  run  under  operating  systems  other  than  DEC-VMS, 
such  as  the  UNIX  operating  system.  Different  graphics  devices,  other  than  Tektronix, 
simply  need  a  different  graphics  driver  module  containing  the  graphics  machine  depen¬ 
dent  calls.  There  are  many  drivers  available,  and  new  ones  are  easily  developed. 

Mass  properties  calculator. 

The  capability  was  developed  to  calculate  the  mass  properties  (volume,  mass, 
center  of  gravity,  and  moments  of  inertia)  of  any  NAVGRAPH  solid  model. 

A  NAVGRAPH  session  is  run  to  create  the  desired  solids  and  to  group  them  to 
define  the  object  being  modeled.  An  object  is  defined  as  a  group  of  solids  that  will 
rotate  and  translate  together.  In  preparation  for  mass  properties  calculation,  each  solid 
in  the  object  is  meshed  with  linear  (i.e.  no  midside  nodes)  solid  finite  elements.  Differ¬ 
ent  linear  solid  finite  elements  (hexahedrons,  pentahedrons,  wedges,  and  tetrahedrons) 
are  generated  automatically  depending  on  the  degree  of  degeneracy  of  the  solid  being 
meshed.  If  there  are  no  degeneracies,  all  elements  will  be  hexahedrons.  If  there  is  one 
surface  degenerating  to  a  line,  there  will  be  a  row  of  wedge  elements  along  the  degenerate 
surface.  If  a  surface  degenerates  to  a  point  then  there  will  be  a  layer  of  pentahedrons. 
MASSPROP  was  written  to  handle  any  of  these  cases. 
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Once  the  object  has  ail  of  its  solids  meshed  with  solid  finite  elements,  MASSPROP 
calculates  the  mass  properties  element  by  element.  MASSPROP  uses  the  Gauss  Quadra¬ 
ture  method  of  integration  to  evaluate  the  volume  (t\),  mass  (m,)  and  centroid  (xt,  y,,  z,), 
of  each  ith  element.  Once  this  information  is  known,  the  composite  mass  properties  of 
the  whole  object  are  simply  found  with  the  proper  summation  formula  for  each  mass 
property.  The  following  equations  show  how  this  summation  is  performed  with  respect 
to  the  global  coordinate  system. 


V 

= 

(i) 

m 

= 

(2) 

X 

= 

m 

(3) 

y 

= 

Ei  m,  z, 

m 

(4) 

z 

= 

Hi  mxy, 
m 

(5) 

= 

Em.(y.2  +  *?) 

1 

(6) 

Iyy 

= 

Em.(x?  +  z2) 

$ 

(<) 

hz 

= 

+  y.2) 

(8) 

4v 

= 

-  X)  m.(x.y.) 

(9) 

/« 

= 

(10) 

4* 

= 

-£«*.-($*) 

(11) 

The  values  of  the  moments  of  inertia  with  respect  to  the  local  coordinate  system 
with  an  origin  at  (x,y,  z)  are  simply  found  with  the  parallel  axis  theorem, 


IX'X> 

=  Iix  -  rn(y2  +  z2) 

(12) 

4'jf' 

=  Iyy  ~  m(X 2  +  Z2) 

(13) 

Iz'z' 

=  Izz  ~  m(x2  +  y2) 

(14) 

IX'V' 

=  Ixv  +  m(xy) 

(15) 

1  x'z' 

=  Ixz  +  m(xz) 

(16) 

h'z' 

=  Iy,  +  m(yz ) 

(IT) 

Figure  2  shows  a  an  object  that  was  modeled,  meshed  and  had  its  mass  properties 
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Figure  2:  The  mass  properties  of  a  model  calculated  with  MASSPROP  and  compared 
with  the  closed  form  solution. 

computed  by  MASSPROP.  The  mass  properties  from  MASSPROP  are  compared  to 
those  obtained  by  closed  form  solution.  The  values  calculated  by  MASSPROP  axe 
exact  for  models  with  straight  sides,  and  are  approximate  for  models  with  curved  sides. 
As  more  solid  elements  are  used  along  a  curved  side  the  solution  approaches  the  exact 
solution.  As  a  general  rule,  10  element  subdivisions  along  a  curved  side  will  result  in  less 
than  one  percent  deviation  from  the  exact  solution.  More  complex  models  are  solved 
by  simply  creating  more  solids,  meshing  the  solids  into  finite  elements,  and  including 
the  solids  into  the  desired  object. 

Once  MASSPROP  has  calculated  the  mass  properties  for  an  object  they  are 
printed  out  to  the  screen  and  written  to  an  output  file  on  disk. 

The  mass  properties  are  then  available  to  PACKER  for  use  as  optional  constraints 
and  objective  functions  in  the  optimal  packaging  problem,  as  discussed  later  in  the 
optimal  packaging  with  mass  properties  section. 

Highlight  interferences. 

The  NAVGRAPH  code  was  enhanced  to  graphically  highlight  the  interferences 
between  objects.  This  was  accomplished  by  adding  two  commands  to  the  NAVGRAPH 
global  group  (GLOB  GROU)  menu.  While  in  the  global  group  menu,  one  needs  to 
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type  "CALF  for  CALIPER  to  detect  and  calculate  the  location  and  magnitude  of  all 
interferences.  Once  CALIPER  is  finished  with  its  calculations,  the  user  types  ’’HIGH” 
and  defines  the  desired  color  of  the  highlight,  then  a  flashing  vector  is  drawn  that 
shows  the  location,  direction  and  magnitude  of  all  interferences.  If  the  user  translates, 
rotates,  or  scales  any  group  in  the  data  base,  CALIPER  would  then  need  to  be  rerun 
for  HIGHLIGHTS  to  be  displayed  again. 

Geometric  data  base  methodologies. 

The  ability  to  use  NAVGRAPH  to  create  solid  models  of  specific  shapes,  such  as 
ellipsoids  and  tangent  ogive  bodies  of  revolution,  was  studied.  The  methodologies  used 
to  create  such  models  are  documented. 

The  solid  entity  is  a  volume  bounded  by,  at  most,  six  surface  entities  or  twelve 
edge  lines.  A  surface  entity  in  NAVGRAPH  is  a  bi-linear  blended  Coon’s  patch  from 
its  four  edge  lines,  thus  the  midspan  of  a  curved  surface  may  not  represent  the  desired 
surface  exactly.  The  degree  to  which  a  NAVGRAPH  surface  represents  a  real  surface 
depends  on  the  ability  of  the  bi-linear  interpolation  to  give  points  on  the  NAVGRAPH 
surface  that  do  not  vary  from  where  the  points  would  be  on  the  read  surface.  To  get  a 
good  representation  of  a  real  surface  with  parametric  surfaces,  one  may  need  to  provide 
more  information  by  breaking  up  the  surface  into  more  surface  segments  and  providing 
more  lines  in  the  definition.  This  gives  the  model  more  information  and  less  distance 
between  lines  in  which  to  interpolate  the  surface.  One  must  keep  this  in  mind  as  he 
creates  surfaces  from  lines  and  solids  from  surfaces  and  lines,  especially  for  bodies  of 
revolution.  The  trade-off  here  is  more  calculations  for  increased  accuracy.  That  is, 
more  separation/interference  calculations  are  performed  where  more  solids  are  needed 
to  make  complex  solid  shapes. 

An  example  is  given  here  of  the  modeling  of  a  solid  ellipsoid.  Let  us  model  an  el¬ 
lipsoid  that  is  5.0  inches  long  along  the  global  x  axis,  3.0  inch  tall  along  the  global  y  axis, 
and  7.0  inches  wide  along  the  global  z  axis.  The  lines  defining  a  surface  segment  of  the 
solid  ellipsoid  are  created  with  the  conic  line  option.  This  option  requires  three  points 
(two  end  points  and  third  point  defining  the  plane)  and  a  parameter  defining  its  relative 
altitude  (0.5  gives  a  true  parabola,  .41375  gives  a  very  close  ellipse  approximation).  The 
steps  are  listed  here. 

1.  Define  9  points  in  order  at  (5,0,0),  (5,3,0),  (0,3,0),  (0,3,7),  (0,0,7),  (5,0,7).  (5,0,0), 

(5,0,0)  and  (5,0,0) 

2.  Define  line  1  as  a  conic  with  points  1,2,3 

3.  Define  line  2  as  a  conic  with  points  3,4,5 

4.  Define  line  3  as  a  conic  with  points  5,6,1 

5.  Define  line  4  by  mirroring  line  1  about  zx  plane. 
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Figure  3:  A  parametric  solid  representation  of  an  Ellipsoid. 


6.  Define  line  5  by  mirroring  line  3  about  xy  plane. 

7.  Define  line  6  by  mirroring  line  2  about  xy  plane. 

8.  Define  line  7  by  mirroring  line  2  about  zx  plane. 

9.  Define  line  8  by  mirroring  line  6  about  zx  plane. 

10.  Define  lines  9,10,11,12  as  a  degenerate  straight  lines  between  point  1  and  point  7, 
point  7  and  point  8,  point  8  and  point  9,  point  9  and  point  1. 

11.  Define  solid  1  from  lines  6, 8, 7, 2  on  the  front  face,  lines  9,10,11,12  on  the  back  face 
and  lines  5, 4, 3,1  connecting. 

12.  Define  solid  2  by  mirroring  solid  1  about  the  yz  plane. 

13.  Put  solids  1  and  2  into  a  group  and  use  the  GROUP  ROTATE  and  GROUP 
TRANSLATE  commands  to  put  ellipsoid  in  the  desired  orientation  (see  Figure  3) 


Notice  that  the  ellipsoid  was  modeled  with  2  solids.  This  is  a  good  example  of 
using  more  solids  to  represent  the  shape  more  exactly.  The  steps  needed  to  create  a 
tangent  ogive  body  of  revolution  and  other  shapes  are  contained  in  the  tutorial. 

The  general  methodologies  studied  and  learned  as  a  result  of  the  work  done  in  the 
data  base  creation,  CALIPER  examples  and  PACKER  examples  have  been  documented 
and  supplied  in  the  users  manual/tutorial.  Basically,  however,  the  following  items  should 
be  kept  in  mind. 
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1.  Look  ahead  and  decide  how  the  desired  model  can  be  broken  up  into  piecewise 
solids  each  bound  by  12  edges. 

2.  Using  data  from  drawing,  notions,  etc.,  create  the  points  that  will  make  the  lines 
(edges)  necessary  to  make  the  solids. 

3.  Be  familiar  with  all  of  NAVGRAPH’s  command  such  as  scale,  rotate,  sweep, 
translate,  etc.. 

4.  Don’t  be  afraid  to  ’jump  in’  and  learn  by  using  NAVGRAPH. 

5.  If  something  does  not  work  the  way  it  seems  that  it  should,  try  a  ’work  around’ 
method  to  arrive  at  the  same  results. 

CALIPER:  Separation  and  Interference  Distances  for  Physical 
Fit  Study. 

A  fundamental  need  exists  for  a  robust  capability  to  calculate  the  minimum  sep¬ 
aration  distance  between  two  parametrically  defined  solid  objects  and,  if  interference  is 
detected  (e.g.  the  minimum  distance  is  less  than  or  equal  to  zero),  calculate  the  inter¬ 
ference  distance  between  the  two  solid  objects.  This  need  has  expressed  itself  in  many 
different  disciplines,  such  as  robotics  path  planing,  animation,  and  now  the  physical  fit 
compatibility  and  optimal  packaging  problems.  The  problem  is  similar  in  all  cases,  and 
has  received  some  attention  in  the  recent  past. 

Most  of  the  work  has  been  with  polyhedral  or  polygonal  geometry  where  the 
minimum  distance  is  determined  with  an  exhaustive  search  of  distances  between  the 
nodes  of  one  object  and  the  nodes  of  another,  and  between  the  nodes  of  one  object 
and  the  polygonal  sides  of  another  object,  etc..  This  leads  to  a  massive  amount  of 
data  and  computations  for  each  pair  of  objects,  especially  if  arbitrarily  shaped,  curved 
sided  objects  are  used.  Accuracy  is  limited  in  this  case  since  polyhedrons  are  flat  sided 
approximations  to  the  shape  being  modeled.  The  closer  the  objects  get  to  each  other 
the  larger  the  percent  error  in  the  distance  calculation.  With  discrete  nodal  points  and 
flat  sided  facets,  a  smooth  function  for  interference  and  separation  is  not  possible. 

Other  work  has  been  done  where  the  shapes  are  mapped  into  pixels  and  if  a  pixel 
is  occupied  by  more  that  one  object,  collision  is  detected.  The  measure  of  interference 
is  done  by  counting  pixels  of  overlap  and  a  scalar  value  of  area  for  the  2-D  problem  is 
supplied.  This  method  has  several  limitations,  1)  it  requires  a  large  amount  of  storage 
for  each  object  to  be  rasterized  into  pixels,  2)  it  has  no  vector  data  for  direction  of  the 
interference,  and  3)  it  has  no  smooth  function  of  interference  and  separation  distance. 

Another  area  that  has  received  a  lot  of  attention  is  the  area  of  computing  the 
intersections  of  two  parametrically  defined  solids,  surfaces  and/or  lines.  This  area  of 
study  is  exciting  and  provides  important  information  for  solid  geometric  modeling,  but 
the  result  of  such  methods  does  not  give  the  information  needed  for  a  physical  fit  study 
or  an  optimal  packaging  problem.  For  example,  where  two  solids  intersect  a  geometric 
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volume  or  surface  entity  is  produced,  where  two  surfaces  intersect  a  line  is  produced, 
and  where  two  lines  intersect  points  are  produced.  The  resulting  entities  do  not  give  the 
necessary  measure  of  interference  or  separation  distances  that  axe  needed  for  studying 
physical  fit  compatibility  and  performing  optimal  packaging  problems. 

All  methods  discussed  heretofore  have  similar  limitations  when  applied  to  physical 
fit  compatibility  and  optimal  packaging  problems,  that  is  they  do  not  give  a  smooth 
function  containing  the  direction  and  magnitude  of  the  separation  and/ or  interference 
between  two  solid  objects. 

A  new  and  innovative  method  to  calculate  separation  and  interference  distances 
between  parametrically  defined  complex  solids  is  presented  here.  The  algorithms  have 
been  programmed  and  incorporated  into  the  software  module  named  CALIPER. 

The  Solid  Modeling  Convention. 

A  brief  summary  of  the  solid  modeling  convention  used  is  presented  here. 

The  point  entity  is  the  basic  geometric  entity  and  is  defined  by  its  three  global 
coordinates  (x,  y,  z). 

A  line  is  defined  by  its  cubic  parametric  coefficients  and  can  represent  straight 
lines,  quadratic  parametric,  cubic  parametric,  and  spline  segments.  A  line  can  also  be 
an  arc  segment  and  a  complex  line  created  by  merging  simple  lines.  The  Cartesian 
coordinates  of  any  point  on  a  cubic  parametric  line  are  given  by  the  cubic  equations: 


3 


*(*)  = 

t'sO 

3 

(18) 

y(«)  = 

1=0 

3 

(19) 

*(*)  = 

•art) 

(20) 

where: 


3 

a' s,  Vs,ds 


Cartesian  coordinates  of  a  point, 
parametric  coordinate  between  0  and  1  along  the  line, 
(could  also  be  parameters  t  or  u). 
parametric  coefficients. 


Cartesian  coordinates  on  arcs  are  calculated  directly  from  the  closed  form  expression 
for  an  arc.  For  complex  lines,  each  line  segment  is  mapped  into  the  toted  complex  line 
according  to  its  length  so  that  the  complex  line  also  has  parametric  limits  between  0 
and  1.  With  this  mapping  it  is  possible  to  calculate  the  (x,y,r)  coordinates  for  any 
parameter  value  s  (between  0  and  1)  along  emy  line  definition. 
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Figure  4:  A  surface  is  mapped  parametrically  into  a  unit  square,  and  a  solid  into  a  unit 
cube. 


A  surface  entity  is  defined  parametrically  with  a  bilinear  blended  Coon’s  patch 
between  four  lines  connected  in  the  proper  order  and  direction.  The  lines  forming  the 
boundary  of  the  surface  can  be  any  of  the  types  above  including  the  complex  type.  The 
four  lines  are  mapped  into  a  unit  square  with  parameters  a  and  t  (see  Figure  4).  A 
point  on  the  surface  is  calculated  with  the  Coon’s  interpolation  formula.  Eight  points 
are  needed  for  the  formula  (four  midside  points  and  four  surface  corner  points). 

x(s,t)  =  (1  —  t)/i  +  s/2  +  (1  —  s)f3  +  s/4  (21) 

— x(0,0)(l  -  a)(  1  -t)-  x(l,0)s(l  -  t)  -  x(0, 1)(1  -  s)t  -  x(l,  1  )at 

(and  similar  equations  for  y(s,t)  and  z(s,t))  where: 

/,  =  x  on  line  t  (t  =  1  to  4)  as  calculated  in  (1). 

x(0,0)  =  1  at  the  a  —  0  and  t  =  0  corner. 

x(l,0)  =  x  at  the  a  =  1  and  t  =  0  comer. 

x(0, 1)  =  x  at  the  a  =  0  and  t  =  1  corner. 

x(  1 , 1 )  =  1  at  the  a  =  1  and  t  =  1  corner. 

A  solid  entity  is  defined  parametrically  in  (s,f,u)  by  bounding  the  volume  it 
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represents  with  6  surfaces  or  with  12  line  edges  mapped  into  a  unit  cube  (see  Figure  4). 
Calculation  of  any  Cartesian  coordinate  within  the  solids  parametric  limits  is  done 
in  a  similar  fashion  to  that  of  surfaces  only  with  a  tri-linear  blended  Coon’s  patch 
interpolation  between  20  points  (eight  corners  points  and  12  midside  points).  The 
equation  is  given  here: 


x(M,u)  =  (1  -  0(1  -u)/l  +  (1  -  t)uft  +  tu/a  +  t(l  -  u)f4  +  (1  -s)(l  -  u)/s 
+(1  -  s)uf$  +  suf7  +  5(1  -  u)fs  +  (1  -  s)(l  -  t)f9  +  (1  -  s)tfw 
+5t/n  +  s(l  —  t)/12  +  C(s,t,u)  (22) 


where: 


C(s,t,u)  =  — 2{(1  —  s)(l  —  t)(l  —  u)x(0, 0, 0)  +  (1  —  s)(l  —  t)ux(0, 0, 1)  + 

(1  -  s)f(l  -  u)x(0, 1,0)  +  (1  -  s)tux(0, 1, 1)  +  s(l  -  f)(l  -  u)x(l, 0, 0) 
3(1  -  t)ux(l,  0, 1)  +  st(  1  —  u)x(l,  1,0)  +  stux(l,  1,  1) 


(and  similar  equations  for  y(s,t,u)  and  z(s,t,u))  and  where: 


/, 

*(0,0,0) 

*(1,0,0) 

*(0,1,0) 

*(0,0,1) 

*(1,1,0) 

*(1,0,1) 

*(0,1,1) 

*(1,1,1) 


x  on  line  i  (i  =  1  to  12)  as  calculated  in  (1). 
x  at  the  s  =  0,  t  =  0  and  u  =  0  corner. 

x  at  the  s  =  1,  t  =  0  and  u  =  0  corner. 

x  at  the  3  =  0,  t  =  1  and  u  =  0  corner. 

x  at  the  3  =  0,  t  =  0  and  u  =  1  corner. 

x  at  the  s  =  1,  t  =  1  and  u  =  0  corner.. 

x  at  the  3  =  1,  t  =  0  and  u  =  1  corner. 

x  at  the  s  =  0,  t  =  1  and  u  =  1  corner. 

x  at  the  s  =  1,  t  =  1  and  u  =  1  corner. 


This  representation  of  solid  geometry  is  commonly  referred  to  as  the  a  boundary 
representation  of  geometric  solids.  Note  that  the  interpolation  is  a  smooth  continuous 
function  in  most  cases  defined  with  any  simple  line  types  and  at  least  continuous  in 
all  cases  defined  with  any  complex  line  types.  The  continuous  property  is  important 
for  the  numerical  method  used  in  calculating  the  minimum  separation  and  maximum 
interference  distances. 


Finally,  an  object  is  made  by  grouping  parametrically  defined  solids  together  so 
they  rotate  and  translate  as  a  single  entity.  This  is  commonly  referred  to  as  a  boolean 
add  in  Constructive  Solid  Geometry  (CSG).  This  ability  significantly  increases  the  abil¬ 
ity  to  model  very  complex,  non-convex,  arbitrarily  shaped  objects. 
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Separation  Distance  Calculation  and  Interference  Detection. 

In  the  separation  distance  calculation,  two  objects  are  defined  by  grouping  one  or 
more  parametrically  defined  solids  together  for  each  object  then  the  minimum  distance 
between  them  is  calculated.  The  given  information  is  the  geometric  information  as 
defined  above  (i.e.  the  mapping  between  parametric  space  (a,  t,  u)  and  Cartesian  space 
(x,y,  z)  for  each  solid).  The  following  optimization  problem  is  solved  for  each  solid  in 
the  first  object  with  respect  to  each  solid  in  the  second  object. 

Find: 


that  minimizes: 


D  =  ((*1  -  x2)2  +  (#1  -  J/2)2  +  (21  -  *2)2)’ 


such  that: 

0  <  ai  <  1 

0  <  U  <  1 
0  <  <  1 
0  <  s2  <  1 
0  <  t2  <  1 
0  <  u2  <  1 

where: 


$2i*2>  «2 

*ii  Vii  zi 
x2i  y2i  z2 

D 


=  parametric  coordinates  of  solid  1 
=  parametric  coordinates  of  solid  2 
=  Cartesian  coordinates  of  solid  1  from  Si,tx,ui 
=  Cartesian  coordinates  of  solid  2  from  a2)t2,u2 
=  the  distance  between  the  two  Cartesian  points 


The  optimization  problem  is  simply  the  search  for  the  two  points  that  define  the 
minimum  distance  between  them.  The  constraints  force  the  parametric  coordinates  to 
remain  in  their  respective  solids.  If  an  object  has  more  that  one  solid  then  this  process 
is  repeated  for  all  combinations  of  solids  in  object  1  with  respect  to  all  solids  in  object  2. 
The  minimum  of  all  separation  distances  for  all  combinations  is  kept  as  the  minimum 
separation  between  the  two  objects.  This  definition  is  the  same  as  that  used  in  the 
preceding  Phase  I  SBIR  study  except  for  the  grouping  of  solids  into  objects.  The  use 
of  many  solids  per  objects  significantly  increases  the  ability  to  model  complex  objects 
and  still  get  separation  distances.  If  the  minimum  distance  between  any  pair  of  solids 
is  computed  to  be  zero  or  less,  then  interference  has  been  detected.  The  calculation  of 
the  interference  distance  is  then  carried  out  as  discussed  in  the  next  section. 
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The  generalized  reduced  gradient  method  was  employed  to  solve  the  distance 
minimization  problem.  This  solution  has  been  found  to  be  very  robust  and  accurate. 
There  have  been  literally  thousands  of  test  cases  run  with  no  failures. 

Interference  Distance  Calculation 

Once  the  separation  distance  has  been  detected  to  be  less  than  or  equal  to  zero, 
the  interference  can  be  calculated.  A  new  definition  of  interference  is  presented  here. 

Given  two  objects  that  are  interfering,  interference  is  defined  as  the  distance 
needed  to  translate  the  second  object  away  from  the  first  object  along  some  vector 
until  the  two  objects  just  touch  (i.e.  incipient  separation).  This  definition  implies  both 
direction  and  magnitude.  The  direction  of  interference  may  be  defined  arbitrarily  by 
the  user  or  by  default  as  the  direction  between  the  centroids  of  the  two  objects. 

The  user  defined  direction  of  interference  provides  the  direction  and  the  magnitude 
of  the  translation  necessary  to  remove  the  interference  between  two  objects.  The  default 
direction  between  the  centroids  of  the  objects  is  best  used  in  the  optimal  packaging 
problem. 

The  interference  calculation  is  also  defined  as  an  optimization  problem.  For  any 
pair  of  solids  that  axe  detected  to  interfere  in  the  separation  calculation,  the  following 
problem  is  solved: 

Find: 


suthUi,S2,t2,U2,T 


that  maximizes: 


such  that: 


T  =  translation  of  object  2  with  respect  to  object  1 


0  <  sa  <  1 

0  <  fj  <  1 
0  <  ui  <  1 

0  <  s2  <  1 
0  <  t2  <  1 
0  <  u2  <  1 

Xj  =  x2 

y\  -  2/2 


z2 
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Where: 

sutuui  =  parametric  coordinates  of  solid  1 
s2,  ^21  u2  —  parametric  coordinates  of  solid  2 
xl,y1,z1  =  Cartesian  coordinates  of  solid  1  from  Si,ti,ui 
x2i^2iz2  =  Cartesian  coordinates  of  solid  2  from  S2,t2,u2 

Note  that  the  additional  equality  constraints  insure  that  the  second  solid  cannot 
translate  away  from  the  first  solid  any  further  than  that  amount  needed  to  reach  incipient 
separation.  The  distance  T  is  the  magnitude  of  the  translation  as  well  as  the  negative 
magnitude  of  interference.  Therefore,  interference  is  always  <  0. 

This  definition  of  interference  is  an  important  improvement  over  the  definition 
that  was  presented  in  the  Phase  I  work.  It  gives  the  same  interference  whether  solid 
1  is  translated  with  respect  to  solid  2  or  vice  versa.  It  also  handles  the  case  where 
there  is  full  penetration  of  one  object  into  another  (even  when  identical  objects  are 
spatially  coincident).  It  allows  a  continuous  function  from  separation  through  incipient 
separation  to  interference. 

There  were  some  interesting  alternative  approaches  studied  before  the  above  defi¬ 
nition  was  decided  upon.  One  method  involved  shrinking  the  two  objects  with  a  scaling 
transformation  until  incipient  separation  was  reached.  This  method  failed  in  handling 
the  full  penetration  situation.  Another  method  translated  the  second  object  far  away 
from  the  first  object  (along  the  same  translation  vector  as  described  above)  then  it¬ 
eratively  translated  the  object  back  along  the  translation  vector  (with  the  repeated 
separation  distance  calculations)  until  the  objects  just  touched.  This  method  gave  the 
same  results  as  the  chosen  method  but  required  much  more  cpu  time. 

The  generalized  reduced  gradient  method  was  also  employed  to  solve  the  trans¬ 
lation  distance  maximization  problem.  This  solution  has  also  been  shown  to  be  very 
robust  and  accurate.  Here  again,  many  test  cases  have  been  run  with  no  failures. 

CALIPER  Examples 

An  example  is  included  here  to  show  how  the  calculation  of  separation  and  inter¬ 
ference  distances  with  CALIPER  can  help  the  concept  missile  designer  solve  physical 
fit  compatibility  problems.  The  use  of  CALIPER  to  help  optimally  place  objects  into 
containers  is  discussed  in  the  next  section. 

It  is  desired  to  find  the  orientation  of  a  conceptual  bomb  on  an  F-15  wing/pylon- 
/rack  assembly  such  that  there  are  no  interferences  during  stowage,  take-off  or  landing. 
It  is  not  desirable  to  built  a  proto-type  just  for  compatibility  study.  First  the  bomb 
is  modeled  with  NAVGRAPH.  Next,  DBMERGE  is  used  to  combine  the  NAVGRAPH 
data  base  files  containing  models  of  the  F-15  wing  with  swept’  gear,  the  runway,  the 
pylons,  the  rack,  and  the  conceptual  bomb  (simulated  here  with  an  MK84)  into  one 
NAVGRAPH  file  (see  Figure  5).  A  NAVGRAPH  session  is  run  and  the  CALIPER 
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Figure  5:  CALIPER  calculates  the  separation  and  interference  distances  between 
weapon  and  wing/pylon/landing  strip  assembly. 

command  issued.  Giving  the  two  object  names  (the  bomb  and  the  wing/pylon/rack 
assembly  in  this  case)  and  the  direction  of  translation  for  the  interference  calculation, 
CALIPER  returns  the  distance  needed  to  move  the  bomb  such  that  there  is  no  more 
interference.  Figure  5  shows  the  distance  and  location  of  the  interference  (a  flashing 
vector  is  actually  shown  on  the  color  monitor).  The  user  simply  translates  the  bomb 
along  the  direction  defined  and  with  the  magnitude  calculated  to  remove  the  interfer¬ 
ence.  Then  CALIPER  can  be  run  again  to  make  sure  there  are  no  other  interferences 
caused  by  the  bomb’s  recent  translation.  The  process  is  repeated  for  each  weapon  and 
piece  of  equipment  to  be  stowed  on  the  aircraft. 

This  process  can  also  be  thought  of  as  representing  the  volume  available  for  stores 
on  an  aircraft  by  showing  the  ’volume  of  exception’  available  to  the  stores  placement 
design  engineer. 

The  accuracy  obtained  by  CALIPER  in  predicting  the  separation  and  interference 
distances  is  within  ss  .001  units  for  the  worst  case  observed,  however,  most  cases  are 
much  better.  This  is  well  within  the  limits  set  by  the  Air  Force  for  this  development, 
that  is,  within  |  inch  for  submunition  to  submunition  distances  and  J  inch  for  stores  to 
aircraft  distances. 
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PACKER  with  rotations  and  translation. 

A  fundamental  need  exists  for  the  ability  to  automatically  package  objects  into 
containers.  The  optimal  packaging  problem  arises  constantly  in  our  lives.  Anyone  who 
has  ever  packed  a  trunk  or  suitcase  understands  the  difficulty  of  placing  as  many  objects 
as  possible  into  a  limited  amount  of  space.  Packaging  is  generally  performed  by  trial 
and  error  which  is  time  consuming  and  inefficient.  In  industry  packaging  problems  are 
becoming  increasingly  important  as  companies  seek  to  minimize  their  costs  and  increase 
efficiency.  These  problems  are  diverse  and  range  from  part  layout  in  the  textiles  and 
sheet-metal  industries  to  packaging  freight  in  the  transportation  industry.  The  Air  Force 
sees  this  problem  arise  when  packaging  interned  components  and  submunitions  into 
missile  cargo  bays.  Consequently,  automation  of  the  packaging  problem  has  received  a 
great  deal  of  interest.  The  latest  related  work  in  this  field  is  reviewed  here. 

Two-dimensional  problems  have  received  considerable  attention  from  researchers. 
These  generally  involve  the  placement  of  planar  shaped  objects  inside  some  bounding 
area  with  the  objective  of  minimizing  the  distances  between  objects.  For  instance  in  the 
textile  industry  this  would  apply  to  cutting  material  so  as  to  waste  as  little  as  possible 
Another  application  is  the  design  of  integrated  circuits  where  individual  components  are 
arranged  to  produce  the  smallest  design  possible.  In  general,  two  approaches  have  been 
used  in  solving  the  optimal  packaging:  numerical  optimization  and  heuristic  algorithms. 

Most  all  work  done  in  this  area  has  been  done  in  the  two-dimensional  layout 
problems  where  the  objects  packaged  are  represented  with  polygons  which  can  have  an 
arbitrary  number  of  edges  including  non-convex  shapes.  The  design  variables  axe  the  x 
and  y  translations  and  the  one  rotation  of  each  polygon.  The  objective  was  usually  to 
minimize  the  area  needed  to  bound  the  polygons  without  overlap. 

Very  little  was  found  in  the  literature  on  work  done  in  the  three-dimensional  realm 
of  optimal  packaging.  One  project  has  addressed  whether  a  given  shape  would  fit  in  a 
given  size  box  or  and  to  find  the  minimum  size  box  that  would  fit  around  a  given  object. 
The  work  was  limited  to  one  object  in  a  container. 

The  Phase  I  work  preceding  this  Phase  II  work  showed  the  feasibility  of  using 
optimization  techniques  and  solid  geometric  modeling  software  to  package  3-D  convex 
shapes  into  3-D  containers.  Even  though  feasibility  was  shown,  some  limitations  were 
experienced.  One  fundamental  problem  was  with  the  use  of  Euler  angles  for  rotational 
degrees  of  freedom.  It  was  found  that  the  optimization  routines  were  often  unable  to 
rotate  the  objects  and  the  process  would  stop  prematurely.  Another  problem  was  that 
an  object  could  only  be  represented  with  one  arbitrary  solid  which  limited  the  freedom 
to  model  very  complex  shapes. 

The  Phase  II  work  has  done  much  to  improve  the  robustness  and  functionality 
of  PACKER  while  at  the  same  time  keeping  the  necessary  computing  resources  at  a 
reasonable  limit.  The  state  of  the  software  is  detailed  here  by  describing  how  the  many 
tasks  completed  during  Phase  II  have  been  incorporated  into  the  software  program 
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PACKER.  We  first  start  by  outlining  the  properties  of  PACKER  and  then  discuss  the 
steps  taken  to  optimally  package  objects  into  container  with  spatial  relationships.  The 
next  main  section  will  add  the  use  of  mass  properties. 

Geometric  modeling  convention 

The  method  of  representing  3-D  arbitrarily  shaped  solid  objects  is  the  same  as  the 
geometry  modeling  convention  as  used  for  CALIPER  described  in  the  previous  main 
section.  The  ability  to  group  arbitrarily  shaped  solids  together  to  form  objects  of  very 
complex,  non-convex,  arbitrary  shapes  is  one  of  the  most  important  improvements  over 
the  Phase  I  work.  Essentially,  the  user  now  can  model  any  shape  needed  and  is  only 
limited  by  his  own  imagination  and  ability  to  use  the  modeling  software. 

Rotational  degrees  of  freedom  and  the  quaternion 

The  question  of  how  to  rotate  objects  properly  in  the  optimization  process  is  an 
important  one.  This  next  section  discusses  the  work  done  on  the  general  3-D  packaging 
problem  addressing  the  rotational  degrees  of  freedom. 

A  fundamental  component  of  the  solution  to  the  packaging  problem  is  the  ability 
to  change  the  position  and  orientation  of  the  objects  being  packaged.  It  was  proven 
by  Euler  that  the  position  of  a  rigid  body  may  be  described  by  a  translation  and  one 
rotation  about  am  arbitrary  axis.  However,  as  shown  in  Appendix  I,  Euler  derived  the 
rotation  formula  in  terms  of  three  independent  rotations.  These  rotations  know  as  Euler 
angles  are  unsuitable  for  use  in  numerical  optimization. 

Euler  angles  suffer  from  singularities.  The  most  severe  is  known  as  “gimbal  lock”, 
where  in  certain  orientations,  a  rotational  degree  of  freedom  is  lost.  Euler  angles  suffer 
from  gimbal  lock  because  they  ignore  the  interdependence  of  the  axes  of  rotation.  When 
applying  a  series  of  translations  to  an  object  the  final  position  of  the  object  may  be 
described  by  one  translation.  Clearly  a  series  of  translations  along  one  axis  of  the 
coordinate  system  are  totally  independent  of  any  translation  along  any  of  the  other  two 
axes.  Rotations  are  not  like  translations;  translations  add  while  rotations  are  multiplied 
and  the  axes  of  a  series  of  rotations  are  not  independent.  For  example,  the  x  axis  is  not 
a  linear  combination  of  the  y  and  z  axis.  The  fact  that  rotation  axes  multiply  with  a 
cross  product  is  called  the  confounding  of  axes.  Euler  Jingles  ignore  this  confounding 
of  axes  when  trying  to  achieve  an  orientation  by  rotation  about  three  independent 
axes.  Ignoring  this  cross-product  interdependence  of  rotation  axes  causes  gimbal  lock 
by  aligning  two  of  the  three  axes. 

Singularities  are  not  unique  to  Euler  angles.  It  has  been  shown  that  any  three 
independent  parameters  can  not  describe  orientation  and  be  both  global  and  nonsin¬ 
gular.  A  minimum  of  four  variables  are  needed  to  uniformly  describe  orientation  and 
be  nonsingular.  So  normally  when  using  Euler  angles  a  general  rotation  matrix  is  used 
which  has  the  singularity  at  an  orientation  which  hopefully  is  seldom  needed. 
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These  problems  with  Euler  angles  make  them  poorly  suited  for  use  with  numerical 
optimization.  Another  approach  to  describing  rotations  is  the  use  of  quaternion  calculus. 
Quaternions  have  been  around  for  about  150  years  and  were  virtually  forgotten  until 
they  came  to  be  used  in  the  spacecraft  dynamics  and  in  robotics  for  trajectory  planning. 

Quaternions  solve  the  problems  associated  with  Euler  angles  in  representing  orien¬ 
tations.  Quaternions  can  achieve  any  orientation  as  a  single  rotation  around  an  arbitrary 
axis  so  that  rotations  are  independent  for  any  orientation.  The  quaternion  describes  the 
object’s  new  orientation  from  its  original  position.  If  another  change  in  orientation  is 
desired  from  the  intermediate  position,  a  quaternion  may  again  be  applied  to  the  object 
with  the  final  quaternion  being  given  by  the  product  of  the  first  and  second  quaternions. 
Quaternions  are  able  to  do  this  because  the  cross  product  is  preserved  in  the  quaternion 
product  just  as  they  are  in  rotations. 

The  optimal  packaging  problem 

The  packaging  problem  is  posed  as  an  optimization  problem  in  the  form: 

Find:  rotations  and  translations  for  each  object 

Minimize  or  Maximize:  Some  user  defined  objective  function 

Subject  to:  Constraints  that  no  object  interfere  with  each  other 

and  that  no  object  interferes  with  container  walls. 

The  design  variables  in  the  packaging  problem  are  the  three  rotational  and  the 
three  translational  degrees  of  freedom  per  object.  All  design  variables  are  applied  to 
the  object  with  respect  to  its  original  position.  Objects  are  rotated  via  the  quaternion 
about  the  quaternion’s  direction  vector  passing  through  the  objects  parametric  center. 
This  in  effect  rotates  the  object  0  degrees,  about  an  axis  through  the  object’s  local 
centroid. 

Quaternions  were  applied  to  the  packaging  problem  by  setting  the  components  of 
the  axis  vector  (nlt  n2,n 3)  as  design  variables.  A  quaternion  has  four  degrees  of  freedom 
(0,  nl,  n2,  n3)  ,  so  in  order  to  eliminate  one  degree  of  freedom,  the  angle  of  rotation  0  was 
defined  as  the  magnitude  of  the  axis  vector  in  radians.  The  rotation  angle  is  therefore 
restricted  to  always  being  a  positive  quantity,  however,  this  is  not  a  problem.  With 
four  variables,  a  quaternion  is  able  to  represent  all  possible  orientations  twice,  since 
a  positive  angle  and  a  positive  axis  vector  (0,ni,n2,n3)  are  equivalent  to  a  negative 
angle  and  a  negative  vector  (—0,  — nt,  —  n2,  —  n3).  By  constraining  the  angle  to  remain 
positive  the  quaternion  is  still  able  to  represent  any  possible  orientation.  The  other 
design  variables  are  the  x ,y,z  translations  for  each  object. 

The  constraint  functions  are  the  distances  between  objects  and  the  distances  be¬ 
tween  objects  and  container  walls,  which  may  be  either  a  separation  or  an  interference. 
If  the  objects  interfere,  the  distance  is  negative,  and  if  the  objects  are  separated,  the 
distance  is  positive.  Therefore,  in  the  packaging  problem  all  distances  between  objects 
and  between  objects  and  container  walls  ire  constrained  to  be  positive  so  that  there  are 
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no  interferences. 


The  objective  function,  in  many  cases,  is  to  minimize  or  maximize  the  separa¬ 
tion  distances  between  objects  and  a  container  wall,  or  the  distances  between  the  ob¬ 
ject  centroids  and  a  container  wall.  The  separation  distance  is  defined  as  described  in 
CALIPER.  If  several  objects  are  to  be  packaged  to  the  bottom  of  a  container  the  objec¬ 
tive  function  would  be  computed  as  the  sum  of  the  separation  distances  of  the  objects 
with  the  bottom  container  surface  of  the  container.  The  optimization  procedure  then 
minimizes  the  sum  of  the  separation  distances  by  translating  and  rotating  the  objects. 

The  following  steps  are  taken  to  solve  a  packaging  problem. 

Step  1:  Create  and  initializing  a  starting  design. 

The  user  creates  the  initial  design  in  NAVGRAPH  by  modeling  the  objects  and 
the  container,  and  then  by  translating  and  rotating  the  objects  in  NAVGRAPH  into 
their  intended  initial  position.  In  fact,  one  might  try  to  solve  the  packaging  problem 
without  the  aid  of  the  optimization  routines  simply  by  trial  and  error  within  NAV¬ 
GRAPH.  This  can  be  accomplished  by  moving  the  objects  until  they  appear  to  be  in 
an  optimal  position.  The  CALIPER  command  in  the  NAVGRAPH  GROUP  menu  is 
the  interference  detector.  If  one  did  not  have  the  interference  detection  capability,  he 
would  need  to  try  to  detect  the  interference  visually  by  viewing  the  model  from  various 
positions.  This  is  virtually  impossible  for  a  model  of  any  real  size  with  a  large  number 
of  objects. 

The  most  efficient  solution  to  the  packaging  problem  is  to  have  the  user  and 
PACKER  complement  each  other.  The  major  advantages  of  PACKER  is  that  it  guar¬ 
antees  that  interferences  do  not  occur:  an  extremely  difficult  task  for  humans.  One 
disadvantage  of  numerical  optimization  is  that  the  algorithms  may  get  trapped  in  lo¬ 
cal  minima  and  terminate  the  search  before  the  global  minimum  is  reached,  whereas 
a  human  can  provide  heuristics  and  avoid  most  local  minimum.  A  simple  example  of 
how  optimization  algorithms  can  get  caught  in  local  minima  is  shown  in  Figure  6.  The 
objective  function  of  this  problem  was  to  minimize  the  distance  between  the  centroid  of 
the  rectangle  and  the  bottom  of  the  container.  The  constraints  are  that  the  rectangle 
is  not  allowed  to  penetrate  the  sides  of  the  container.  It  is  obvious  that  the  rectangle 
may  be  rotated  90  degrees  and  translated  down  to  the  bottom  of  the  container  but 
the  optimization  routines  are  unable  to  find  this  design  because  the  objective  function 
would  initially  increase  or  constraint  functions  would  be  violated  with  any  rotation  or 
translation. 

Even  a  simple  packaging  problem  may  contain  many  local  minima.  Therefore,  it 
is  apparent  that  the  packaging  routines  become  much  more  efficient  if  they  are  given  a 
good  starting  design.  The  better  the  starting  design,  the  better  the  final  design  and  the 
quicker  the  optimal  design  is  found.  It  is  also  apparent  that  if  the  optimization  problem 
is  guided  by  expertise  and  common  sense,  many  of  the  local  minima  will  be  avoided. 
This  means  that  even  though  this  tool  is  powerful  and  useful,  if  the  user  can  couple 


Figure  6:  Example  of  local  optima  in  a  packaging  problem. 

his  expertise,  common  sense,  and  creativity  with  the  program,  he  may  solve  even  more 
complex  problems  than  either  program  or  user  could  solve  separately. 

Once  the  model  has  been  created,  it  needs  to  be  initialized  by  running  a  version  of 
CALIPER  called  PREPAC.  PREPAC  simply  calculates  all  separation  distances  between 
all  objects  and  between  all  objects  and  container  walk,  then  writes  an  OPTDES  problem 
file  containing  all  initial  analysis  variables  and  analysis  functions  for  input  into  SETUP 
and  eventually  PACKER. 

Step  2:  The  Setup. 

To  setup  the  packaging  problem,  one  executes  the  SETUP  module  of  OPTDES 
that  reads  in  the  problem  file  from  PREPAC  then  defines  design  variables,  constraint 
functions,  and  the  objective  functions  for  the  desired  optimization  problem. 

The  design  variables  are  defined  as  a  subset  of  the  analysis  variables.  If  design 
variables  are  not  defined,  the  object  to  which  the  analysis  variables  belong  will  not 
move  in  the  direction  of  the  design  variables  omitted.  The  user  can  use  this  to  his 
advantage  in  problem  setup  to  keep  objects  from  moving  with  respect  to  certain  degrees 
of  freedom  (e.g.  let  objects  translate  but  not  rotate,  etc.)  Upper  and  lower  bounds  are 
placed  on  the  variables.  These  bounds  are  important  as  they  greatly  effect  the  scaling 
of  the  problem.  By  varying  the  bounds  on  the  rotations  and  translations  it  is  possible 
to  favor  one  over  the  other.  This  becomes  important  in  highly  constrained  problems 


where  there  are  tight  fits  between  objects  and  container  wall. 

Constraint  functions  are  defined  as  inequality  functions  by  restricting  the  functions 
to  be  either  greater,  or  less  than  some  value.  In  the  packaging  problem,  the  separation 
distances  were  constrained  to  be  greater  than  some  very  small  positive  value.  This  is 
because  it  was  found  that  by  putting  a  small  buffer  zone  around  the  objects,  (of  say, 
less  that  .01  inches)  less  interference  calculations  were  needed  near  the  optimal  solution. 
This  practice  saves  execution  time. 

The  objective  function  is  chosen,  and  then  minimized  or  maximized  depending  on 
the  desired  optimization  problem.  More  than  one  function  may  be  mapped  to  one  ob¬ 
jective  function  so  the  objective  function  may  be  defined  as  the  sum  of  all  the  separation 
distances  between  all  objects  and  a  certain  container  wall.  It  has  been  discovered  that  if 
only  spatial  relationships  are  available  for  design  functions  and  if  one  wishes  to  package 
the  objects  tighfly  to  one  end  (or  side)  of  the  container,  it  is  much  better  to  maximize 
the  sum  of  all  separation  distances  from  an  opposite  wall  rather  than  minimize  the  sum 
of  all  separation  distances  with  respect  to  the  side  desired.  This  is  because  during  a 
minimization  problem,  the  objective  function  is  not  as  sensitive  to  changes  in  the  sepa¬ 
ration  distances  (when  objects  rotate  or  translate)  as  it  is  when  the  objective  function 
is  maximized.  The  mass  properties  as  design  functions  (discussed  later)  behave  much 
better  for  this  type  of  packaging. 

Step  3:  Optimal  search  for  packaging  problem. 

PACKER  is  executed  to  search  for  the  optimal  placement  of  all  objects  in  the 
container.  This  is  where  the  problem  is  actually  solved.  The  user  reads  in  the  problem 
file  created  by  SETUP  and  geometry  file  created  in  NAVGRAPH.  This  step  is  basically 
the  DESIGN  module  from  the  OPTDES  package  with  the  analysis  function  call  being 
CALIPER.  An  optimization  algorithm  is  chosen  and  the  optimization  software  begins  a 
search  for  the  optimal  design.  The  generalized  reduced  gradient  has  been  found  to  be  the 
best  since  it  handles  highly  constrained  problems  with  the  best  results.  The  algorithm 
runs  a  specified  number  of  iterations  or  until  an  optimum  is  reached  or  assumed. 

While  in  PACKER  the  user  may  also  perform  trial  and  error  packaging  with  the 
SET  UDV  (set  unsealed  design  variables)  command.  Here  the  rotations  and  transla¬ 
tions  may  be  modified  for  each  object  and  PACKER  will  automatically  inform  the  user 
if  constraints  were  violated  by  the  operation.  Even  if  constraints  are  violated,  the  gener¬ 
alized  reduced  gradient  algorithm  can  find  feasible  designs  from  a  non-feasible  starting 
position. 

Step  4:  Post  processing. 

At  any  iteration  during  the  search  for  the  optimal  solution,  the  user  may  view  the 
current  state  of  the  problem.  The  PRO  POS  command  while  in  OPTDES  puts  the  user 
into  a  NAVGRAPH  menu  where  a  DRAW  command  may  be  issued.  All  commands  in 
the  GLOBAL  menu  are  available  to  the  user  (except  those  that  modified  the  location 
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and  orientation  of  the  objects).  At  any  given  iteration,  the  analysis  variables  are  used 
to  update  the  data  base  so  that  the  correct  view  is  drawn.  When  finished  drawing,  the 
user  may  exit  or  return  to  the  OPTDES  menu  for  further  packaging  work.  If  the  user 
exits  the  NAVGRAPH  menu,  he  can  save  the  current  design  for  a  starting  design  of 
another  optimization  process. 

Step  5:  Further  optimization. 

The  user  can  exit  PACKER  at  any  time  and  all  design  iterations  are  stored  in  a 
history  file  for  later  retrieval  or  continuation  of  an  optimization  process.  Actually,  the 
user  has  all  options  available  that  are  detailed  in  the  OPTDES  users  manual  supplied 
with  this  report. 

Even  though  the  optimal  packaging  problem  presented  here  in  this  section  is  useful, 
it  is  limited  to -packaging  with  the  spatial  relationships  provided  by  the  separation 
and/or  interference  distances  from  CALIPER.  Many  of  the  packaging  problems  would 
stop  prematurely  in  local  minima.  This  is  mainly  due  to  the  multiple  functions  in  the 
objective  function.  To  remedy  this,  there  needs  to  be  analysis  functions  for  optimal 
packaging  that  do  not  require  multiply  defined  objective  functions,  such  as,  the  mass 
properties  as  discussed  in  the  next  section. 


PACKER  with  spatial  relationships  and  mass  properties. 

In  the  previous  section  optimal  packaging  with  just  spatial  relationships  was  dis¬ 
cussed.  This  section  discusses  the  optimal  packaging  software  as  enhanced  to  include 
the  mass  properties  as  possible  design  functions.  There  was  no  literature  found  in  the 
engineering  and  computer  graphics  journals  showing  any  related  work  that  has  been 
done  in  the  area  of  optimally  packaging  objects  into  containers  with  mass  properties  as 
design  functions.  This  section  also  discusses  how  to  setup  and  run  packaging  problems 
with  mass  properties  and  compares  this  to  packaging  with  just  spatial  relationships. 

Much  of  the  general  discussion  from  the  previous  section  applies  here.  The  use 
of  mass  properties  gives  additional  capabilities  and  options  to  the  user  for  packaging 
objects  into  containers. 

Geometric  modeling  convention 

The  method  of  representing  3-D  arbitrarily  shaped  solid  objects  is  the  same  as 
the  geometry  modeling  convention  which  was  used  for  CALIPER  and  PACKER,  with 
spatial  relationships  as  described  in  the  previous  two  sections,  with  one  additional  step. 
Each  object  needs  its  mass  properties  calculated  with  MASSPROP  in  its  initial  design 
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configuration. 

Rotational  degrees  of  freedom  and  the  quaternion 


The  design  variables  remain  the  same,  that  is  the  translations  and  rotations  (via 
the  quaternion)  for  each  object  are  defined  as  before.  The  discussion  on  rotations  with 
quaternions  is  found  in  the  previous  section. 

Mass  properties  for  the  packaging  problem 

The  mass  properties  for  each  object  and  container  in  the  initial  design  configu¬ 
ration  are  stored  in  a  mass  properties  input  file  named  massprop.MAS.  The  packaging 
software  needs  each  object’s  v,p,x,  y,  z,  Jvv,  I22, 1^,  IX2,  and  I2y  (where  p  =  density). 
The  format  of  the  massprop.MAS  file  is  discussed  in  the  users  manual.  Upon  running 
PREP  AC,  the  design’s  composite  mass  properties  are  calculated  with  the  appropriate 
summation  formula  (as  in  equations  1  through  11  above  except  that  the  index  i  here 
represents  the  ith  object  rather  than  the  ith  element).  The  additional  mass  proper¬ 
ties  functions  are  then  written  to  the  OPTDES  problem  file  for  use  in  SETUP  and 
PACKER. 

As  objects  are  translated  and  rotated  during  the  optimal  design  search,  the  design 
functions  change.  The  separation  distances  change  as  objects  are  moved  (i.e.  the 
design  variables  are  changed)  by  OPTDES  or  the  user.  When  an  object  is  moved 
the  design  functions  are  recalculated.  The  separation  distances  are  recalculated  with 
another  CALIPER  call.  The  mass  properties  are  updated  by  recalculation  as  follows. 

The  volume  and  mass  remain  constant  because  objects  are  not  added  and  the 
container  and  objects  do  not  change  shape  during  am  optimal  packaging  session. 

The  x,  y  and  z  are  recalculated  with  the  following  equations: 


T;  TTliXi 
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Rotations  do  not  affect  x,  y  and  z  because  the  quaternion  rotates  an  object  about  a 
vector  through  the  object’s  centroid. 

The  moments  of  inertia  and  products  of  inertia  define  a  tensor  (a  3X3  moment 
of  inertia  matrix).  Therefore,  the  inertia  terms  are  interdependent  with  respect  to 
rotations.  A  rotation  transformation  matrix  is  needed  to  update  each  object’s  inertia 


tensor  in  the  object’s  local  coordinate  system  with  respect  to  the  global  coordinate 
system.  Each  object’s  local  coordinate  system  is  initially  parallel  to  the  global  coordinate 
system  with  its  origin  at  the  object’s  centroid.  But  as  objects  axe  rotated  each  object’s 
inertia  tensor  with  respect  to  the  global  coordinate  system  is  changed.  The  rotations 
are  in  terms  of  a  quaternion,  therefore,  we  need  a  rotation  transformation  matrix  from 
the  quaternion  design  variables.  The  requisite  transformation  matrix  R  is  given  as: 

R  = 

2(eJ  +  ei2)  —  1  2(eie2  —  eoe3)  2(eie3  +  eoe2) 

2(e!e2  +  eoe3)  2(e§  +  e22)  -  1  2(6363-606!) 

.  2(eie3  —  eoe2)  2(e2e3  +  eoej)  2(e£  +  e32)  —  1 

where 

Co 

ej 
C3 

and 

ej  +  e\  +  c|  +  e|  =  1. 

The  moment  of  inertia  matrix  for  each  object  in  its  new  orientation  is  given  with: 

I  =  RIR-1. 

The  global  inertia  terms  for  the  composite  model  are  the  sum  of  the  individual  trans¬ 
formed  inertia  terms  via  the  parallel  axis  theorem. 

The  optimal  packaging  problem 

The  packaging  problem  is  posed  again  as  an  optimization  problem  in  the  form: 

Find:  rotations  and  translations  for  each  object 

Minimize  or  Maximize:  Some  user  defined  objective  function 

Subject  to:  Constraints  that  no  object  interfere  with  each  other 

and  that  no  object  interferes  with  container  walls. 

Also  other  optional  constraints  on  mass  properties. 

The  design  variables  in  the  packaging  problem  are  still  the  three  rotational  and 
the  three  translational  degrees  of  freedom  per  object.  All  design  variables  are  applied  to 
the  object  with  respect  to  its  original  position.  Objects  are  rotated  via  the  quaternion 
about  the  quaternion’s  direction  vector  passing  through  the  object’s  centroid. 


=  cos  | 

=  ni  sin  | 
=  n2  sin  | 
=  n3  sin  | 
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The  spatial  constraint  functions  are  still  that  the  distance  between  objects  and 
the  distances  between  objects  and  container  walls  remain  >  0.  Additionally,  the  mass 
properties  allow  more  options  to  constrain  the  design  to  simulate  reed  problems.  For 
example,  one  may  set  equality  constraints  on  the  location  of  the  center  of  gravity  coor¬ 
dinates,  or  constrain  the  roll  moment  of  inertia  to  be  less  than  some  value,  etc.. 

The  objective  function  now  can  be  defined  to  minimize  y  to  package  objects  to 
the  bottom  of  a  container.  This  has  proven  to  be  a  much  better  definition  since  there 
is  not  the  multiple  function  to  one  objective  function  mapping.  Other  possible  objec¬ 
tive  functions  are:  1)  minimize  or  maximize  x  or  z  or  some  combination,  2)  minimize 
roll,  pitch  or  yaw  moments  of  inertia  (this  tends  to  package  the  object  toward  the  cen¬ 
troid),  3)  minimize  the  products  of  inertia  (this  tends  to  place  the  objects  in  a  fashion 
that  approaches  symmetry)  and  4)  some  combination  of  the  above  with  care  to  apply 
proper  coefficients  in  the  multiple  objective  function  mapping  for  proper  scaling  in  the 
optimization. 

The  following  steps  are  taken  to  solve  a  packaging  problem  with  mass  properties. 
Step  1:  Create  and  initialize  a  starting  design. 

The  user  creates  the  initial  design  in  NAVGRAPH  by  modeling  the  objects  and 
the  container,  and  then  by  translating  and  rotating  the  objects  in  NAVGRAPH  into 
their  intended  initial  position. 

Step  2:  Calculate  mass  properties. 

Each  object  needs  its  mass  properties  supplied.  This  can  be  done  with  several 
methods.  One  way  is  to  mesh  each  object  with  solid  finite  elements  and  run  MASSPROP 
for  each  object.  Another  way  is  to  use  existing  data  if  measured  experimentally,  or 
previously  with  MASSPROP.  Each  object’s  mass  properties  need  to  be  input  into  the 
file  massprop.MAS.  The  massprop.MAS  file  is  read  by  PREPAC  and  PACKER  at  a 
user  prompt. 

Once  the  model  has  been  created,  it  needs  to  be  initialized  by  running  PREPAC. 
PREPAC  calculates  all  separation  distances  between  all  objects  and  between  all  objects 
and  container  walls,  and  then  initializes  the  values  of  total  volume,  total  mass,  composite 
x,  y,  z,  and  the  moments  of  inertia  about  the  local  coordinate  system  defined  by  a 
parallel  axis  system  with  its  origin  at  (x,  y,  z).  PREPAC  then  writes  an  OPTDES 
problem  file  containing  all  initial  analysis  variables  and  analysis  functions  for  input  into 
SETUP  and  eventually  PACKER. 

Step  3:  The  Setup. 

To  setup  the  packaging  problem,  one  executes  the  SETUP  module  of  OPTDES 
that  reads  in  the  problem  file  from  PREPAC  then  defines  design  variables,  constraint 
functions,  and  the  objective  functions  for  the  desired  optimization  problem. 
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The  design  variables  are  defined  as  a  subset  of  the  analysis  variables.  If  design 
variables  are  not  defined,  the  object  to  which  the  analysis  variables  belong  will  rot  move 
in  the  direction  of  the  design  variables  omitted.  The  user  can  use  this  to  his  advantage  in 
problem  setup  to  keep  objects  from  moving  with  respect  to  certain  degrees  of  freedom 
(e.g.  let  objects  translate  but  not  rotate,  etc.).  Upper  and  lower  bounds  are  placed 
on  the  variables.  These  bounds  are  important  as  they  greatly  effect  the  scaling  of  the 
problem.  By  varying  the  bounds  on  the  rotations  and  translations  it  is  possible  to  favor 
one  over  the  other.  This  becomes  important  in  highly  constrained  problems  where  there 
are  tight  fits  between  objects  and  the  container  wall. 

Again,  the  separation  distances  should  be  constrained  to  be  greater  than  some 
very  small  positive  value  to  save  interference  distance  calculations  when  the  solution 
is  near  the  optimal  design.  If  packaging  material  is  needed  around  the  objects,  the 
constraint  limits  should  reflect  the  necessary  thickness  of  the  material. 

The  objective  function  is  chosen  here  from  the  design  functions  defined  by  the  mass 
properties.  It  has  been  discovered  that  if  one  wishes  to  package  the  objects  as  tightly 
to  one  end  (or  side)  of  the  container,  it  is  much  better  to  maximize  or  minimize  the 
coordinates  of  the  global  composite  centroid  (z,  y,  z,  or  some  combination)  instead  of  the 
sum  of  all  separation  distances  from  a  container  wall.  This  is  because  the  optimization 
function  is  at  most  a  mapping  of  only  several  mass  properties  functions  rather  than  a 
multitude  of  separation  distances.  An  objective  function  defined  by  minimization  or 
maximization  of  x ,  for  example,  is  much  more  sensitive  to  translations  and  rotations  of 
objects  than  any  objective  function  defined  with  spatial  relationships  alone. 

Step  4:  Optimal  search  for  packaging  problem. 

PACKER  is  executed  to  search  for  the  optimal  placement  of  all  objects  in  the  con¬ 
tainer.  This  is  where  the  problem  is  actually  solved.  The  user  reads  in  the  problem  file 
created  by  SETUP,  the  geometry  file  created  in  NAVGRAPH,  and  the  massprop.MAS 
file.  This  step  is  basically  the  DESIGN  module  from  the  OPTDES  package  with  the 
analysis  function  module  being  CALIPER.  An  optimization  algorithm  is  chosen  and  the 
optimization  software  begins  a  search  for  the  optimal  design.  The  generalized  reduced 
gradient  has  been  found  to  be  the  best  since  it  handles  highly  constrained  problems 
with  the  best  results.  The  algorithm  runs  a  specified  number  of  iterations  or  until  an 
optimum  is  reached  or  assumed. 

While  in  PACKER  the  user  may  also  perform  trial  an  error  packaging  with  the 
SET  UDV  (set  unsealed  design  variables)  command.  Here  the  rotations  and  translations 
may  be  modified  for  each  object  and  PACKER  will  automatically  inform  the  user  if 
constraints  were  violated  by  the  operation.  The  DIS  DV  and  DIS  DF  command  can  be 
used  to  view  the  values  of  the  design  variables  and  design  functions  at  any  iteration. 
Even  if  constraints  are  violated,  the  generalized  reduced  gradient  algorithm  can  find 
feasible  designs  from  a  non-feasible  starting  position. 

Step  5:  Post  processing. 
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Again,  at  any  iteration  during  the  search  for  the  optimal  solution,  the  user  may 
view  the  current  state  of  the  problem.  The  PRO  POS  command  while  in  OPTDES 
puts  the  user  into  a  NAVGRAPH  menu  where  a  DRAW  command  may  be  issued.  All 
commands  in  the  GLOBAL  menu  are  present  to  the  user  (except  those  that  modify  the 
location  and  orientation  of  the  objects).  At  any  given  iteration,  the  analysis  variables 
are  used  to  update  the  data  base  so  that  the  correct  view  is  drawn.  When  finished 
drawing,  the  user  may  exit  or  return  to  the  OPTDES  menu  for  further  packaging  work. 

Step  6:  Further  optimization. 

The  user  can  exit  PACKER  at  any  time  and  all  design  iterations  are  stored  in  a 
history  file  for  later  retrieval  or  continuation  of  an  optimization  process.  Actually,  the 
user  has  all  options  available  that  are  detailed  in  the  OPTDES  users  manual  supplied 
with  this  report. 

The  mass  properties  give  much  improvement  over  just  spatial  relationships  in  the 
packaging  process.  They  help  the  search  routines  search  for  the  optimal  solution  faster 
by  supplying  a  function  that  does  not  need  multiple  mapping  in  the  objective  function. 
It  also  gives  the  user  much  more  latitude  in  objective  function  definition. 

Examples  of  packaging  problems  solved  with  PACKER 

Several  examples  are  presented  here  showing  optimal  packaging  problems  that 
were  solved  with  PACKER. 

The  first  example,  shown  in  Figure  7,  has  a  container  with  a  flat  bottom  and 
straight  sides.  Three  cylinders  are  placed  in  a  feasible  starting  design  where  there  are  no 
interferences  initially.  The  middle  cylinder  is  oriented  to  a  position  not  parallel  with  the 
other  cylinders.  The  objective  function  was  defined  to  minimize  y.  The  design  variables 
were  so  that  each  cylinder  could  rotate  and  translate  freely  in  all  directions.  The  spatial 
constraints  were  defined  so  that  there  would  be  at  least  .01  inches  between  all  cylinders 
and  between  all  cylinders  and  container  wails.  The  GRG  algorithm  was  employed,  and 
after  just  seven  iterations,  the  optimal  positions  were  found  (see  Figure  7). 

A  second  example  is  seen  in  Figure  8.  This  shows  a  hemi- spherical  shell  as  a 
container  wall.  Several  objects,  namely  two  cylinders,  two  pyramids  and  one  long  box, 
are  placed  in  positions  defining  an  initial  feasible  design.  Each  object  is  free  to  rotate 
and  translate  in  all  directions.  The  objective  function  was  to  minimize  y.  The  optimal 
solution  was  reached  after  16  GRG  iterations  with  no  user  intervention. 

The  above  problems  would  have  been  solved  in  less  iterations  if  better  starting 
designs  were  supplied,  however,  the  robustness  of  PACKER  is  apparent  when  optimal 
solutions  are  found  regardless  of  the  starting  design. 


Figure  8:  An  example  of  PACKER  packaging  several  objects  into  a  hemi-spherical 
container.  The  objective  was  to  minimize  y. 
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Capabilities  added  to  increase  speed,  robustness  and  flexibil¬ 
ity. 

There  were  several  tasks  (in  addition  to  those  discussed  in  the  previous  sections) 
that  were  completed  to  increase  the  robustness,  flexibility  and  efficiency  of  the  overall 
CALIPER  and  PACKER  software.  A  discussion  of  the  additional  tasks  is  given  here. 

Explicit  gradients. 

The  algorithms  in  OPTDES  require  derivatives  of  analysis  functions  with  respect 
to  analysis  variables.  Without  user  supplied  derivative  software  explicitly,  OPTDES 
calculates  derivatives  numerically  using  the  finite  forward  difference  approach: 

dfi_  _  /,( Axj  +X])~  ft{xj) 
dxj  Axj 

where  x;  is  the  jth  design  variable  and  /,  is  the  ith  design  function.  With  this  method, 
each  design  variable  is  perturbed  in  turn  by  a  small  amount  in  order  to  evaluate 

fi(Axj  +  Xj). 

This  requires  n  analysis  calls  to  evaluate  one  gradient  call  at  the  beginning  of  each 
design  iteration  where  n  is  equal  to  the  number  of  design  variables  in  the  problem. 

In  the  optimal  packaging  problem,  the  finite  forward  difference  would  require  an 
analysis  call  for  every  rotation  and  translation  design  variable  for  every  object  in  the 
model  to  get  just  one  gradient  call.  It  was  readily  seen  that  if  we  had  explicit  gradient 
formulae  we  could  save  an  enormous  amount  of  computations  (analysis  calls).  Recall 
that  each  analysis  call  involves  computing  all  distances  between  all  objects  and  between 
all  objects  and  all  container  walls  (a  cpu  intensive  call)  along  with  the  mass  properties  of 
the  composite  model.  Explicit  gradients  of  the  analysis  functions  involved  in  the  optimal 
packaging  problem  were  derived.  The  resulting  algorithms  were  programmed  into  the 
ANAGRA  subroutine  in  OPTDES  for  PACKER.  The  following  is  a  discussion  of  these 
formulae.  The  calculation  begins  by  getting  the  separation  distances  and  interference 
distances  between  all  objects  and  between  all  objects  and  container  wadis  (just  one 
analysis  call)  and  stores  the  end  points  of  each  distance  calculated.  Object  one  is  defined 
as  the  object  with  the  design  variable  being  considered  in  the  gradient  calculation  (i.e. 
the  object  that  would  be  perturbed  if  forward  difference  was  applied)  and  object  two  is 
defined  as  the  object  with  which  the  separation  distance  is  calculated. 

Separation/Interference  distances  with  respect  to  translation. 

A  unit  vector  v  is  defined  by  the  direction  of  the  design  variable  begin  considered 
(i.e.  for  x  translation  v  =  (1,0,0)).  The  end  points  of  each  distance  defines  a  vector 
s.  The  direction  of  s  points  away  from  object  one  and  toward  object  two.  The  value  of 
the  gradient  for  separation  distance  between  object  one  and  object  two  with  respect  to 
translation  is  found  by 

v  •  s. 
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Separation/ Interference  distances  with  respect  to  rotation. 


A  vector  r  is  defined  as  the  vector  pointing  from  the  centroid  of  object  one  to 
the  separation  distance  end  point  on  the  surface  of  object  one.  The  vector  v  is  defined 
by  the  direction  of  the  design  variable  considered  (i.e.  for  rotation  about  the  x  axis 
v  =  (1,0,0)).  Using  r,  v  and  s,  the  gradient  for  separation  with  respect  to  rotation  is 
found  by 

(v  x  r)  •  s. 

Center  of  gravity  coordinates  with  respect  to  translation. 

The  gradients  of  the  three  coordinates  of  the  center  of  gravity  with  respect  to 
translations  are  calculated  by  differentiating  the  summation  formula  for  each  c.g.  coor¬ 
dinate  separately. 

For  the  gradient  of  x  with  respect  to  an  x  translation  of  an  object,  we  start  with 

expanded  it  becomes 

x  —  — (m ixx  -F  m2x2  H - h  m*x*  H - h  m.x,). 

mtot 

The  partial  derivative  with  respect  to  x*  (where  k  represents  the  the  design  variable 
being  considered)  we  get 

dx  mk 
dik  mtot 

The  partial  derivative  of  x  with  respect  to  translation  in  the  y  or  z  directions  is 
zero.  This  means  that  translation  in  y  or  z  has  no  effect  on  the  gradient  of  x. 

Similar  expressions  are  obtained  for  y  and  z. 

Center  of  gravity  coordinates  with  respect  to  rotation. 

The  gradients  of  the  three  coordinates  of  the  center  of  gravity  with  respect  to 
any  rotation  is  zero  in  all  cases  since  we  rotate  about  the  centroid  of  each  object.  The 
formulae  for  x,  y  and  z  do  not  depend  on  the  rotational  degrees  of  freedom. 

Moments  of  Inertia  with  respect  to  translation. 

The  gradients  of  the  moments  of  inertia  with  respect  to  translation  are  calcu¬ 
lated  by  differentiating  the  summation  formula  for  each  moment  of  inertia  component 
separately  with  respect  to  the  design  variable  being  considered. 
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For  the  gradient  of  I„  with  respect  to  a  y  translation,  we  start  with 

4r  =  +  l\) 

i 


expanded  it  becomes 

4*  =  +  z\)  +  •  •  •  +  mk(y2k  +  z\)  +  •  •  -  +  m,(y?  +  z2). 

The  partial  derivative  with  respect  to  yk  (where  k  represents  the  design  variable  being 
considered)  we  get 


With  respect  to  zk  we  get 


With  respect  to  xk  we  get 


Similar  expressions  are  obtained  for  Iyy,  /«,  4*,  4*  and  V 

Moments  of  Inertia  with  respect  to  rotation. 

Obtaining  the  explicit  gradients  of  the  moments  of  inertia  with  respect  to  rotation 
(where  the  rotations  are  in  terms  of  the  quaternion)  is  very  complex  and  requires  many 
computations.  It  was  decided  to  use  the  finite  forward  difference  method  for  this  case 
with  a  Ax,  Ay  or  A z  =  .01  perturbation.  The  only  term  in  the  summation  formulae 
for  I  that  changed,  with  respect  to  a  rotation  design  variable,  is  the  term  involving  the 
design  variable  being  considered.  Only  the  term  involving  the  design  variable  is  modified 
according  to  the  quaternion  transformation  matrix  associated  with  the  perturbation. 
The  gradient  is  then  found  by  the  usual  formula: 

dlzx  _  4»(A yk  +  yk)  -  4r(y*) 

dyk  ±yk 

Similar  equations  for  the  other  terms  cam  be  derived.  This  formula  turned  out  to  be  a 
rather  elegant  algorithm  as  programmed  in  the  ANAGRA  subroutine  and  has  shown  to 
work  well. 

NAVGRAPH  group  commands. 

Several  capabilities  were  necessarily  added  to  NAVGRAPH  to  facilitate  the  oper¬ 
ations  of  CALIPER  and  PACKER.  Most  of  the  added  capabilities  were  in  the  GLOBAL 


dlxx 

dyk 


=  2  mkyk. 


dlxx 

dzk 


=  2mkzk. 


dlxi 

dzk 


=  0. 
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GROUP  menu  to  allow  the  user  the  freedom  to  manipulate  the  solid  objects  into  the  de¬ 
sired  models  for  physical  fit  studies  and  optimal  packaging  problems.  The  users  manual 
shows  the  syntax  of  each  command.  These  commands  are  discussed  here. 

GLOBAL  GROUP  TRANSLATE  allows  the  user  to  translate  an  object  to  a  new 
position  or  replicate  an  object  with  the  translation.  The  user  is  queried  on  group  names 
and  the  x,  y,  z  translation. 

GLOBAL  GROUP  ROTATE  allows  the  user  to  rotate  an  object  to  a  new  position 
or  replicate  an  object  with  the  rotation.  The  user  is  queried  on  group  names  and  the 
x,  y,  z  rotation. 

GLOBAL  GROUP  SCALE  allows  the  user  to  scale  an  object  to  a  new  size  or  to 
replicate  am  object  with  the  scale  values  given.  The  user  is  queried  on  group  names  and 
x,  y,  z  scale  values. 

GLOBAL  GROUP  QROTATE  allows  the  user  to  rotate  an  object  to  a  new  posi¬ 
tion  or  replicate  an  object  with  the  quaternion  representation  of  rotation.  The  user  is 
queried  on  group  names,  vector  of  rotation,  and  magnitude  of  rotation  about  the  vector. 

GLOBAL  GROUP  SNAP  allows  the  user  to  translate  an  object  to  a  new  position 
with  a  snap-to  command.  The  user  is  queried  for  two  group  names  and  two  points  (one 
for  each  group).  The  second  group  is  snapped  to  the  first  group  so  that  the  two  given 
points  occupy  the  same  point. 

The  GLOBAL  GROUP  CALIPER  command  calculates  the  separation  and  inter¬ 
ference  distances  between  all  objects  or  optionally  between  just  two  objects.  The  group 
names  are  printed  out  along  with  the  separation  and/or  interference  distances. 

Once  the  separation  and  interferences  have  been  calculated  the  GLOBAL  GROUP 
HIGHLIGHT  will  display  all  interferences  as  vectors  with  flashing  lines.  The  color  of 
the  flashing  lines  is  supplied  by  the  user. 

GLOBAL  GROUP  DISTANCE  allows  the  user  to  get  the  distance  between  any 
two  points  on  any  two  objects.  The  user  is  queried  on  group  names  and  point  numbers. 

The  ability  to  use  the  local  hardware  options  on  the  Tektronix  4107  (and  4109) 
has  been  made  available  with  the  GLOBAL  DISPLAY  SEGMENT  command.  This 
command  toggles  on  and  off  for  segments  to  be  captured  during  a  NAVGRAPH  display. 
Once  the  segment  has  been  captured  the  local  2-D  zoom  and  pan  may  be  used  on  the 
Tektronix  4107  (and  4109). 

The  ability  to  use  the  local  hardware  options  on  the  Tektronix  4129  has  been 
made  available  with  the  GLOBAL  DISPLAY  HARDWARE  command.  This  command 
toggles  on  and  off  for  segments  to  be  captured  during  a  NAVGRAPH  display.  Once  the 
segment  has  been  captured  the  local  3-D  zoom,  pan,  rotate,  etc.  may  be  used  on  the 
Tektronix  4129. 
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The  capability  to  use  Tektronix  4957-4958  digitizers  was  added  to  the  GLOBAL 
DIGITIZER  menu.  This  allows  easy  point  input  and  other  entity  picking  and  delet¬ 
ing.  The  digitizer  hook-up  and  initialization  needs  the  following  Tektronix  SET  UP 
mode  commands  made  to  insure  proper  function  of  the  digitizer  board  (the  4958). 
PASSIGN  PX:4958 
PBAUD  PX:9600 
PBITS  PX:1  7 
PPARITY  PX:ODD 
PFLAG  PX:NONE 
PEOF  PX:NONE 

NAVGRAPH  ployfil  text  font. 

The  text  font  available  in  NAVGRAPH  was  simple  line  text  only,  and  the  user 
had  no  control  over  location  and  orientation.  NAVGRAPH  was  enhanced  to  include  a 
polyfil  font.  The  user  now  has  control  over  size,  slant,  location,  orientation  and  color  of 
the  text  in  NAVGRAPH.  This  new  font  is  much  more  readable  and  allows  the  user  the 
freedom  to  label  the  various  parts  of  a  model  as  needed. 

Methodologies  to  create  AEROHEAT  shapes  with  NAVGRAPH. 

The  NUMERICAL  FLOW  FIELD  PROGRAM  FOR  AERODYNAMIC  HEAT¬ 
ING  ANALYSIS  (AEROHEAT)  as  described  in  AFFDL-TR-85-3054  has  the  ability 
to  solve  aero  heating  predictions  to  several  exact  geometry  aero-shapes  as  well  as  sur¬ 
faces  generated  from  curve  fitting  routines.  It  was  discovered  that  the  geometric  input 
data,  required  by  the  AERO-HEATING  program,  was  so  different  from  the  geometric 
data  produced  by  NAVGRAPH,  that  it  made  the  development  of  a  translator  between 
NAVGRAPH  and  AEROHEAT  highly  impractical.  However,  most  of  the  exact  ge¬ 
ometry  shapes  defined  in  the  AEROHEAT  users  manual  (namely,  the  tangent  ogive, 
ogive  cylinder,  cone  with  nose  cap,  cylinder  with  nose  cap,  and  hyperboloid)  were  mod¬ 
eled  in  NAVGRAPH.  Figure  9  shows  several  of  these  as  modeled  in  NAVGRAPH.  The 
methodologies  used  to  create  these  shapes  are  documented  in  the  tutorial  in  Part  2. 

Fast  surface  meshing  capability. 

A  command  to  mesh  many  surfaces  in  one  command  was  included  in  the  MESH 
CREATE  ELEMENT  menu.  This  allows  the  user  to  easily  create  a  finite  element  mesh 
on  all  surfaces  in  the  data  base  for  use  in  visualization  with  solid  smooth  renderings. 

MOVIE. BYU  software  supplied. 

The  MOVIE.BYU  software  is  supplied  with  a  Tektronix  4129  driver  that  was 
modified  to  load  the  video  lookup  table  with  a  256  color  ramp  in  one  color.  This 
gives  the  user  the  ability  to  display  a  very  smooth  rendering  of  solid  polygonal  models 
that  look  realistic.  The  basic  color  may  be  changed  without  a  need  to  redraw.  The 
user  is  also  prompted  for  background  and  text  colors  that  may  be  different  than  the 
ramp  color.  MOVIE.BYU  offers  transparencies,  multiple  light  sources,  animation  and 
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other  capabilities  not  yet  available  in  NAVGRAPH.  However,  NAVGRAPH  can  produce 
smooth  renderings  of  finite  element  models. 

GRG  algorithm  modifications. 

The  generalized  reduced  gradient  (GRG)  algorithm  used  in  CALIPER  was  mod¬ 
ified  to  be  more  robust  and  efficient  for  the  specific  problem  of  solving  separation  and 
interference  calculations.  The  modifications  are  given  here. 

The  resetting  of  the  Hessian  matrix  after  n  (n  being  the  number  of  design  vari¬ 
ables)  iterations  was  removed.  This  change  allows  the  GRG  algorithm  to  save  function 
and  constraint  information  to  help  search  for  the  solution  near  the  optimum.  This 
change  helped  to  improve  the  robustness  for  problems  with  very  shallow  and  narrow 
objective  functions. 

The  value  of  EPSOPT  was  changed  from  0.001  to  0.0001.  This  allows  for  an 
acceptable  convergence  in  some  separation  calculations  that  previously  were  not  de¬ 
tecting  interferences  properly.  This  increases  computation  time  slightly  but  increases 
robustness  a  great  deal. 

The  parabolic  fit  routine  was  modified  to  find  the  minimum  point  in  a  separa¬ 
tion  calculation  when  the  minimum  point  is  between  the  first  and  second  point  on  the 
parabola. 

A  last  ditch  mode  was  incorporated  for  the  interference  calculation  to  see  if  the 
equality  constraints  axe  violated  and  sets  a  flag  if  they  are.  This  enables  GRG  to 
enter  a  Newton- Raphson  method  to  satisfy  the  constraints  and  start  the  optimization 
process  with  a  valid  starting  design.  This  eliminates  a  few  function  calls  each  time  an 
interference  iteration  is  started  at  an  infeasible  design  point. 

Solid  groupings  into  objects. 

As  discussed  throughout  this  report,  a  very  important  and  beneficial  capability 
was  added  to  NAVGRAPH  to  allow  the  user  to  group  a  number  of  solids  together  to 
model  complex,  convex  and  non-convex  objects.  This  is  simply  the  boolean  add  function 
as  found  in  CSG. 

Packaging  by  layers. 

The  capability  to  package  by  layers  was  studied  during  this  effort.  There  was 
no  method  implemented  in  software  to  automatically  package  by  layer,  but  a  straight 
forward  heuristic  approach  is  offered  here. 

Step  one:  Model  the  desired  container  wall  with  solids  and  group  them  into  one 
object. 

Step  two:  Model  the  desired  objects  by  creating  the  necessary  solids.  Group  the 
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solids  into  the  desired  objects.  Limit  the  number  of  objects  to  a  reasonable  number 
(e.g.  the  number  that  would  fill  the  first  layer). 

Step  three:  Initialize  the  design  with  PREPAC  and  setup  the  desired  optimal 
packaging  problem  with  SETUP  such  that  the  objects  are  packaged  toward  the  desired 
container  wall. 

Step  four:  Run  PACKER  and  optimally  package  the  objects.  This  can  be  done  as 
outlined  in  the  PACKER  section  of  this  report.  Use  the  GRG  algorithm  or  manually  set 
the  unsealed  design  variables  such  that  the  objects  are  packaged  as  tightly  as  possible 
toward  the  container  wall. 

Step  five:  Once  the  objects  are  packaged  tightly,  exit  PACKER  with  a  NAV- 
GRAPH  exit  via  the  PRO  POS  display  option.  This  will  save  the  current  design  as  a 
NAVGRAPH  database  with  updated  object  positions. 

Step  six:  Execute  NAVGRAPH  with  the  updated  data  base  file  from  step  five. 
Ungroup  all  of  the  solids  in  the  objects  (leaving  the  container  alone)  with  the  GLOBAL 
GROUP  DELETE  command.  Then  group  all  solids  (previously  defined  in  the  original 
objects)  into  one  large  group.  This  group  may  then  be  replicated  with  the  GLOBAL 
GROUP  TRANSLATE  command  with  the  appropriate  translation.  The  appropriate 
translation  would  come  from  inspection  or  by  using  the  CALIPER  command  to  detect 
the  interference  distance  between  the  first  new  object  and  the  second.  Repeat  this  step 
until  the  container  is  full  or  until  the  desired  number  of  objects  are  packaged. 

Other  heuristic  methods  and  trouble  shooting. 

During  the  development  and  testing  of  CALIPER,  PACKER  and  the  data  base 
creation,  there  were  some  activities  found  in  each  area  that  worked  better  than  others. 
A  list  of  heuristical  ’rules  of  thumb’  and  trouble  shooting  guides  were  complied  and 
included  in  the  users  manual/tutorial. 
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Translator  Software  Development 


In  order  for  the  software  package  to  be  most  useful,  it  must  be  able  to  generate 
various  data  base  files  which  can  be  read  by  other  existing  programs.  This  capability  is 
called  forward  translation  and  allows  a  user  to  transfer  finite  element  and/or  geometry 
information  to  other  programs  for  further  development  and/or  analysis.  Similarly,  the 
software  package  must  be  able  to  read  information  generated  by  other  programs.  This 
capability  is  called  reverse  translation  and  allows  a  user  to  read  data  that  was  created, 
altered,  or  analyzed  by  another  program. 

Several  forward  and  reverse  translators  were  developed  under  this  work.  A  user 
interface  was  written  to  facilitate  the  execution  of  these  translators  in  a  fashion  that 
resembles  the  NAVGRAPH  user  interface.  All  translators  that  were  developed  are 
described  here. 

TEKNICAD  To  NAVGRAPH  Translator. 

A  reverse  translator  was  written  to  format  TEKNICAD  data  files  into  a  NAV¬ 
GRAPH  data  base  file.  TEKNICAD  is  a  2-D  drafting  program  marketed  by  Tektronix. 
The  TEKNICAD  to  NAVGRAPH  translator  converts  all  points,  straight  lines  (those 
with  one  line  segment  per  line)  and  circular  arcs  into  like  NAVGRAPH  geometry  en¬ 
tities.  The  user  then  may  use  the  points,  lines  and  arcs  in  a  NAVGRAPH  session  for 
further  model  development. 

Translators  Between  NAVGRAPH  and  VUMADM. 

The  Missile  Aerodynamic  Design  Method  (MADM)  system  calculates  aerody¬ 
namic  characteristics  of  complex  configurations  at  supersonic  and  hypersonic  speeds. 
This  is  the  next  generation  supersonic/hypersonic  aerodynamic  analysis  code  with  ma¬ 
jor  upgrades  to  SHABP  (MARK  IV).  The  pre-  and  post-processor  software  (VUMADM) 
can  prepare  geometry  data  for  use  in  analysis  as  well  as  display  analysis  results. 

The  VUMADM  forward  translator  developed  here  writes  geometry  and  finite  el¬ 
ement  information  generated  in  NAVGRAPH  (including  nodes,  elements,  points,  para¬ 
metric  cubic  lines  and  circular  arcs)  to  a  VUMADM  format  file.  This  file  is  read  by  the 
VUMADM  module  and  can  then  be  displayed,  modified  further,  or  otherwise  prepared 
for  analysis  by  the  MADM  program. 

The  VUMADM  reverse  translator  reads  a  VUMADM  format  file  containing  ge¬ 
ometry  and  finite  element  information  and  writes  it  to  a  NAVGRAPH  data  base  file. 
The  user  is  prompted  for  the  VUMADM  file  name  and  the  program  takes  a  few  sec¬ 
onds  to  several  minutes  to  translate  the  file  depending  on  the  model  size.  A  translated 
VUMADM  file,  containing  the  points  and  lines  defining  the  F-15  nose  and  canopy,  is 
shown  in  Figure  10.  The  finite  element  mesh,  also  translated  is  shown.  A  diagnostic 
file  containing  information  about  the  entities  translated,  entities  unsupported  and  any 
errors  that  occurred  is  also  written  and  is  found  on  disk  with  a  .DIA  extension  to  the 
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Figure  10:  A  VUMADM  model  of  the  F-15  nose  geometric  and  finite  element  models 
read  into  the  NAVGRAPH  data  base. 

file  name  input  above. 

IGES  Translators  Between  NAVGRAPH  and  ANVIL-5000. 

ANVIL  is  a  program  that  can  develop  numerical  controlled  milling  machine  in¬ 
structions  from  geometric  surface  definitions.  ANVIL  can  read  IGES  (Initial  Graphics 
Exchange  Specifications)  neutral  files  containing  geometric  points,  lines,  splines  and 
circular  arcs  developed  by  other  CAD/CAM  systems.  Because  ANVIL  (like  other 
CAD/CAM  programs)  uses  the  IGES  neutral  file  for  data  translation,  having  an  IGES 
translator  in  NAVGRAPH  is  important.  Both  forward  and  reverse  IGES  translators 
were  developed. 

The  forward  translator  developed  here  writes  geometric  information  generated 
in  NAVGRAPH  such  as  points,  lines,  splines  and  circular  arcs  to  IGES  Version  3.0 
format.  This  translator,  made  up  of  several  subroutines,  represents  a  major  modification 
to  the  previous  IGES  translator  that  was  available  in  NAVGRAPH.  This  translator 
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Figure  11:  NAVGRAPH  model  via  the  IGES  translator. 

requests  the  name  of  the  IGES  file  to  be  output  and  the  units  of  the  model  (i.e.  inches, 
meters,  etc.).  A  list  of  current  global  section  parameters  then  appears  on  the  screen  for 
acceptance  by  default  or  for  selective  editing.  The  global  section  of  an  IGES  file  contains 
information  describing  the  preprocessor  and  information  needed  by  the  postprocessor 
to  handle  the  file.  Upon  acceptance  of  the  global  section  parameters,  the  IGES  file  is 
written  and  the  program  returns  to  the  forward  translation  menu.  If  an  error  occurs  in 
opening  the  IGES  file,  an  error  message  is  given  and  returns  to  the  forward  translation 
menu. 


Note  that  the  IGES  parametric  spline  entity  (entity  112)  definition  written  by 
this  translator  is  slightly  modified  in  order  that  the  program  ANVIL-5000  can  interpret 
it  correctly.  The  definition,  however,  is  sufficiently  general  and  can  be  read  by  other 
programs  with  IGES  translators. 

The  IGES  reverse  translator  developed  reads  geometric  information  (points,  lines, 
splines  and  circular  arcs)  from  an  IGES  Version  3.0  file  to  a  NAVGRAPH  data  base 
file.  The  user  inputs  the  name  of  the  IGES  file  to  be  read  when  prompted.  The  time 
required  to  translate  the  file  is  dependent  on  the  system  used  and  the  size  of  the  model. 
Large  models  can  take  several  minutes  to  translate.  Figure  11  is  a  NAVGRAPH  plot  of 
several  points,  lines,  parametric  cubic  lines,  and  arcs  translated  from  an  IGES  file. 

The  program  keeps  track  of  the  type  and  number  of  entities  translated,  entities 
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unsupported  in  NAVGRAPH  and  any  errors  that  occurred  in  translation.  This  infor¬ 
mation  is  written  to  a  diagnostic  file  with  the  extension  .IGR  on  the  input  file  name. 
Examination  of  the  diagnostic  file  can  aid  in  reproducing  lost  entities  or  checking  for 
errors. 

Note  that  points  associated  with  straight  lines  axe  written  twice  to  the  NAV¬ 
GRAPH  data  base  since  the  IGES  definition  of  a  line  is  given  by  its  endpoints.  The 
endpoints  of  spline  segments  and  points  that  make  up  arc  segments  axe  only  written 
once. 


Once  a  NAVGRAPH  file  has  been  translated  to  an  ANVIL  file,  one  can  use  the 
data  defining  the  lines  to  create  surfaces  and  instruct  ANVIL  to  create  the  numerical 
machining  instructions  to  manufacture  the  desired  aero-shape. 

Translator  for' AFATL/ ASE 

The  AFATL/ ASE  format  is  the  required  format  for  input  files  to  such  codes  as 
radar  signature  prediction  programs  used  by  the  U.S.  Air  Force.  The  forward  transla¬ 
tor  developed  writes  the  finite  element  information  generated  in  NAVGRAPH  (nodal 
coordinates  and  element  connectivity)  to  the  AFATL/ASE  format.  The  resulting  file 
contains  1)  target  components  (logical  group  of  triangular  facets  or  elements),  2)  nodal 
point  coordinates,  and  3)  element  connectivity  (or  link- lists).  The  subroutine  named 
NAV2ASE  writes  two  output  files.  One  file  is  an  unformatted  (or  binary)  file  containing 
the  above  mentioned  information.  The  name  of  this  file  is  entered  by  the  user  when 
prompted  for  by  the  program.  The  second  file  is  a  formatted  (ASCII)  file  which  contains 
a  list  of  component  numbers  and  the  beginning  and  ending  element  numbers  associated 
with  each  component.  The  name  of  this  file  is  the  same  as  the  entered  file  name  except 
with  a  .TRG  extension  which  identifies  it  as  a  target  component  list  file.  If  an  error 
is  encountered  in  opening  either  file,  an  error  message  is  sent  to  the  screen  and  the 
program  returns  to  the  forward  translation  menu. 

Note  that  for  every  skip  in  element  numbers  found  in  the  NAVGRAPH  data  base 
a  new  component  is  formed  and  written  to  the  AFATL/ASE  format  file.  Only  linear 
triangular  surface  elements  are  accepted  by  AFATL/ASE. 

The  AFATL/ASE  format  reverse  translator  reads  an  AFATL/ASE  format  file 
(containing  nodal  and  triangular  element  information)  into  storage  arrays  from  which 
component  information  is  selectively  extracted  and  rewritten  to  a  NAVGRAPH  data 
base  file.  Component  information  may  be  seen  in  the  default  filename.TRG  file  which 
is  associated  with  the  AFATL/ASE  format  data  base  file.  The  program  requests  the 
name  of  the  file  to  be  read  and  then  prompts  the  user  for  the  beginning  and  ending 
component  numbers  to  be  translated.  A  message  is  sent  to  the  screen  indicating  that 
a  particular  component  is  being  translated.  As  with  the  forward  translator  appropriate 
error  messages  are  sent  to  the  screen  if  errors  are  found  in  opening  files.  Since  a  single 
node  can  be  referenced  several  times  for  a  group  of  adjacent  elements,  a  simple  technique 
was  devised  to  write  each  node  only  once.  Each  time  a  node  is  placed  in  the  NAVGRAPH 
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data  base  its  number  is  placed  in  an  array.  Before  the  next  node  is  written,  it  is  checked 
against  existing  nodes  in  the  array.  If  present,  the  node  is  not  written  to  the  data  base, 
otherwise  it  is.  This  saves  considerable  time  because  the  data  base  routines  sort  the 
existing  data  each  time  they  are  called 

It  is  important  to  note  that  the  coordinates  defining  nodes  are  also  written  as 
points  to  the  NAVGRAPH  data  base.  This  allows  a  user  to  access  these  points  (which 
are  geometric  entities)  in  NAVGRAPH  for  the  creation  of  surface  and  solid  entities. 
The  surfaces  can  then  be  remeshed  with  different  elements  or  the  solid  meshed  with 
solid  hexahedrons  elements.  The  surfaces,  solids  and  new  meshes  axe  then  available  for 
a  wide  range  of  analysis  and/or  further  development. 

The  U.S.  Air  Force  typically  has  models  of  aircraft  and  missiles  in  the  AFATL/  ASE 
data  base.  These  models  can  now  be  translated  to  NAVGRAPH  using  this  reverse 
translator.  The  user  may  define  geometry  entities  in  NAVGRAPH.  Once  the  geometry 
entities  are  defined  in  NAVGRAPH,  they  may  be  used  in  compatibility  studies  and/or 
optimal  packaging  problems.  This  is  the  method  used  in  generating  the  F-15  model  in 
the  next  section.  Note  that  this  does  require  user  interaction,  and  may  require  much 
time  and  expertise  in  using  the  points  supplied  by  the  AFATL/ASE  data  base  to  define 
the  necessary  lines,  surfaces  and  solid  entities.  Several  ’tricks-of-the-trade’  and  other 
guidelines  found  to  work  using  this  procedure  are  given  in  the  tutorial. 


Data  Base  Generation  with  NAVGRAPH. 

A  data  base  consisting  of  an  F-15  aircraft  with  equipment,  several  weapons  and 
submunitions  was  generated  with  NAVGRAPH  for  test  cases  in  1)  optimal  packaging  of 
submunitions  in  missiles,  2)  studying  the  physical  fit  compatibility  of  stores  on  aircraft. 
3)  finite  element  model  generation  and  4)  data  base  format  translation. 

Aircraft  Data  Base  Generation. 

After  some  discussion  of  which  aircraft  to  model,  and  how  to  model  it,  it  was 
jointly  decided  by  Eglin  AFB  and  APTEK  that  the  preferred  method  would  be  to 
use  an  F-15  model  with  data  reformatted  from  the  AFATL/ASE  data  base.  This  was 
accomplished  using  the  ASE  to  NAVGRAPH  translator  and  reading  in  all  nodes  and 
elements  from  the  F-15  ASE  data  file  into  NAVGRAPH.  Each  node  was  also  defined  as 
a  point  for  use  in  defining  lines,  surface  and  solids  in  NAVGRAPH. 

It  was  necessary  to  fit  the  points  to  parametrically  defined  lines  and  surfaces  by 
hand  with  NAVGRAPH  line  and  sur<ace  commands.  Solids  were  then  defined  with  lines 
and  surfaces  where  appropriate. 
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Figure  12:  F-15  ’clean’  and  ’swept’  configurations. 


NAVGRAPH  with  CALIPER  can  now  be  used,  with  the  F-15  model  and  weapons, 
to  perform  physical  fit  studies  of  loadouts. 

Two  models  of  the  F-15  were  created,  one  with  its  equipment  and  gear  in  a 
’clean’  nominal  flight  configuration  and  another  with  its  equipment  and  gear  in  a  ’swept’ 
configuration  (see  Figure  12).  The  ’swept’  configuration  has  the  flaps  and  landing  gear 
modeled  with  NAVGRAPH  by  sweeping  the  flap  and  landing  gear  solids  through  the 
deployed  articulation. 

Because  the  F-15  aircraft  was  modeled,  the  following  equipment  was  also  modeled. 

•  MAU-12  Bomb  Rack 

•  LAU-117  Rail  Launcher 

•  610-gal.  Fuel  Tank 

Figure  13  displays  several  of  these  models. 
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Weapon  Data  Base  Generation. 

NAVGRAPH  was  used  to  generate  a  geometric  model  of  the  external  geometry 
of  each  of  the  following  weapons.  This  data  base  was  created  in  sufficient  detail  to 
determine  store-to-store  or  store- to- aircraft  separation  or  interference. 


•  MK-82  1380901 

•  MK-84  ASIM 

•  BLU-109  8463195 

•  MAVERICK  AGM65  ASIM 

•  SUU-64/65  777108 

•  Sparrow  ASIM 

•  AIM-9  ASIM 


Figure  14  shows  the  geometric  models  of  several  of  the  above  weapons. 

The  SUU-64/65  submunition  dispenser  was  modeled  with  its  internal  geometry 
also,  in  sufficient  detail  for  optimal  packaging  of  submunitions.  Models  containing  fins 
both  deployed  and  stowed  were  created. 

NAVGRAPH  was  used  to  generate  two  models  of  the  external  geometry  of  each  of 
the  following  weapons.  One  file  contains  the  weapon  with  its  wings  extended,  and  the 
other  file  contains  the  weapon  with  its  wings  retracted.  This  data  base  was  also  created 
with  sufficient  detail  to  determine  the  store  to  store  or  store  to  aircraft  separation  and/or 
interference.  Since  the  following  weapons  are  modular  in  nature,  their  modularity  was 
exploited  for  economy. 


•  GBU-10C/B 

•  GBU-10G/H/J/ 

•  GBU-12B/B 

•  GBU-22 

•  GBU-24 

•  GBU-24A/B 
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;  GBU-10C/B  and  the  GBU-24-A/B 


Figure  16:  NAVGRAPH  models  of  the  GATOR  and  CEB  submunitions. 

Figure  15  shows  the  geometric  models  of  several  of  the  weapons  mentioned  above. 


Submunition  Data  Base  Generation. 

NAVGRAPH  was  also  used  to  generate  a  data  base  of  the  following  submunitions. 
The  data  base  was  created  to  provide  sufficient  detail  to  determine  the  submunition  to 
8ubmunition  and  submunition  to  dispenser  separation  and  interference  distances  for 
optimal  packaging.  Two  models,  one  detailed  and  one  simple,  of  the  GATOR,  the  CEB 
and  the  BLU  63  are  provided.  The  detailed  models  contain  all  of  the  curved,  complex 
geometry  (see  Figure  16).  The  simple  models  contain  a  minimum  of  detail  needed  for 
the  optimal  packaging  problem. 


•  BLU-92  GATOR  9299710 

•  BLU-97  CEB  X809420 

•  BLU-63  68C6001 

•  BLU-106  BKEP  162-72168 

•  HTW  PENETRATOR  SK5433186 


50 


Summary. 

The  work  completed  and  documented  here  has  described  the  software  tools  that 
were  developed  and  delivered  to  help  the  Air  Force  weapon  designer  take  a  weapon  from 
concept  to  test  specimen  quickly  and  accurately.  The  capabilities  include: 

1.  Creation  of  geometric  models  of  aircraft,  weapons,  submunitions,  pylons,  racks, 
etc.  to  the  necessary  detail. 

2.  Recall  geometric  models  from  the  an  existing  data  base  to  define  complex,  convex 
and  non-convex  shapes  for  a  desired  model. 

3.  Study  the  physical  fit  compatibility  of  stores  on  aircraft. 

4.  Calculate  'the  mass  properties  (including  surface  area,  volume,  mass,  center  of 
gravity  and  moments  of  inertia)  of  the  desired  model. 

5.  Optimally  package  submunitions  and/or  internal  components  into  a  dispenser  mis¬ 
sile  cargo  bay  with  constraints  on  mass  properties. 

6.  Generate  finite  element  models  of  geometric  models  for  use  in  prediction  software 
codes. 

7.  Format  geometric  data  for  use  by  numerically  controlled  milling  machines  to  ma¬ 
chine  test  specimens. 

8.  Accept  geometric  and  finite  element  data  from  other  Air  Force  data  base  formats. 

9.  And  generate  computer  color  graphic  displays  (line  drawings  or  smooth  rendering) 
of  the  models  including  interference  highlighting. 


This  section  has  discussed  the  technical  aspects  of  the  work  completed  under 
the  contract.  The  next  section  contains  the  users  manual  and  tutorial  for  all  software 
modules  developed  under  the  contract. 
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Part  2 

The  Users  Manual  and  Tutorial 


52 


The  Users  Manual 


This  section  contains  the  specific  command  syntax  for  all  new  capabilities  added 
to  NAVGRAPH  and  the  Eglin  data  base  translators.  These  commands  supplement 
the  users  manuals  for  both  NAVGRAPH  and  OPTDES  (which  accompany  this  users 
manual). 

The  file  named  MSTAR.COM  (a  VMS  DCL  command  file)  is  the  user  friendly 
menu  selector  for  the  many  executable  modules  of  the  software  package.  Note  that 
MSTAR  stand  for  MOVIESTAR  which  is  another  name  (the  commercial  name)  for 
NAVGRAPH.  Upon  issuing  the  ©MSTAR.COM  command,  the  following  will  appear 
on  the  screen: 


Select 

a  NAVGRAPH  option: 

1. 

NAV09 

Run  NAVGRAPH  TEK4109  (4208) 

2. 

NAV29 

Run  NAVGRAPH  TEK4129 

3. 

FOR (MAT) 

Run  NAVGRAPH  Eglin  Data  Formatter 

4. 

PRE(PAC) 

Run  PREPAC  (CONVEN) 

5. 

SETUP 

Run  SETUP 

6. 

PACO  9 

Run  PACKER  TEK4109 (DESIGN) 

7. 

PAC29 

Run  PACKER  TEK4129 (DESIGN) 

8. 

CAL(IPER) 

Run  CALIPER 

9. 

MASS (PROP) 

Run  MASS  PROPERTIES 

10. 

DBM(ERGE) 

Run  DATABASE  MERGE 

11. 

DIR(ECTORY) 

Local  Directory  Listing 

12. 

DCL 

Enter  DCL  Command 

13. 

Q(UIT) 

Quit  NAVGRAPH 

Option 

One  simply  enters  the  desired  option  and  the  chosen  module  is  executed.  Every 
module  is  available  to  execute  at  this  time.  You  can  run  DCL  command  from  inside 
MSTAR.COM  as  well  as  look  at  the  directory. 

Upon  exiting  MSTAR.COM,  the  following  banner  is  displayed.  Care  was  taken 
programming  the  numerous  algorithms  that  make  up  this  software  package.  Also,  there 
were  many  test  cases  run  to  verify  correct  execution  of  all  programs.  However,  there  is 
no  assurance  '  rat  there  are  no  bugs,  but  we  feel  that  they  are  minor.  If  there  are  any 
bugs  found,  please  report  them. 
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**************************************** 


*  * 

*  Report  bugs  to:  * 

*  Hark  Landon  * 

*  APTEK  Inc.  * 

*  1257  Lake  Plaza  Dr.  * 

*  Colorado  Springs,  CO  80906-3578  * 

*  (719)  576-8100  * 

*  * 


**************************************** 


The  new  NAVGRAPH  commands  are  given  in  the  next  21  pages.  They  axe  fol¬ 
lowed  by  the  Eglin  data  base  translator  commands. 
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GLOBAL-COMMANDS  DIGITIZER 


MENU: 

GLOB 

1.  CLEAR 

2.  COLOR  => 

3.  DEFAULTS  => 


4.  DIGITIZER  => 


5. 

DISPLAY-OPTIONS 

=  > 

6. 

DRAW 

7. 

ERASE 

=> 

8. 

EXIT 

9. 

GROUP 

=  > 

10. 

HARDCOPY 

=  > 

11. 

HELP 

12. 

LABEL 

=  > 

13. 

LIST 

=  > 

14. 

MENU 

15. 

PLOT 

=> 

16. 

RECALL 

17. 

RESET 

18. 

TIME 

19. 

VIEW-OPTIONS 

=> 

ENTER:  ->  4  <CR>  or 

DIG  <CR> 

KEYBOARD:  GLOB  DIGI 


DESCRIPTION: 

The  DIGITIZER  command  selects  the  DIGITIZER  menu. 


2-xx  GLOBAL  COMMANDS 
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GLOBAL-COMMANDS  .DIGITIZER.DETECT 


MENU:  GLOB.D1GI 

ENTER:  ->  1  <CR>  or  DET  <CR> 


KEYBOARD:  GLOB  DIGI  DETE 


DESCRIPTION: 

The  DETECT  command  toggles  ON  and  OFF  to  enable  or  disable  graphic  input  from  a  device 
such  as  a  digitizing  tablet,  mouse,  or  thumbwheels.  The  default  value  is  OFF. 


PROMPTS: 

Enter  on  or  off  (default  =  off)> 

EXAMPLES: 

Enter  on  or  off  (default  =  off)>  ON 
Enables  graphic  input. 

Enter  on  or  off  (default  =  off)>  OFF 
Disables  graphic  input. 


2-xx 
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GLOBAL-COMMANDS. DIGITIZER.ORIENT 


MENU:  GLOB.DIGI 

ENTER:  ->  2  <CR>  or  OR  <CR> 


KEYBOARD:  GLOB  DIGI  ORIE 


DESCRIPTION: 

The  ORIENTATION  command  creates  a  transformation  between  screen  coordinates  and  world 
coordinates.  This  is  needed  to  properly  orient  a  drawing  on  a  digitising  tablet  or  some  other  graphic 
input  device.  The  user  will  be  prompted  to  detect  any  three  points  using  the  graphic  input  device 
and  to  enter  the  corresponding  world  coordinates  (x  and  y)  for  each  point. 


PROMPTS: 

Detect  point  1  on  digitizer 
Enter  pointl  (x,y)> 

Detect  point2  on  digitizer 
Enter  pomt2  (x,y)> 

Detect  point3  on  digitizer 
Enter  point3  (x,y)> 


EXAMPLES: 

Detect  pointl  on  digitizer 
Enter  pointl  (x,y)>  0.0  0.0 

Detect  point2  on  digitizer 
Enter  point2  (x,y)>  10.0  0.0 

Detect  point3  on  digitizer 
Enter  point3  (x,y)>  5.0  20.0 

Creates  2-d  transformation  between  screen  and  drawing  using  the  three  point  0.0, 0  0  10.0,0.0 
5.0,20.0. 


2- xx 
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GLOBAL-COMMANDS.DIGITIZER.Z- VALUE 

MENU:  GLOB.DIGI 

ENTER:  ->  3  <CR>  or  Z-  <CR> 


KEYBOARD:  GLOB  DIGI  Z-VA 


DESCRIPTION: 

The  Z- VALUE  command  is  used  to  enter  a  constant  z- value  for  points  enter  with  a  graphic  input 
device. 


PROMPTS: 


Enter  constant  z  value 
(default  =  0.00000E+00)> 


EXAMPLES: 

Enter  constant  z  value 
(default  =  0.00000E+00)>  5.5 

A  z  coordinate  of  5.5  will  be  given  to  every  point  entered  with  a  graphic  input  device  until  a  new 
z  value  is  given. 


2- xx 
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GLOBAL-COMMANDS.DIGITIZER.DEVICE 


MENU: 

GLOB.DIGI 

1.  DETECT 

2.  ORIENTATION 

3.  Z-VALUE 


4.  DEVICE  => 


ENTER:  ->  4  <CR>  or  DEV  <CR> 


KEYBOARD:  GLOB  DIGI  DEVI 


DESCRIPTION: 

The  DEVICE  command  selects  the  DEVICE  menu. 


2- xx 
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GLOBAL-COMMANDS.DIGITIZER.DEVICE.MOUSE 


MENU:  GLOB.DIGI.DEVI 
ENTER:  ->  1  <CR>  or  MO  <CR> 


KEYBOARD:  GLOB  DIGI  DEVI  MOUS 


DESCRIPTION: 

The  MOUSE  command  identifies  a  mouse  as  the  graphic  input  device. 


2- xx 
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GLOBAL-COMMANDS.DIGITIZER.DEVICE.TABLET 

MENU:  GLOB.DIGI.DEVI 
ENTER:  ->  2  <CR>  or  TA  <CR> 


KEYBOARD:  GLOB  DIGI  DEVI  TABL 


DESCRIPTION: 

The  TABLET  command  identifies  a  digitizer  tablet  as  the  graphic  input  device. 


2- xx 
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GLOB  AL-COMMANDS.DIGITIZER.DEVICE.THUMB  WHEELS 

MENU:  GLOB.DIGI.DEVI 
ENTER:  ->  3  <CR>  or  TH  <CR> 


KEYBOARD:  GLOB  DIGI  DEVI  THUM 


DESCRIPTION: 

The  THUMBWHEELS  command  identifies  the  thumbwheels  or  joydisk  as  the  graphic  input  de¬ 
vice. 


2- xx 
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GLOBAL-COMMANDS  .DISPLAY-OPTIONS. TEXT 


MENU:  GLOB.DISP 

ENTER:  ->  26  <CR>  or  TEX  <CR> 


KEYBOARD:  GLOB  DISP  TEXT 


DESCRIPTION: 

The  TEXT  command  prompts  the  user  for  text  information.  Text  can  be  used  for  comments 
dimensions,  etc..  There  are  two  text  fonts,  polygonal  and  line.  The  text  is  two  dimensional  but  can 
be  manipulated  in  three  dimensions.  When  maniplating  text  in  2-d,  screen  coordinates  (0  to  1)  are 
used  for  positioning.  When  manipulating  text  in  3-d,  world  coordinates  are  used.  Only  one  type  of 
text  (polygonal  or  line)  can  be  displayed. 


PROMPTS: 

Polygonal  or  Line  Text?  (  P  or  (L)  ] 

> 

Text  Message  of  Text  String  X 

()•' 

> 

Width  (1.0)  and  Height  (1.0)  :> 

2-D  or  3-D  Text  Position  (2):> 
if  2-D 

Position  X  (  0.00),  Y  (  0.00),  Angle  (  0.00)  :> 
if  3-D 

Position  Origin  X  (  0.00),  Y  (  0.00),  Z  (  0.00)  :> 

Position  Vector  X  (  0.00),  Y  (  0.00),  Z  (  0.00)  :> 

Angle  to  Slant  Text: 

(default  =  O.OOOOOE+O0)> 

Additional  Prompts: 

Color  Red  (10),  Blue  (1.0),  Green  (1.0)  :> 

Modify  Text  String  X  ?  (default  =  no)> 

Add  Another  Text  String  ?  (default  =  no)> 

Delete  Existing  Text  Strings  ?  (default  =  no)> 

Modify  Existing  Text  Strings  ?  (default  =  no)> 

Show  Summary  of  Existing  Text  Strings  ?  (default  =  yes)> 
Text  String  Number  to  be  Modified  :> 


2- xx 
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GLOBAL-COMMANDS.DISPLAY-OPTIONS.ENTITY.TEXT 


MENU:  GLOB.DISP.ENTI 

ENTER:  ->  17  <CR>  or  TEX  <CR> 


KEYBOARD:  GLOB  DILP  ENTI  TEXT 


DESCRIPTION: 

The  TEXT  command  toggles  ON  and  OFF  and  turns  the  display  of  text  on  or  off. 


PROMPTS: 

Enter  on  or  off  (default  =  off)> 


EXAMPLES: 

Enter  on  or  off  (default  =  off)>  ON 

Turns  display  of  text  on. 

Enter  on  or  off  (default  =  off)>  OFF 

Turns  display  of  text  off. 


2- xx 
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GLOBAL-COMMANDS. LIST.TEXT 


MENU:  GLOB. LIST 

ENTER:  ->  16  <CR>  or  TEX  <CR> 


KEYBOARD:  GLOB  LIST  TEXT 


DESCRIPTION. 

The  TEXT  lists  ail  existing  text  strings  on  the  screen. 


2- xx 
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GLOBAL-COMMANDS.DISPLAY-OPTIONS. HARDWARE-GRAPHICS 


MENU:  GLOB.DISP 

ENTER:  ->  25  <CR>  or  HA  <CR> 


KEYBOARD:  GLOB  DISP  HARD 


DESCRIPTION: 

Tlie  HARDWARE-GRAPHICS  turns  the  local  hardware  graphics  capabilities  on  and  off. 


PROMPTS: 

Enter  on  or  off  (default  =  off)> 

EXAMPLES: 

Enter  on  or  off  (default  =  off)>  ON 
Turns  hardware  graphics  on. 

Enter  on  or  off  (default  =  off)>  OFF 
Turns  hardware  graphics  off. 


2-xx 
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GLOBAL-COMMANDS.DISPLAY-OPTIONS. SEGMENT 


MENU:  GLOB.DISP 

ENTER:  ->  27  <CR>  or  SE  <CR> 


KEYBOARD:  GLOB  DISP  SEGM 

DESCRIPTION: 

The  SEGMENT  turns  the  local  segment  graphics  capabilities  on  and  off. 

PROMPTS: 

Enter  on  or  off  (default  =  off)> 

EXAMPLES: 

Enter  on  or  off  (default  =  off)>  ON 
Turns  segments  on. 

Enter  on  or  off  (default  =  off)>  OFF 
Turns  segments  off. 


2- xx 
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GLOBAL-COMMANDS.GROUP.CALIPER 


MENU:  GLOB.GROU 

ENTER:  ->  9  <CR>  or  CA  <CR> 


KEYBOARD:  GLOB  GROU  CALI 


DESCRIPTION: 

The  CALIPER  command  executes  the  CALIPER  distance  calculator. 


PROMPTS: 

Calculate  all  interferences  ?  (default  =  no)> 

Input  the  first  group  name  (<CR>  to  quit):> 

Input  the  second  group  name:> 

Do  you  want  to  specify  interference  direction  ?  (default  =  no)> 
Input  x,y,z  components  of  interference  direction  vector:  > 


EXAMPLES: 

Calculate  all  interferences  ?  (default  =  no)> 

AVAILABLE  GROUP  NAMES  ARE: 

G1 

G2 

Input  the  first,  group  name  (<CR>  to  quit):>  G1 
Input  the  second  group  name:>  G2 

Do  you  want  to  specify  interference  direction  ?  (default  =  no)>  Y 
Input  x,y,z  components  of  interference  direction  vector:>  10  0 
Distance  information  will  be  listed  here. 


2- xx 
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GLOBAL-COMMANDS. GROUP-DISTANCE 


MENU:  GLOB.GROU 

ENTER:  ->  10  <CR>  or  DIS  <CR> 


KEYBOARD:  GLOB  GROU  DIST 


DESCRIPTION: 

The  DISTANCE  command  calculates  the  distance  between  two  points. 


PROMPTS: 

Enter  the  first  point  > 
Enter  the  second  point  > 


EXAMPLES: 

Enter  the  first  point  >  5 
Enter  the  second  point  >  10 

Distance  between  two  points  is  :  20.0 


2- xx 
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GLOBAL-COMMANDS.GROUP.HIGHLIGHT-INT 


MENU:  GLOB.GROU 

ENTER:  ->  11  <CR>  or  HI  <CR> 


KEYBOARD:  GLOB  GROU  HIGH 


DESCRIPTION: 

The  HIGHLIGHT-INT  command  turn*  the  interference  highlights  on  and  off.  This  feature  is 
availible  on  Tektronix  machines.  A  highlight  consists  of  a  blinking  line  along  the  vector  of  interfer¬ 
ence.  The  CALIPER  command  must  be  executed  to  calculate  interferences  whenever  changes  are 
made  to  existing  groups. 


PROMPTS: 

Highlight  line  color  components  (  red,  green,  blue)  ?  > 
Do  you  want  to  list  interferences?  (default  =  no)> 


EXAMPLES: 

HiglUiglit  line  color  components  (  red,  green,  blue)  ?  >  1  1  0 
Interferences  will  be  drawn  on  the  screen. 

Do  you  want  to  list  interferences?  (default  =  no)>  y 
Interference  information  will  be  listed. 


2-xx 
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GLOBAL-COMMANDS.GROUP.ROTATE 


MENU:  GLOB.GROU 

ENTER:  ->  12  <CR>  or  RO  <CR> 


KEYBOARD:  GLOB  GROU  ROTA 


DESCRIPTION: 

The  ROTATE  command  rotates  a  group  using  z,  y,  and  z  rotations.  A  group  can  be  rotated 
with  or  without  creating  a  new  group. 


PROMPTS: 

Enter  new  group  name 
(default  =  )> 

Enter  enter  group  to  rotate 
(default  =  )> 

Enter  rotation  angle  (x,y,z)  > 


EXAMPLES: 

To  create  a  new  group: 

Enter  new  group  name 
(default  =  groupl  )>  group2 
Enter  enter  group  to  rotate 
(default  —  group2  )>  groupl 
Enter  rotation  angle  (x,y,z)  >  10  15  20 
To  rotate  a  group: 

Enter  new  group  name 
(default  =  groupl  )>  groupl 
Enter  enter  group  to  rotate  > 

(default  =  groupl  )>  groupl 
Enter  rotation  angle  (x,y,z)  >  20  5  20 


2- xx 
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MENU:  GLOB.GROU 

ENTER:  ->  13  <CR>  or  QR  <CR> 


GLOBAL-COMMANDS.GROUP.QROTATE 


KEYBOARD:  GLOB  GROU  QROT 


DESCRIPTION: 

The  QROTATE  command  rotates  a  group  using  a  quaternion.  A  quaternion  consists  of  a  vector 
and  a  rotation  angle  about  that  vector.  A  group  can  be  rotated  with  or  without  creating  a  new  group. 


PROMPTS: 

Enter  new  group  name 
(default  =  )> 

Enter  enter  group  to  rotate 
(default  =  )> 

Enter  the  center  point  of  the  group  (x,y ,z)  > 
Enter  the  rotation  vector  (x,y,z)  > 

Enter  rotation  angle  (alpha)  > 


EXAMPLES: 

To  create  a  new  group: 

Enter  new  group  name 

(default  =  groupl  )>  group2 

Enter  enter  group  to  rotate 

(default  •=  group2  )>  groupl 

Enter  the  center  point  of  the  group  (x,y,z)  >  5,10,0 

Enter  the  rotation  vector  (x,y,z)  >  1,0,0 

Enter  rotation  angle  (alpha)  >  20 

To  rotate  a  group: 

Enter  new  group  name 

(default  =  groupl  )>  groupl 

Enter  enter  group  to  rotate 

(default  =  groupl  )>  groupl 

Enter  the  center  point  of  the  group  (x,y,z)  >  3,0,1 

Enter  the  rotation  vector  (x,y,z)  >  1,1,1 

Enter  rotation  angle  (alpha)  >  30 
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GLOBAL-COMMANDS. GROUP.SNAP 


MENU:  GLOB.GROU 

ENTER:  ->  14  <CR>  or  SN  <CR> 


KEYBOARD:  GLOB  GROU  SNAP 


DESCRIPTION: 

Tlie  SNAP  command  is  another  method  for  translating  a  group.  A  point  is  selected  in  each  of 
two  groups.  The  second  group  is  then  translated  so  that  the  world  coordinates  of  the  point  in  the 
second  group  are  the  same  as  the  point  in  the  second  first. 


PROMPTS: 

Enter  first  group  name  > 

Enter  a  point  on  first  group  > 
Enter  second  group  name  > 
Enter  a  point  on  second  group  > 


EXAMPLES: 

Enter  first  group  name  >  group  1 
Enter  a  point  on  first  group  >  10 
Enter  second  group  name  >  group2 
Enter  a  point  on  second  group  >  50 


2- xx 


GLOBAL  COMMANDS 
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GLOBAL-COMMANDS.GROUP.TRANSLATE 


MENU:  GLOB.GROU 

ENTER:  ->  15  <CR>  or  TR  <CR> 


KEYBOARD:  GLOB  GROU  TRANS 


DESCRIPTION: 

The  TRANSLATE  command  translates  a  group  using  z,  y,  and  t  translations.  A  group  can  be 
translated  with  or  without  creating  a  new  group. 


PROMPTS: 

Enter  new  group  name 
(default  =  )> 

Enter  enter  group  to  translate 
(default  =  )> 

Enter  translation  (x,y,z)  > 


EXAMPLES: 

To  create  a  new  group: 

Enter  new  group  name 
(default  =  )>  group 2 
Enter  enter  group  to  translate 
(default  =  )>  group  1 
Enter  translation  (x,y,z)  >  5,0,10.5 
To  translate  a  group: 

Enter  new  group  name 

(default  =  )>  group  1 

Enter  enter  group  to  translate 

(default  =  )>  group  1 

Enter  translation  (x,y,z)  >  3.0, 0,0 


2-xx 


GLOBAL  COMMANDS 
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GLOBAL-COMMANDS.GROUP.SCALE 


MENU:  GLOB.GROU 

ENTER:  ->  16  <CR>  or  SC  <CR> 


KEYBOARD:  GLOB  GROU  SCAL 


DESCRIPTION: 

The  SCALE  command  scales  a  group  using  x,  y,  and  z  scale  factors.  A  group  can  be  scaled  with 
or  without  creating  a  new  group. 


PROMPTS: 

Enter  new  group  name 
(default  as  )> 

Enter  enter  group  to  scale 
(default  =  )> 

Enter  scale  factors  (x,y,z)  > 


EXAMPLES: 

To  create  a  new  group: 

Enter  new  group  name 
(default  =  )>  group2 
Enter  enter  group  to  scale 
(default  =  )>  groupl 
Enter  scale  factors  (x,y,z)  >  .9,1,1 
To  scale  a  group: 

Enter  new  group  name 

(default  =  )>  groupl 

Enter  enter  group  to  scale 

(default  =  )>  groupl 

Enter  scale  factors  (x,y,z)  >  .75, .5, .3 


2- xx 


GLOBAL  COMMANDS 
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Eglin  AFB  Translators 

The  Eglin  AFB  tranlators  module  of  NAVGRAPH  translates  data  base  files  to 
and  from  AFATL/ASE,  IGES,  TEKNICAD,  and  VUMADM  formats.  (NOTE:  The 
AFATL/ASE  format  is  used  for  radar  signature  prediction;  the  VUMADM  format  is  re¬ 
quired  for  use  in  the  Missile  Aerodynamic  Design  Method  (formerly  S/HABP)  system.) 

The  VUMADM  translator  provides  for  both  geometry  and  finite  element  entity 
translation,  the  AFATL/ASE  translator  for  finite  element  entity  translation,  and  the 
IGES  and  TEKNICAD  translators  for  geometry  entity  translation  only.  The  IGES 
forward  translator  features  selective  editing  of  the  global  section  parameters  the  control 
file  processing. 
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EGLIN  AFB  TRANSLATORS . FORWARD 


MENU: 

EGLIN  AFB  TRANSLATORS 

1 .  FORWARD  => 

2 .  REVERSE  => 

ENTER:  ->  1  <CR  or  FO  <CR> 


KEYBOARD :  FORW 


DESCRIPTION: 

The  FORWARD  command  selects  the  FORWARD  translation  menu. 


eglin  afb  translators .forward. afatl/ as e 


MENU: 

FORWARD 

1 .  AFATL /ASE 

2 .  IGES 

3 .  TEKNICAD 

4 .  VUMADM 

ENTER:  ->  1  <CR>  or  AF  <CR> 


KEYBOARD:  FORW  AFAT 


DESCRIPTION: 

The  AFATL/ASE  command  creates  an  AFATL/ASE  format  file  from 
the  finite  element  information  (nodal  coordinates  and  element 
connectivity)  of  the  NAVGRAPH  database.  Target  components,  node 
coordinates,  and  element  link-lists  are  written  to  an  AFATL/ASE 
compatible  file. 


PROMPTS : 

Enter  output  file  name 
(default  *  problem  name) > 


EXAMPLES: 

Enter  output  file  name 
(default  =  problem  name)>  PYLON. ASE 

Creates  the  AFATL/ASE  format  file  PYLON. ASE  from  the  finite 
element  information  in  the  NAVGRAPH  database. 


EGLIN  AFB  TRANSLATORS . FORWARD . IGES 


MENU: 

FORWARD 

1 .  AFATL/ASE 

2.  IGES 

3 .  TEKNICAD 

4 .  VUMADM 

ENTER:  ->  2  <CR>  or  IG  <CR> 


KEYBOARD:  FORW  IGES 


DESCRIPTION: 

The  IGES  command  creates  an  IGES  file  from  the 
geometry  of  the  NAVGRAPH  database.  Points,  lines,  splines,  and 
circular  arcs  are  written  to  an  IGES  compatible  file. 


PROMPTS : 

Enter  units 
(default  -  INCH) > 

Enter  output  file  name 
(default  =  problem  name)> 


EXAMPLES : 

Enter  units 
(default  ~  INCH) >  MM 

Enter  output  file  name 

(default  =  problem  name)>  NEWFILE.IGE 

Creates  the  IGES  file  NEWFILE.IGE  from  the  geometric 
information  in  the  NAVGRAPH  database  and  gives  it  units  of 
millimeters . 


EGLIN  AFB  TRANSLATORS . FORWARD . TEKNICAD 


MENU: 

FORWARD 

1 .  AFATL/ASE 

2 .  IGES 

3.  TEKNICAD 

4 .  VUMADM 

ENTER:  ->  3  <CR>  or  TE  <CR> 


KEYBOARD :  FORW  TEKN 


DESCRIPTION: 

The  TEKNICAD  command  creates  a  TEKNICAD  file  from  the 
geometry  of  the  NAVGRAPH  database.  2-D  points,  lines,  splines, 
and  circular  arcs  are  written  to  a  TEKNICAD  compatible  file. 


PROMPTS: 

Enter  output  file  name 
(default  =«  problem  name)> 


EXAMPLES : 

Enter  output  file  name 

(default  *  problem  name)>  PLATE. TEK 

Creates  the  TEKNICAD  file  PLATE. TEK  from  the  geometric 
information  (points,  lines,  splines  and  circular  arcs)  in  the 
NAVGRAPH  database. 


EGLIN  AFB  TRANSLATORS . FORWARD . VUMADM 


MENU: 

FORWARD 

1 .  AFATL/ASE 

2 .  IGES 

3 .  TEKNICAD 

4.  VUMADM 

ENTER:  ->  4  <CR>  or  VU  <CR> 


KEYBOARD:  FORW  VUMA 


DESCRIPTION: 

The  VUMADM  command  creates  a  VUMADM  format  file  from  the 
geometry  and  finite  element  information  of  the  NAVGRAPH  database. 
Points,  lines,  splines,  circular  arcs,  nodes,  and  elements  are 
written  to  a  VUMADM  compatible  file. 


PROMPTS: 

Enter  output  file  name 
(default  -  problem  name)> 


EXAMPLES: 

Enter  output  file  name 
(default  =  problem  name)>  MISSLE.MAD 

Creates  the  VUMADM  file  MISSLE.MAD  from  the  geometric  and 
finite  element  information  in  the  NAVGRAPH  database. 
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EGLIN  AFB  TRANSLATORS . REVERSE 

MENU: 

EGLIN  AFB  TRANSLATORS 

1 .  FORWARD  •> 

2 .  REVERSE  »> 

ENTER:  ->  2  <CR  or  RE  <CR> 

KEYBOARD :  REVE 

DESCRIPTION: 

The  REVERSE  command  selects  the  REVERSE  translation  menu. 


EGLIN  AFB  TRANSLATORS . REVERSE . AFATL/ASE 


MENU : 

REVERSE 

1 .  AFATL/ASE 

2 .  IGES 

3 .  TEKNICAD 

4 .  VUMADM 

ENTER:  ->  1  <CR>  or  AF  <CR> 


KEYBOARD :  REVE  AFAT 


DESCRIPTION: 

The  AFATL/ASE  command  translates  an  AFATL/ASE  format  file 
containing  finite  element  information  to  a  NAVGRAPH  database. 
Nodal  coordinates  and  element  link-lists  are  translated  from 
AFATL/ASE  format  to  NAVGRAPH  format . 


PROMPTS : 

Enter  input  file  name 
(default  *  problem  name)> 


EXAMPLES: 

Enter  input  file  name 

(default  *  problem  name)>  PYLON. ASE 

Copies  the  finite  element  and  nodal  coordinate  information 
from  the  AFATL/ASE  file,  PYLON. ASE,  to  the  current  NAVGRAPH 
database . 


EGLIN  AFB  TRANSLATORS . REVERSE . IGES 


MENU: 

REVERSE 

1 .  AFATL/ASE 

2.  IGES 

3 .  TEKNICAD 

4 .  VUMADM 

ENTER:  ->  2  <CR>  or  IG  <CR> 


KEY30ARD :  REVE  IGES 


DESCRIPTION: 

The  IGES  command  translates  an  IGES  file  to  a  NAVGRAPH 
database.  Points,  lines,  splines,  and  circular  arcs  are  trans¬ 
lated  from  IGES  to  NAVGRAPH  format. 


PROMPTS: 

Enter  input  file  name 
(default  -  problem  name)> 


EXAMPLES: 

Enter  input  file  name 

(default  -  problem  name)>  NEWFILE . IGE 

Copies  the  geometric  information  (points,  lines,  splines  and 
circular  arcs)  from  the  IGES  file,  NEWFILE. IGE,  to  the  current 
NAVGRAPH  database. 
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EGLIN  AFB  TRANSLATORS . REVERSE . TEKNICAD 


MENU: 

REVERSE 

1 .  AFATL/ASE 

2 .  IGES 

3.  TEKNICAD 

4 .  VUMADM 

ENTER:  ->  3  <CR>  or  TE  <CR> 


KEYBOARD :  REVE  TEKN 


DESCRIPTION: 

The  TEKNICAD  command  translates  a  TEKNICAD  file  to  a  NAVGRAPH 
database.  2-D  Points,  lines,  splines,  and  circular  arcs  are  trans¬ 
lated  from  TEKNICAD  format  to  NAVGRAPH  format. 


PROMPTS: 

Enter  input  file  name 
(default  -  problem  name)> 


EXAMPLES: 

Enter  input  file  name 

(default  -  problem  name)>  PLATE. TEK 

Copies  the  geometric  information  (2-D  points,  lines,  splines 
and  circular  arcs)  from  the  TEKNICAD  file,  PLATE. TEK,  to  the 
current  NAVGRAPH  datbase  . 
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EGLIN  AFB  TRANSLATORS . REVERSE . VUMADM 


MENU: 


REVERSE 

1 .  AFATL/ASE 

2 .  IGES 

3 .  TEKNICAD 

4.  VUMADM 


ENTER:  ->  4  <CR>  or  VU  <CR> 


KEYBOARD :  REVE  VUMA 


DESCRIPTION: 

The  VUMADM  command  translates  a  VUMADM  format  file  to  a 
NAVGRAPH  database.  Points,  circular  arcs,  parametric  cubic  lines, 
nodes,  and  elements  are  translated  from  VUMADM  format  to  NAVGRAPH 
format . 


PROMPTS : 

Enter  input  file  name 
(default  «  problem  name)> 


EXAMPLES: 

Enter  input  file  name 

(default  =  problem  name)>  MISSLE.MAD 

Copies  the  geometric  (points,  circular  arcs,  and  parametric 
cubic  lines) ,  nodal  coordinate,  and  finite  element  information 
from  the  VUMADM  file,  MISSLE.MAD,  to  the  current  NAVGRAPH 
database . 
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The  Tutorial 


The  tutorial  presented  here  is  a  collection  of  actual  runs  with  the  various  modules 
of  the  software  package. 

First,  we  give  several  NAVGRAPH  modeling  examples.  Because  every  module 
needs  a  geometric  model  to  work  with,  it  is  important  that  the  user  becomes  very  good 
at  using  NAVGRAPH.  The  more  the  user  learns  about  the  commands  and  capabilities  of 
NAVGRAPH,  the  better  and  faster  he  will  be  able  to  solve  the  problem  at  hand.  Among 
the  examples  given  are  the  basic  exact  geometry  shapes  as  found  in  the  AERO  HEAT 
users  manual.  Also,  the  steps  to  create  a  MAVERICK  missile  are  provided.  These 
examples  are  in  the  navgraph.RUN  format.  They  can  be  typed  line-by-line  interactively. 
They  can  also  be  run  automatically  with  the  RECALL  command  in  NAVGRAPH. 

There  are  other  modeling  examples  throughout  that  the  user  may  glean  some 
pointers  on  how  to  model  better  and  faster. 

The  third  section  contains  an  example  with  CALIPER  in  the  physicad  fit  com¬ 
patibility  problem.  In  this  example,  the  user  is  given  models  of  an  F-15  wing,  pylon, 
runway  and  an  MK84  weapon.  It  is  desired  to  place  the  weapon  on  the  wing/pylon 
aissembly  so  that  there  are  no  interferences. 

The  fourth  section  contains  several  optimal  packaging  examples.  These  examples 
take  the  user  through  model  creation,  meshing,  mass  properties  calculation,  design 
setup  and  optimal  packaging  with  PACKER.  There  are  many  heuristic  pointers  and 
discussions  in  the  right  margin  of  the  example  outputs. 

Finally,  there  is  a  generad,  useful,  modeling  guideline  providing  general  rules  of 
thumb  aind  pointers  for  malting  good  geometric  models.  A  trouble  shooting  guide  is  also 
provided  that  suggests  some  proven  ’work  arounds’. 
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This  file  creates  a  cylinder  with  a  nose  cap.  The  model  is  shown  here. 


DEFI 
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10 

270  0  0 

3 

SWEE 

14 

4 
3 
90 

ROTA 

15 
14 

270  0  0 
3 

ENTE 
18 
1  1 


SOLI 

LINE 

1 

17  16  15  14 
13  12  11  10 
2  9  8  7 
LINE 

2 

13  12  11  10 

18  18  18  18 
16  5  4 
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This  file  creates  a  tangent  ogive  body  of  revolution.  The  model  is  shown. 


DEFI 

POIN 

ENTE 

1 

0  0  0 
2 

.177124  .82288  0 

3 

.67712  1.5  0 

4 

2.524431  2.612486  0 

5 
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6 
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0 

LINE 
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1 
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2 
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3 
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1 
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3 
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7 
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3 
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10 

3 
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The  following  are  the  MOVIESTAR  commands  for  creating  a  solid  model  of  a 
MAVERICK  missile. 

DISP  AUTO-CEN  ON 
DEFI 
POIN 
ENTE 
1 

0  0  0 
2 

2.5  0  4 

3 

2.5  2.8284  2.8284 

4 

2.5  4  0 

5 

.6764  0  2.3585 
0 

LINE 
ENTE 
1 

1  1 
1  1 
1  1 
1  1 

ARC 

5 

2  3  4 

ROTA  Lines  for  first  solid. 

6 
5 

270  0  0 

3 

ARC 

9 

1  5  2 
ROTA 

10 
9 

270  0  0 
3 


Turn  on  auto-center. 


Points  for  first  solid. 
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SOLI 
LINE 
1 

12  3  4 
5  6  7  8 
9  10  11  12 

POIN 
ENTE 

24 

14  0  6 

25 

14  4.2426  4.2426  Points  for  second  solid. 

26 

14  6  0 
27 

8.1637  0  5.4963 
0 

LINE 

ARC 

13 

24  25  26 
ROTA 

14 
13 

270  0  0 

3  Lines  for  second  solid. 

ARC 

17 

2  27  24 
ROTA 

18 
17 

270  0  0 

3 

SOLI 
LINE 
2 

5  6  7  8 
13  14  15  16 


Second  solid  made  with 
twelve  lines. 


First  solid  created  from 
twelve  lines. 
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17  18  19  20 


SURF 

TRAN 

13 

12 

83  0  0 

SOLI 

PARA 

3 

12  13 

POIN 

ENTE 

54 

31.5  -.5  6 

55 

31.5  .5  6 

56 

84.16  -.5  6 

57 

84.16  .5  6 

58 

84.16  .1  14.26 

59 

84.16  -.1  14.25 
0 

SOLI 

ENTE 

4 

54  55  55  54 
56  57  58  59 

ROTA 

5 

4 

270  0  0 
3 

SURF 

TRAN 


Translate  surface  to  be 
used  in  third  solid. 


Third  solid  made  by  parametric 
fit  between  two  surfaces. 


Points  for  large  fin. 


Fin  solid  from  points. 


Rotate  fin  solid  to  create 
the  other  three  large  fins. 
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42 
23 

4.63  0  0 
TRAN 

43 
42 

6  0  0 

SOLI 

PARA 

8 

42  43 
ROTA 
9 
8 

270  0  0 
3 


Translate  surfaces  to  be 
used  for  small  fins. 


Create  small  fin  with  a 
parametric  fit  between 
two  surfaces.  Rotate  solid 
to  create  other  fins. 


GLOB 

GROU 

NEW 

ADD 

SOLI 

123456789  10  11 

SAVE 
MAV 
1  0  1 


Place  all  solids  into  a  group. 


DISP  ENT  ALL  OFF 

DISP  ENT  POST  GALL  ON  Set  . 

DISP  ENT  GRQU  ON  display  options  to  display 

DRAW  groups  only. 

EXIT 

Y 

MAV. DAT 

Y 
N 
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The  following  commands  and  procedure  may  be  used  to  create  and  evaluate  a 
physical-fit  compatability  of  a  weapon  on  an  aircraft. 

Begin  by  executing  DBMERGE  to  get  the  following  groups  from  their  original 
files  to  create  the  file  tcasel.dat. 


GROUP 


FILE 


tclving 

sweep 

leinbpyl 

runway 

mk84 


wing. geo 
wingswp.geo 
inbpylon. geo 
runway . geo 
mk84.geo 


Now  execute  MOVIESTAR,  open  the  file  teasel  .dat  and  execute  the  following 
commands. 


DISP  ENT  ALL  OFF 
DISP  ENT  GROU  ON 
GLOB  COLO 
GROU 
TC1WING 
1  0  0 
GROU 
SWEEP 
0  1  0 
GROU 

LEINBPYL 
1  1  0 
GROU 
MK84 
0  0  1 
GROU 
RUNWAY 
1  1  1 


Set  display-options  to 
display  only  groups  and 
give  each  group  a  different 
color  so  they  will  be  easy 
to  see. 


VIEW 

ROT  Orient  the  model  to  a 
SCREE  side-on  configuration. 
-90  0  0 
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DRAW 


GLQB 
GROU 
ROTA 
MK84 
MK84 
45  0  180 


Orient  the  weapon  so  it  can 
be  placed  on  the  pylon. 


LIST 

GROU 

MK84 


GROU 

LENINBPYL  Get  hook  point  numbers . 

GROU 

LEINBPYL 


GROU 

SNAP 

LEINBPYL 

6088 

MK84 

6255 

DRAW 


After  weapon  is  properly 
oriented  we  can  use  the 
group  snap  command  to  place 
it  in  the  proper  position. 


TRANS 
MK84 
MK84 
0  0  8 
NEW 
ADD 
SOLI 
53 

SAVE 
MID 
0  0  1 
DRAW 


Now  we'll  translate  it  up 
so  it  interferes  with  the 
pylon. 


Create  a  new  group  of  one  of 
the  weapon  solids  for  quick 
interf erence  calculation. 
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CALIPER 


Execute  caliper. 


MID 

LEINBPYL 

Y 

0  0  1 


HIGH 
10  1 

Y 

EXIT 
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Y 
N 


Display  interference 
highlight . 


103 


Optimal  packaging  examples  with  PACKER. 

Several  different  simple  examples  Me  given  here  of  packaging  objects  to  user 
defined  objective  functions  and  constraints. 

The  problems  given  here  have  only  a  few  objects  and  a  few  container  walls 
for  simplicity.  More  complex  problems  have  been  solved  and  documented.  These 
problems  here  Me  for  demonstration  and  tutorial  and  do  not  reflect  the  type  of 
complex  problems  that  can  be  solved. 

The  first  example  is  one  where  there  is  a  single  container  wall  and  two  objects 
and  we  wish  to  minimize  the  x  centroid  of  the  composite  model.  We  setup  to  only 
translate,  and  put  constraints  on  the  problem  so  that  there  will  be  a  sepMation 
between  the  objects  and  between  the  objects  and  the  container  wall  of  .1  inches. 
The  example'  shows  the  input  (and  some  output)  through  the  complete  problem. 

The  second  example  takes  the  model  from  the  first  and  simply  changes  the 
objective  function  to  minimize  Ixx  with  an  equality  constraint  on  the  x  centroid 
(=  1.2).  The  solution  shows  how  PACKER  can  find  a  feasible  solution  when  the 
starting  design  is  not  feasible. 

The  last  example  is  a  simple  one  with  rotational  degrees  of  freedom.  Here  an 
object  (cylinder)  is  oriented  Mbitrarily  at  the  starting  design.  The  objective  is  to 
minimize  /„,  only  now  the  rotations  will  be  found  that  orient  the  object  in  such  a 
way  that  the  object’s  axis  of  symmetry  is  lined  up  with  the  globed  axis  (the  optima] 
solution  for  a  cylinder). 


104 


We  begin  with  a  NAVGRAPH  session.  The  following  command  create  the  desired 
model  shown. 
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2 
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10 
9 

315 
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ENTE 

13 

25  26  27  28 

TRAN 

14 
13 

0  0  2 

SOLI 

PARA 

3 

13  14 


a  box. 


GLOB 

GROU 

NEW 

ADD 

SOLI 

1 


SAVE 

CONT 

1 

NEW 

ADD 

SOL 

2 


define  the  container  group. 


SAVE 

CYL1 

®  *  define  the  cylinder  group. 

NEW 

ADD 

SOLI 

3 


SAVE 
B0X1 
0  0  1 
DR 


define  the  hex  group. 


VIEW  ROTA  SCREE  0  -90 
DRAW 
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GROUP 
TRAN 
CYL1 
CYL1 
0-2  2 


translate  the  cylinder. 


DRAW 

GENE 

CREA 

ELEN 

SOLI 

1 


mash  tha  container  with  one 
hexahedron  finite  element  for 
mass  properties  calculation. 


Ill 


SI 

SOLI 

2  mash  the  cylinder. 


10  10  1  10  elements  along  the  curved 

sections  of  the  cylinder. 

S2 

SOLI 

3  mesh  the  box. 


Ill  one  hexahedron  is  enough  for  a  box. 

S3 

DRAW 

EXIT 


108 


The  model  is  created  and  saved,  we  call  it  PACK.DAT.  Now  we  run  MASSPROP 
and  calculate  the  mass  properties  for  each  object.  We  need  the  volume,  the  density,  the 
center  of  gravity  and  the  local  (centroidal)  moments  of  inertial  for  each  object.  The  output 
files  from  MASSPROP  for  each  group  is  shown  here.  The  arrow  shows  the  information 
needed  for  input  into  the  filename.MAS  (PACK.MAS  in  this  case). 


NAVGRAPH  DATA  FILE  ■  PACK 
GROUP  NAME  -  CONT 

#  OF  HEXAHEDRON  ELEMENTS  -  1 

#  OF  HEDGE  ELEMENTS  -  0 

#  OF  PENTAHEDRON  ELEMENTS  *  0 

#  OF  TETRAHEDRON  ELEMENTS  ■  0 


SPECIFIC  HEIGHT 
GRAVITY  CONST. 
DENSITY 
TOTAL  VOLUME 
TOTAL  MASS 
GLOBAL  XBAR 
GLOBAL  YBAR 
GLOBAL  ZBAR 
GLOBAL  IXX 
GLOBAL  IYY 
GLOBAL  IZZ 
GLOBAL  IXY 
GLOBAL  IXZ 
GLOBAL  IYZ 
CENTROID  IXX 
CENTROID  IYY 
CENTROID  IZZ 
CENTROID  IXY 
CENTROID  IXZ 
CENTROID  IYZ 


490.0000 
32.20000 
15.21739 
64.00000 
973.9130 
0.5000000 
6 . 2670026E-08 
9 . 4005038E-08 
10388.41 
5518.840 
5518.840 
-3.0517578E-0S 
-3 . 0517S78E-05 
-6 . 1035156E-05 
10388.41 
5275.362 
5275 . 362 
-3.6379788E-12 
1 . 525878SE-05 
-6 . 1035149E-05 


< 

< 

< 

< 

< 


< 

< 

< 

< 

< 

< 


NAVGRAPH  DATA  FILE  -  PACK 
GROUP  NAME  -  CYL1 

#  OF  HEXAHEDRON  ELEMENTS  «  100 

i  OF  HEDGE  ELEMENTS  *  0 

i  OF  PENTAHEDRON  ELEMENTS  *  0 

•  OF  TETRAHEDRON  ELEMENTS  -  0 


SPECIFIC  HEIGHT 
GRAVITY  CONST. 
DENSITY 


490.0000 

32.20000 

15.21739 
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< 


TOTAL  VOLUME 
TOTAL  MASS 
GLOBAL  XBAR 
GLOBAL  YBAR 
GLOBAL  ZBAR 
GLOBAL  IXX 
GLOBAL  IYY 
GLOBAL  IZZ 
GLOBAL  IXY 
GLOBAL  IXZ 
GLOBAL  IYZ 
CENTROID  IXX 
CENTROID  IYY 
CENTROID  IZZ 
CENTROID  IXY 
CENTROID  IXZ 
CENTROID  IYZ 


9.386064 
142.8314 
5.500000 
-2.000000 
2.000000 
1213.773 
5034.661 
5034.661 
1571.146 
-1571.145 
571.3256 
71.12195 
142.6846 
142.6846 
O.OOOOOOOE+OO 
8.5449219E-04 
-1 . 8310547E-04 


NAVGRAPH  DATA  FILE  -  PACK 
GROUP  NAME  -  B0X1 

#  OF  HEXAHEDRON  ELEMENTS  » 
t  OF  WEDGE  ELEMENTS 

f  OF  PENTAHEDRON  ELEMENTS  - 

#  OF  TETRAHEDRON  ELEMENTS  * 


SPECIFIC  WEIGHT 
GRAVITY  CONST. 
DENSITY 
TOTAL  VOLUME 
TOTAL  MASS 
GLOBAL  XBAR 
GLOBAL  YBAR 
GLOBAL  ZBAR 
GLOBAL  IXX 
GLOBAL  IYY 
GLOBAL  IZZ 
GLOBAL  IXY 
GLOBAL  IXZ 
GLOBAL  IYZ 
CENTROID  IXX 
CENTROID  IYY 
CEN1 KOID  IZZ 
CENTROID  IXY 
CENTROID  IXZ 
CENTROID  IYZ 


490.0000 
32.20000 
15.21739 
4.000000 
60.86956 
5.000000 
2 . 500000 
0.9999999 
466.6667 
1623.188 
1927.536 
-760.8696 
-304.3478 
-152.1739 
25.36227 
40.57971 
25.36230 
O.OOOOOOOE+OO 
-3.0517578E-05 
O.OOOOOOOE+OO 


These  data  are  put  into  the  file  PACK. MAS  in  free  format  in  the  following  manner. 
Each  object  has  two  lines,  the  first  line  has  the  volume,  density,  x  centroid,  y  centroid 
and  z  centroid.  The  second  line  has  the  six  moments  of  inertia  in  the  order  above.  The 
PACK.MAS  file  is  shown  here  (notice  the  values  from  above). 

64.0  15.21739  .5  0.  0. 

10388.41  5275.362  5275.362  0.  0.  0. 

9.386064  15.21739  5.5  -2.0  2.0 
71.12195  142.6846  142.6846  0.  0.  0. 

4.0  15.21739  5.  2.5  1.0 
25.36227  40.57971  25.3623  0.  0.  0. 

We  are  now  ready  to  run  PREPAC.  PREPAC  calculates  all  of  the  initial  analsis 
functions  from  the  initial  analysis  variables.  The  user  is  asked  for  the  NAVGRAPH  file 
name,  the  container  name,  the  mass  properties  file  name  and  a  location  for  the  center  of 
pressure  (optional).  The  user  is  also  asked  what  the  OPTDES  problem  file  name  is.  In 
this  problem  we  choose  PACK. PRO. 

SETUP  is  the  next  step.  Here  we  define  the  design  variables  (tranlations  and  rota¬ 
tions)  as  a  subset  of  the  analysis  variables  and  the  design  function  (the  objective  and  the 
constraints)  as  a  subset  of  the  analysis  functions  (separation  distances  and  mass  proper¬ 
ties). 


The  following  is  what  is  seen  on  the  screen  in  the  SETUP  session. 


WELCOME  TO  THE  OPTDES. B YU  VERSION  4.0  SETUP  PROGRAM 
COPYRIGHT  (C)  1988  BY  BRIGHAM  YOUNG  UNIVERSITY 

PROMPT  TYPES: 

"<...>"  TYPE  RESPONSE 

PRESS  RETURN  KEY 

SPECIAL  CHARACTERS: 

ESCAPE  TO  MAIN  MENU 
"X"  EXECUTE  OPERATING  SYSTEM  COMMAND 
OPEN  AND  CLOSE  INTERACTION  FILES 

INTERACTION  FILE  NAMES: 

"LOG"  STARTING  LOG  FILE  NAME 
"RES”  STARTING  RESPONSE  FILE  NAME 
"NUL"  THERE  IS  NO  STARTING  BATCH  FILE 

<NAME  OF  PROBLEM  FILE>  PACK. PRO  input  the  problem  name 

SETUP  MENU : D IS, VAR, FUN, COM ,QUI, HELP 
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DIS 

DISPLAY  MENU :AV.AF. HELP 


<SELECT>  AV 

AV# 

VALUE 

1 

0 . OOOOOE+OO 

2 

0 . OOOOOE+OO 

3 

0. OOOOOE+OO 

4 

0. OOOOOE+OO 

5 

0. OOOOOE+OO 

6 

0. OOOOOE+OO 

7 

0. OOOOOE+OO 

8 

0. OOOOOE+OO 

9 

0. OOOOOE+OO 

10 

0. OOOOOE+OO 

11 

0. OOOOOE+OO 

12 

0. OOOOOE+OO 

DESCRIPTION 

1QUAT, 

GROUP 

CYL1 

2QUAT, 

GROUP 

CYL1 

3QUAT, 

GROUP 

CYL1 

1TRAN, 

GROUP 

CYL1 

2TRAN, 

GROUP 

CYL1 

3TRAN , 

GROUP 

CYL1 

1QUAT , 

GROUP 

B0X1 

2QUAT , 

GROUP 

B0X1 

3QUAT , 

GROUP 

B0X1 

1TRAN, 

GROUP 

B0X1 

2TRAN, 

GROUP 

B0X1 

3TRAN, 

GROUP 

B0X1 

display  analysis 
variables 


SETUP  MENUiDIS, VAR, FUN, COM, QUI, HELP 

<SELECT>  VAR  edit  design  variables 

EDIT  DESIGN  VARIABLES 

Cl)  (2)  (3)  (4)  (5)  (6) 

DV#  #MP  AV#  COEFF  VALUE  MINIMUM  MAXIMUM  DESCRIPTION 


1) 

<EDIT:  C,G,J,Q,HELP>  S  1  3  2  4  1  G  4  6  2  10  i 

select  design  variables  as 
all  tranlations.  Note,  that 
by  not  defining  any  rotations 
the  objects  will  not  rotate. 


EDIT 

DESIGN 

VARIABLES 

(1) 

(2) 

(3) 

(4)  (5) 

(6) 

DV# 

#MP 

AV# 

COEFF 

VALUE  MINIMUM 

MAXIMUM 

DESCRIPTION 

1) 

1 

1 

4 

1.00 

0 . OOE+OOP  -1.00 

1.00 

1TRAN, 

GROUP 

CYL1 

2) 

2 

1 

5 

1.00 

0 . OOE+OOP  -1.00 

1.00 

2TRAN, 

GROUP 

CYL1 

3) 

3 

1 

6 

1.00 

0. OOE+OOP  -1.00 

1.00 

3TRAN, 

GROUP 

CYL1 

4) 

4 

1 

10 

1.00 

0. OOE+OOP  -1.00 

1.00 

1TRAN, 

GROUP 

B0X1 

6) 

5 

1 

11 

1.00 

0. OOE+OOP  -1.00 

1.00 

2TRAN , 

GROUP 

B0X1 

6) 

6 

1 

12 

1.00 

0. OOE+OOP  -1.00 

1.00 

3TRAN , 

GROUP 

B0X1 

7) 

<EDIT 

:  C,G,J,Q,HELP>  G  1 

6  5  -5  G  1  6  6  S 

edit 

how  far 

the  obj 

ects 
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can  translate  (min/max) 


EDIT  DESIGN  VARIABLES 


(1) 

(2) 

(3) 

(4)  (5) 

(6) 

DV# 

#HP 

AV# 

COEFF 

VALUE  MINIMUM 

MAXIMUM 

DESCRIPTION 

1) 

1 

1 

4 

1.00 

0 . OOE+OOP  -5.00 

5.00 

1TRAN,  GROUP 

CYLI 

2) 

2 

1 

S 

1.00 

O.OOE+OOP  -5.00 

5.00 

2TRAN,  GROUP 

CYL1 

3) 

3 

1 

6 

1.00 

0 . OOE+OOP  -5.00 

5.00 

3 TRAN,  GROUP 

CYL1 

4) 

4 

1 

10 

1.00 

O.OOE+OOP  -5.00 

5.00 

1TRAN,  GROUP 

B0X1 

5) 

5 

1 

11 

1.00 

0. OOE+OOP  -5.00 

5.00 

2TRAN,  GROUP 

B0X1 

6) 

6 

1 

12 

1.00 

O.OOE+OOP  -5.00 

5.00 

3TRAN,  GROUP 

B0X1 

7) 

<EDIT:  C,G,J,Q,HELP>  Q 


SETUP  MENU:DIS, VAR, FUN.COM, QUI, HELP 
<SELECT>  DIS  AF  FUN 


AF# 

VALUE 

DESCRIPTION 

1 

3.0000 

ICON  SEP  ft  CYLi 

2 

3.0000 

ICON  SEP  ft  B0X1 

3 

3.0000 

CYLI  0-0  SEP  B0X1 

4 

5.7446 

ICON  CEN  ft  CYLI 

5 

5.2440 

ICON  CEN  ft  B0X1 

6 

4.6368 

CYLI  0-0  CEN  B0X1 

7 

77.386 

TOTAL  VOLUME 

8 

1177  .6 

TOTAL  MASS 

9 

1.3390 

X  MASS  CENTROID 

10 

-0.11336 

Y  MASS  CENTROID 

11 

0.29427 

Z  MASS  CENTROID 

12 

11952. 

XX  MOMENT  OF  INERTIA 

13 

9963.2 

YY  MOMENT  OF  INERTIA 

14 

10354. 

ZZ  MOMENT  OF  INERTIA 

15 

631.53 

XY  MOMENT  OF  INERTIA 

16 

-1411.5 

XZ  MOMENT  OF  INERTIA 

17 

379.87 

YZ  MOMENT  OF  INERTIA 

18 

4.0000 

X  PRESSURE  CENTER 

19 

0 . OOOOOE+OO 

Y  PRESSURE  CENTER 

20 

0 . OOOOOE+OO 

Z  PRESSURE  CENTER 

display  analysis 
functions 


EDIT  DESIGN  FUNCTIONS 
(1)  (2)  (3) 

DF#  #MP  AF#  COEFF 


1) 


(4)  (5)  (6) 

VALUE  ALLOWABLE  RANGE  DESCRIPTION 
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<EDIT:  C,G, J,Q,HELP>  C129C14v  max  xbar 

as  objective 

EDIT  DESIGH  FUNCTIONS 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

DF«  #MP 

AF* 

COEFF 

VALUE 

ALLOWABLE 

RANGE 

DESCRIPTION 

1)  1  1 

9 

1.00 

1.3 

vP 

0 . OOOE+OO 

1.00 

X  MASS  CENTROID 

2) 

<EDIT:  C,G,J,Q,HELP>  G  2  4  2  1  1  select  contraints 

as  separations 

EDIT  DESIGN  FUNCTIONS 


DF# 

(1) 

#MP 

(2) 

AF# 

(3) 

COEFF 

VALUE 

(4) 

(5) 

ALLOWABLE 

(6) 

RANGE 

DESCRIPTION 

1) 

1 

1 

9 

1.00 

1.3 

vP 

0. OOOE+OO 

1.00 

X  MASS  CENTROID 

2) 

2 

1 

1 

1.00 

3.0 

<P 

0. OOOE+OO 

1.00 

ICON  SEP  ft  CYL1 

3) 

3 

1 

2 

1.00 

3.0 

<P 

0. OOOE+OO 

1.00 

ICON  SEP  ft  B0X1 

4) 

5) 

4 

1 

3 

1.00 

3.0 

<P 

0. OOOE+OO 

1.00 

CYL1  0-0  SEP  B0X1 

<EDIT:  C,G, J,q,HELP>  O  2  4  4  >  G  2  4  5  .1  make  them  inequality 

greater  than  .1 

EDIT  DESIGN  FUNCTIONS 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

DF# 

#MP 

AF# 

COEFF 

VALUE 

ALLOWABLE 

RANGE 

DESCRIPTION 

1) 

1 

1 

9 

1.00 

1.3 

vP 

0. OOOE+OO 

1.00 

X  MASS  CENTROID 

2) 

2 

1 

1 

1.00 

3.0 

>P 

0.100 

1.00 

ICON  SEP  ft  CYL1 

3) 

3 

1 

2 

1.00 

3.0 

>P 

0.100 

1.00 

ICON  SEP  ft  B0X1 

4) 

4 

1 

3 

1.00 

3.0 

>P 

0.100 

1.00 

CYL1  0-0  SEP  B0X1 

S) 

<EDIT :  C,G,J,Q,HELP>  Q 

SETUP  MENU :DIS, VAR, FUN, COM, QUI, HELP 

<SELECT>  qUI  «e  are  ready  to 

optimize  with  PACKER 


<ST0RE  THIS  PR0BLEM?>  Y 
<NAME  OF  PROBLEM  FILE>  PACK. PRO 
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We  can  now  run  PACKER  for  the  optimization  session. 

The  following  is  what  is  seen  on  the  screen  during  the  PACKER  session. 


WELCOME  TO  THE  OPTDES.BYU  VERSION  4.0  DESIGN  PROGRAM 
COPYRIGHT  (C)  1988  BY  BRIGHAM  YOUNG  UNIVERSITY 

PROMPT  TYPES: 

••<...>*•  TYPE  RESPONSE 
"{...}"  PRESS  RETURN  KEY 

SPECIAL  CHARACTERS: 

ESCAPE  TO  MAIN  MENU 
"V  EXECUTE  OPERATING  SYSTEM  COMMAND 
OPEN-AHD  CLOSE  INTERACTION  FILES 

INTERACTION  FILE  NAMES: 

"LOG"  STARTING  LOG  FILE  NAME 
"RES"  STARTING  RESPONSE  FILE  NAME 
"NUL"  THERE  IS  N0  STARTING  BATCH  FILE 

<NAME  OF  PROBLEM  FILE>  PACK. PRO 

****  A  PACKER  BANNER  IS  SHOWN  HERE**** 

The  user  is  asked  for  the  NAVGRAPH  file,  the  container  name,  the 
mass  properties  file,  and  the  center  of  pressure  (optional)  again. 

DESIGN  MENU :  PRO , DIS . SET , RUN , PAR , HIS ,  EXP , TOL .APP.COM, QUI . HELP 
<SELECT>  PAR  GRG 


EDIT  GENERALIZED  REDUCED  GRADIENT  PARAMETERS 
VALUE  DESCRIPTION 


( 

1) 

0 

( 

2) 

0. 100000E-02 

( 

3) 

0 . 100000E-04 

( 

4) 

0. 200000E-01 

( 

5) 

0.100000 

( 

6) 

0. 100000E-02 

ALGORITHM  PRINTOUT  (1-ON  O-OFF) 
CONVERGENCE  TOLERANCE 
MINIMUM  STEPSIZE 
INITIAL  STEPSIZE 
MAX  BOUND  VIOLATION 
CONSTRAINT  TOLERANCE 


<EDIT:  C,G, J,Q,HELP>  C  1  1  Q  edit  to  have  GRG  print 

out  more  info  to  the  screen 


DESIGN  MENU : PRO . D IS , SET , RUN , P  AR , HIS , EXP ,T0L , APP , COM , QUI , HELP 
<SELECT>  PAR  GRA 
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EDIT  GRADIENT  PARAMETERS 

VALUE  DESCRIPTION 


(  l)  1  GRADIENT  TYPE  (1-FINITE  DIFF, 2-EXPLICIT) 

(  2)  0. 100000E-03  SCALED  PERTURBATION  FOR  F-D  GRADIENTS 


<EDIT:  C,G,J,Q,HELP>  C  1  2  Q 


adit  for  explicit 
gradients  to  be  on 


DESIGN  MENU:PRO,DIS, SET, RUN, PAR.HIS, EXP, TOL, APP.COM, QUI, HELP 
<SELECT>  PRO  POS  post  process,  this  gives  a 

NAVGRAPH  menu  for  graphics. 


DESIGN  MENU : PRO , DIS , SET , RUN , PAR , HIS . EXP ,TOL , APP , COM . QUI , HELP 
<SELECT>  RUN  GRG  3  run  3  GRG  iterations 

OBJECTIVE  FUNCTION  GRADIENT 

0.60644  0 . OOOOOE+OO  0. OOOOOE+OO  0.2S844  0. OOOOOE+OO  0 . OOOOOE+OO 

LAGRANGIAN  FUNCTION  GRADIENT 

0 . OOOOOE+OO  0. OOOOOE+OO  0. OOOOOE+OO  0 . OOOOOE+OO  0 . OOOOOE+OO  0. OOOOOE+OO 

LENGTH  OF  SEARCH  DIRECTION  VECTOR-  0.65173 
SEARCH  DIRECTION 

-0.88994  0 . 13808E-05-0 . 17320E-05-0 . 45607  0 . OOOOOE+OO  0 . OOOOOE+OO 

LINE  SEARCH  STEP  1  STEP  LENGTH-  2.0000E-02 
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LINE  SEARCH  STEP  2  STEP  LENGTH-  4.0000E-02 

LINE  SEARCH  STEP  3  STEP  LENGTH-  8.0000E-02 

LINE  SEARCH  STEP  4  STEP  LENGTH-  1.6000E-01 

LINE  SEARCH  STEP  S  STEP  LENGTH-  3.2000E-01 

LINE  SEARCH  STEP  6  STEP  LENGTH-  5.7000E-01 

LINE  SEARCH  STEP  7  STEP  LENGTH-  8.2000E-01 

CONSTRAINT  1  NOT  FORMERLY  BINDING 
NEVTON-RAPHSON  ITERATION 
CONVERGENCE  IN  NEWTON-RAPHSON 

DESIGN  NUMBER  1 

SCALED  STEP  LENGTH  -  0.6517281 

SCALED  OBJECTIVE  (DF*  1)  »  0.9104874 

SCALED  MAX  INEQUALITY  (DF#  2)  -  0 . 1025966E-05 
ANALYSIS  CALLS  -  9  GRADIENT  CALLS  =  1 

UPDATING  HESSIAN 
OBJECTIVE  FUNCTION  GRADIENT 

0.60644  0 . OOOOOE+OO  O.OOOOOE+OO  0.25844  O.OOOOOE+OO  O.OOOOOE+OO 

LAGRANGIAN  FUNCTION  GRADIENT 

0.60644  O.OOOOOE+OO  O.OOOOOE+OO  0.25844  O.OOOOOE+OO  O.OOOOOE+OO 

LENGTH  OF  SEARCH  DIRECTION  VECTOR-  0.28276 
SEARCH  DIRECTION 

0.55326E-06  0.63876E-03  0.27091E-03  -1.0000  0.35697E-04  0.26065E-04 

LINE  SEARCH  STEP  1  STEP  LENGTH-  2.0000E-02 
LINE  SEARCH  STEP  2  STEP  LENGTH-  4.0000E-02 
LINE  SEARCH  STEP  3  STEP  LENGTH-  8.0000E-02  again. 

LINE  SEARCH  STEP  4  STEP  LENGTH-  1.6000E-01 

LINE  SEARCH  STEP  5  STEP  LENGTH-  3.2000E-01  | 

CONSTRAINT  2  NOT  FORMERLY  BINDING 
NEWTON-RAPHSON  ITERATION 
CONVERGENCE  IN  NEWTON-RAPHSON 

DESIGN  NUMBER  2 

SCALED  STEP  LENGTH  -  0.2827642  , 

SCALED  OBJECTIVE  (DF#  1)  =  0.8374085 

SCALED  MAX  INEQUALITY  (DF*  3)  -  0 . 3619690E-06  ’ 

ANALYSIS  CALLS  -  16  GRADIENT  CALLS  -  2  ! 

UPDATING  HESSIAN  I 

OBJECTIVE  FUNCTION  GRADIENT  ! 

0.60644  O.OOOOOE+OO  O.OOOOOE+OO  0.25844  O.OOOOOE+OO  O.OOOOOE+OO  | 

i 

LAGRANGIAN  FUNCTION  GRADIENT 

0.60644  O.OOOOOE+OO  O.OOOOOE+OO  0.25844  O.OOOOOE+OO  O.OOOOOE+OO 
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LENGTH  OF  SEARCH  DIRECTION  VECTOR-  0.00377 
SEARCH  DIRECTION 

-0.51617E-02  0.73796  0.64065  -0.55084E-03-0.57649E-01  0.20404 

LINE  SEARCH  STEP  1  STEP  LENGTH-  3.7690E-03 

LINE  SEARCH  STEP  2  STEP  LENGTH-  7.S380E-03 

LINE  SEARCH  STEP  3  STEP  LENGTH-  1.5076E-02  again. 

LINE  SEARCH  STEP  4  STEP  LENGTH-  3.0152E-02 

LINE  SEARCH  STEP  5  STEP  LENGTH-  6.0304E-02 

NEVTON-RAPHSON  ITERATION 
CONVERGENCE  IN  NEVTON-RAPHSON 
FITTING  PARABOLA 

LINE  SEARCH  STEP  6  STEP  LENGTH-  3.3913E-02 
FITTING  PARABOLA 

LINE  SEARCH  STEP  7  STEP  LENGTH-  2.4549E-02 
FITTING  PARABOLA 

LINE  SEARCH  STEP  8  STEP  LENGTH-  2.8872E-02 
FITTING  PARABOLA 

LINE  SEARCH  STEP  9  STEP  LENGTH-  3.1080E-02 
DESIGN  NUMBER  3 

SCALED  STEP  LENGTH  -  0.3015201E-01 

SCALED  OBJECTIVE  (DF#  1)  -  0.8373099 

SCALED  MAX  INEQUALITY  (DF#  2)  -  0.7766739E-03 
ANALYSIS  CALLS  -  27  GRADIENT  CALLS  -  3 


DESIGN  MENU : PRO , DIS . SET , RUN , PAR , HIS , EXP , TOL , APP , COM , QUI , HELP 
<SELECT>  RUN  GRG  1  run  one  more  GRG  run. 

OBJECTIVE  FUNCTION  GRADIENT 

0.60644  0 . OOOOOE+OO  O.OOOOOE+OO  0.25844  O.OOOOOE+OO  O.OOOOOE+OO 

LAGRANGIAN  FUNCTION  GRADIENT 

0.60644  O.OOOOOE+OO  O.OOOOOE+OO  0.25844  O.OOOOOE+OO  O.OOOOOE+OO 

LENGTH  OF  SEARCH  DIRECTION  VECTOR-  0.00524 
SEARCH  DIRECTION 

0.22105E-01  0 . 89388E-01  0.99486  0.31673E-02  0.11111E-01  0.40634E-01 

LINE  SEARCH  STEP  1  STEP  LENGTH-  5.2357E-03 

DESIGN  NUMBER  4 

SCALED  STEP  LENGTH  -  0 . OOOOOOOE+OO 

SCALED  OBJECTIVE  (DF#  1)  -  0.8373099 

SCALED  MAX  INEQUALITY  (DF*  2)  -  0.7766739E-03 


118 


ANALYSIS  CALLS  -  29  GRADIENT  CALLS  =*  4 

OPTIMUM  ASSUMED  (MINIMUM  STEP SIZE  REACHED)  ve  reached  an 

optimal  solution 

DESIGN  MENU : PRO , DIS , SET , RUN ,PAR , HIS , EXP , TOL , APP , COM , QUI , HELP 
<SELECT>  PRO  POS  draw  it  with  NAVGRAPH  to 

visually  verify 


DESIGN  MENU.-PRO, DIS, SET, RUN, PAR, HIS, EXP, TOL,  APP, COM,  QUI, HELP 
<SELECT>  QUI 
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The  second  example  starts  with  the  same  model  as  the  first.  We  start  with  the 
same  initial  design,  therefore,  we  do  not  need  to  remodel  or  recalculate  the  mass  propeties. 
We  simply  run  SETUP  again  and  redefine  the  objective  and  contraints,  and  maybe  add 
more  design  variables.  In  this  example,  we  will  minimize  Ixx  while  containing  the  global 
composite  x  centroid  to  be  1.2. 


WELCOME  TO  THE  OPTDES.BYU  VERSION  4.0  SETUP  PROGRAM 
COPYRIGHT  (C)  1988  BY  BRIGHAM  YOUNG  UNIVERSITY 

PROMPT  TYPES: 

"<...>"  TYPE  RESPONSE 

PRESS  RETURN  KEY 

SPECIAL  CHARACTERS: 

ESCAPE  TO  MAIN  MENU 
"X"  EXECUTE  OPERATING  SYSTEM  COMMAND 
"«"  OPEN  AND  CLOSE  INTERACTION  FILES 


INTERACTION  FILE  NAMES: 

"LOG'1  STARTING  LOG  FILE  NAME 
"RES"  STARTING  RESPONSE  FILE  NAME 
"NUL"  THERE  IS  NO  STARTING  BATCH  FILE 


<NAME  OF  PROBLEM  FILE>  pack. pro 


start  with  the  previous 
starting  design. 


SETUP  MENU:DIS, VAR, FUN.COM, QUI, HELP 
<SELECT>  dis  af  fun 


AF# 

VALUE 

DESCRIPTION 

1 

3.0000 

ICON  SEP  k  CYL1 

2 

3.0000 

ICON  SEP  k  B0X1 

3 

3.0000 

CYL1  0-0  SEP  B0X1 

4 

5.7446 

ICON  CEN  k  CYL1 

5 

5.2440 

ICON  CEN  k  B0X1 

6 

4.6368 

CYL1  0-0  CEN  B0X1 

7 

77.386 

TOTAL  VOLUME 

8 

1177.6 

TOTAL  MASS 

9 

1.3390 

X  MASS  CENTROID 

10 

-0.11336 

Y  MASS  CENTROID 

11 

0.29427 

Z  MASS  CENTROID 

12 

11952. 

XX  MOMENT  OF  INERTIA 

13 

9963.2 

YY  MOMENT  OF  INERTIA 

14 

10354. 

ZZ  MOMENT  OF  INERTIA 

15 

631.53 

XY  MOMENT  OF  INERTIA 
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16 

-1411.  S 

XZ 

MOMENT  OF 

INERTIA 

17 

379.87 

YZ 

MOMENT  OF 

INERTIA 

18 

4.0000 

X 

PRESSURE  CENTER 

19 

0 . OOOOOE+OO 

Y  PRESSURE 

CENTER 

20 

0 . 00000E+00 

Z 

PRESSURE 

CENTER 

EDIT  DESIGN  FUNCTIONS 


DF# 

(1) 

•MP 

(2) 

AF# 

(3) 

COEFF 

VALUE 

(4) 

(5) 

ALLOWABLE 

(6) 

RANGE 

DESCRIPTION 

1) 

1 

1 

9 

1.00 

1.3 

vP 

0 .  OOOE+OO 

1.00 

X  MASS  CENTROID 

2) 

2 

1 

1 

1.00 

3.0 

>P 

0.100 

1.00 

ICON  SEP  k  CYL1 

3) 

3 

1 

2 

1.00 

3.0 

>P 

0.100 

1.00 

ICON  SEP  ft  B0X1 

4) 

5) 

4 

1 

3 

1.00 

3.0 

>P 

0.100 

1.00 

CYL1  0-0  SEP  B0X1 

<EDIT:  C,G,J,q,HELP>  c  1  2  12  change  the  objective 

function  to  minimize  Izz 

EDIT  DESIGN  FUNCTIONS 


DFt 

(1) 

•MP 

(2) 

AF# 

(3) 

COEFF 

(4) 

VALUE 

(S) 

ALLOWABLE 

(6) 

RANGE 

DESCRIPTION 

1) 

1 

1 

12 

1.00 

0. 12E+05vP 

0. OOOE+OO 

1.00 

XX  MOMENT  OF 

INERT 

2) 

2 

1 

1 

1.00 

3.0 

>P 

0.100 

1.00 

ICON 

SEP  * 

CYL1 

3) 

3 

1 

2 

1.00 

3.0 

>P 

0.100 

1.00 

ICON 

SEP  k 

B0X1 

4) 

5) 

4 

1 

3 

1.00 

3.0 

>P 

0.100 

1.00 

CYL1 

0-0  SEF 

B0X1 

<EDIT :  C,G, J,Q,HELP>  c  5  2  9  c  5  4  >  c  5  5  1.2 

add  equality 
constraint  on 
z  centroid 

EDIT  DESIGN  FUNCTIONS 


DF# 

(1) 

•MP 

(2) 

AF# 

(3) 

COEFF 

VALUE 

(4) 

(S) 

ALLOWABLE 

(6) 

RANGE 

DESCRIPTION 

1) 

1 

1 

12 

1.00 

0 . 12E+05vP 

0. OOOE+OO 

1.00 

XX  MOMENT  OF  INERTIA 

2) 

2 

1 

1 

1.00 

3.0 

>P 

0.100 

1.00 

ICON  SEP  k  CYL1 

3) 

3 

1 

2 

1.00 

3.0 

>P 

0.100 

1.00 

ICON  SEP  k  B0X1 

4) 

4 

1 

3 

1.00 

3.0 

>P 

0.100 

1.00 

CYL1  0-0  SEP  B0X1 

5) 

6) 

5 

1 

9 

1.00 

1.3 

>P 

1.20 

1.00 

X  MASS  CENTROID 

<EDIT:  C,G,J,Q,HELP>  q 
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SETUP  MENU:DIS. VAR, FUN, COM, QUI, HELP 
<SEL£CT>  qui 

<STORE  THIS  PROBLEM? >  y 

<NAME  OF  PROBLEM  FILE>  pack2.pro  call  the  problem  file 

something  new. 


Now  for  the  optimization  in  PACKER. 


WELCOME  TO  THE  OPTDES.BYU  VERSION  4.0  DESIGN  PROGRAM 
COPYRIGHT  (C)  1988  BY  BRIGHAM  YOUNG  UNIVERSITY 

PROMPT  TYPES: 

"<• . .>"  TYPE  RESPONSE 

PRESS  RETURN  KEY 

SPECIAL  CHARACTERS: 

"*"  ESCAPE  TO  MAIN  MENU 

"%"  EXECUTE  OPERATING  SYSTEM  COMMAND 

"•"  OPEN  AND  CLOSE  INTERACTION  FILES 

INTERACTION  FILE  NAMES: 

"LOG"  STARTING  LOG  FILE  NAME 
"RES"  STARTING  RESPONSE  FILE  NAME 
"NUL"  THERE  IS  NO  STARTING  BATCH  FILE 

<NAME  OF  PROBLEM  FILE>  PACK2 . PRO 


****THE  PACKER  BANNER  HERE*****  and  the  user  input  PACK . DAT , 
the  container  name,  PACK. MAS  again. 


DESIGN  NUMBER  0 


SCALED 

OBJECTIVE 

(DF« 

1)  - 

11951.75 

SCALED 

MAX  INEQUALITY 

(DFf 

4)  ■ 

-2.900000 

SCALED 

MAX  EQUALITY 

(DF* 

5)  ■ 

0.1390441 

ANALYSIS  CALLS  * 

0 

GRADIENT  CALLS  - 

DESIGN  MENU : PRO , DIS . SET , RUN , PAR ,HIS , EXP , TOL , APP , COM , QUI , HELP 
<SELECT>  PAR  GRG  C  1  1  Q 
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DESIGN  MENU:PRO,DIS, SET, RUN, PAR, HIS, EXP .TOL.APP, COM, QUI, HELP 
<SELECT>  PAR  GRA  C  1  2  turn  on  explicit 

gradients  again,  this 
is  important,  because 
finite  difference 
gradients  can  tale  a 
very  long  time. 


EDIT  GRADIENT  PARAMETERS 

VALUE  DESCRIPTION 


(  1)  2  GRADIENT  TYPE  (1-FINITE  DIFF,2=EXPLICIT) 

(  2)  0. 100000E-03  SCALED  PERTURBATION  FOR  F-D  GRADIENTS 


<EDIT:  C,G, J,Q,HELP>  Q 

DESIGN  MENU : PRO ,DIS, SET, RUN, PAR, HIS, EXP, TOL.APP, COM. QUI, HELP 
<SELECT>  RUM  GRG  5  run  5  GRG  iterations. 


SEARCHING  FOR  FEASIBLE  STARTING  POINT  note  that  ve  started 

with  a  non-feasible 
design.  The  equality 
contraint  on  the 
x  centroid  is  violated. 
Therefore,  it  becomes 
the  temporary  objective. 

TEMPORARY  OBJECTIVE  IS  TO  MINIMIZE  DESIGN  FUNCTION  5 


OBJECTIVE  FUNCTION  GRADIENT 

0.16864  0 . OOOOOE+OO  O.OOOOOE+OO  0.71870E-01  0 . OOOOOE+OO  O.OOOOOE+OO 

LAGRANGIAN  FUNCTION  GRADIENT 

O.OOOOOE+OO  O.OOOOOE+OO  O.OOOOOE+OO  O.OOOOOE+OO  O.OOOOOE+OO  O.OOOOOE+OO 


LENGTH  OF  SEARCH  DIRECTION  VECTOR-  0.21092 
SEARCH  DIRECTION 

-0.91995  0 . 10337E-03-0 . 12966E-03-0 . 39205  -0.50973E-04  0.96224E-04 


LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
FITTING  PARABOLA 
LINE  SEARCH  STEP 


1  STEP  LENGTH-  2.0000E-02 

2  STEP  LENGTH-  4.0000E-02 

3  STEP  LENGTH-  8.0000E-02 

4  STEP  LENGTH-  1.6000E-01 

5  STEP  LENGTH-  3.2000E-01 

6  STEP  LENGTH-  2.1092E-01 
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FITTING  PARABOLA 


DESIGN  NUMBER  1 


SCALED  STEP  LENGTH 

9 

0.2109229 

SCALED  OBJECTIVE 

(DFt 

1)  - 

11951.61 

SCALED  MAX  INEQUALITY 

CDF* 

2)  • 

-1.929813 

SCALED  MAX  EQUALITY 

(DF* 

5)  - 

O.OOOOOOOE+OO 

ANALYSIS  CALLS  -  6  GRADIENT  CALLS  -  1 

FEASIBLE  POINT  FOUND,  BEGINNING  OPTIMIZATION  OF  PROBLEM 

we  found 
a  feasible 

OBJECTIVE  FUNCTION  GRADIENT 

O.OOOOOE+OO  -2694.6  2436.1  O.OOOOOE+OO  1590.7 

LAGRANGIAN  FUNCTION  GRADIENT 

O.OOOOOE+OO  O.OOOOOE+OO  O.OOOOOE+OO  O.OOOOOE+OO  O.OOOOOE+OO 


LENGTH  OF  SEARCH  DIRECTION  VECTOR-  1.78223 
SEARCH  DIRECTION 

-0.20786E-04  0.56108  -0.56108  0.48774E-04  0.23567 


LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
FITTING  PARABOLA 
LINE  SEARCH  STEP 


1  STEP  LENGTH-  2.0000E-02 

2  STEP  LENGTH-  4.0000E-02 

3  STEP  LENGTH-  8.0000E-02 

4  STEP  LENGTH-  1.6000E-01 

5  STEP  LENGTH-  3.2000E-01 

6  STEP  LENGTH-  5.7000E-01 

7  STEP  LENGTH-  8.2000E-01 

8  STEP  LENGTH-  S.8527E-01 


now  wa  are 
working  on 
minimizing 
Izx. 


DESIGN  NUMBER  2 


SCALED 

STEP  LENGTH 

9 

0.5852728 

SCALED 

OBJECTIVE 

(DF# 

1) 

9 

11148.34 

SCALED 

MAX  INEQUALITY 

(DF# 

2) 

9 

-1.929752 

SCALED 

MAX  EQUALITY 

(DF# 

S) 

9 

0 . 5551115E- 

ANALYSIS  CALLS  -  14  GRADIENT  CALLS  -  2 

UPDATING  HESSIAN 
OBJECTIVE  FUNCTION  GRADIENT 

O.OOOOOE+OO  -684.75  496.63  O.OOOOOE+OO  1867.6 


LAGRANGIAN  FUNCTION  GRADIENT 
O.OOOOOE+OO  -684.75  496.63  O.OOOOOE+OO  1867.6 


LENGTH  OF  SEARCH  DIRECTION  VECTOR-  1.44810 


start 

429.65 

O.OOOOOE+OO 

■0.56111 


-396.95 

-396.95 
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0 . 61626E-04-0 . 67705 


0.39995 


SEARCH  DIRECTION 
-0.26263E-04-0. 40811  -0.46377 


LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
FITTING  PARABOLA 
LINE  SEARCH  STEP 


1  STEP 

2  STEP 

3  STEP 

4  STEP 

5  STEP 

6  STEP 


LENGTH*  2 . 0000E-02 
LENGTH*  4.0000E-02 
LENGTH-  8.0000E-02 
LENGTH*  1 . 6000E-01 
LENGTH*  3.2000E-01 
LENGTH*  5.7000E-01 


7  STEP  LENGTH-  3.3411E-01 


DESIGN  NUMBER  3 


SCALED 

STEP  LENGTH 

m 

0.3341135 

SCALED 

OBJECTIVE 

(DF* 

1) 

m 

10918.79 

SCALED 

MAX  INEQUALITY 

(DF# 

4) 

m 

-1.498272 

SCALED 

MAX  EQUALITY 

(DF* 

S) 

m 

0.8326673E- 

ANALYSIS  CALLS  *  21  GRADIENT  CALLS  *  3 

UPDATING  HESSIAN 
OBJECTIVE  FUNCTION  GRADIENT 

0 . OOOOOE+OO  -1456.9  -525.08  O.OOOOOE+OO  1265.0 


LA GRAN GIAN  FUNCTION  GRADIENT 
O.OOOOOE+OO  -1456.9  -S2S.08  O.OOOOOE+OO  1265.0 


LENGTH  OF  SEARCH  DIRECTION  VECTOR-  0.29450 
SEARCH  DIRECTION 

0.76464E-03  0.38627  0.61437  -0. 17942E-02-0. 63124 


LINE  SEARCH  STEP  1  STEP  LENGTH-  2.0000E-02 
LINE  SEARCH  STEP  2  STEP  LENGTH-  4.0000E-02 
LINE  SEARCH  STEP  3  STEP  LENGTH-  8.0000E-02 
LINE  SEARCH  STEP  4  STEP  LENGTH-  1.6000E-01 
LINE  SEARCH  STEP  5  STEP  LENGTH-  3.2000E-01 
CONSTRAINT  3  NOT  FORMERLY  BINDING 
CONVERGENCE  IN  NEVTON-RAPHSON 


DESIGN  NUMBER  4 

SCALED  STEP  LENGTH  ■  0.2944839 

SCALED  OBJECTIVE  (DF#  1)  -  10631.58 

SCALED  MAX  INEQUALITY  (DF#  4)  -  -0 . 6819350E-04 

SCALED  MAX  EQUALITY  (DF*  5)  -  0.2775558E-16 

ANALYSIS  CALLS  •  27  GRADIENT  CALLS  -  4 

UPDATING  HESSIAN 
OBJECTIVE  FUNCTION  GRADIENT 

O.OOOOOE+OO  -674.48  580.53  O.OOOOOE+OO  686.51 


45.920 

45.920 

0.27364 


711.71 
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LAGRANGIAN  FUNCTION  GRADIENT 

O.OOOOOE+OO  -674.48  580.53  O.OOOOOE+OO  686.51  211.71 


LENGTH  OF  SEARCH  DIRECTION  VECTOR- 
SEARCH  DIRECTION 
-0.39206E-02-0. 24511  -0.21417 


0.45280 

0.91998E-02-0. 24533  -0.91311 


LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
FITTING  PARABOLA 
LINE  SEARCH  STEP 


1  STEP  LENGTH-  2.0000E-02 

2  STEP  LENGTH-  4.0000E-02 

3  STEP  LENGTH-  8.0000E-02 

4  STEP  LENGTH-  1.6000E-01 

5  STEP  LENGTH-  1.0495E-01 


DESIGN  NUMBER  5 

we  have  gone  5 

SCALED  STEP  LENGTH 

■ 

0.1049469  iterations  and  are 

SCALED  OBJECTIVE 

(DF# 

1) 

- 

10614.75  not  at  the  optimum 

SCALED  MAX  INEQUALITY 

(DF* 

4) 

m 

0 . 4523932E-04 

SCALED  MAX  EQUALITY 

(DF# 

5) 

m 

0.2775558E-16 

ANALYSIS  CALLS  -  32  GRADIENT  CALLS  *  5 


DESIGN  MENU :  PRO .  DIS ,  SET , RUN , P AR , HIS  ,EXP  ,TQL ,  APP ,  COM ,  QUI ,  HELP 
<SELECT>  pro  pos  we  go  to  NAVGRAPH  to 

draw  the  current 
configuration. 


DESIGN  MENU :  PRO,  DIS,  SET,  RUN,  PAR,  HIS,  EXP,  TOL.APP,  COM,  QUI,  HELP 

<SELECT>  run  grg  3  run  3  more  iterations. 

OBJECTIVE  FUNCTION  GRADIENT 

O.OOOOOE+OO  -826.41  474.86  O.OOOOOE+OO  621.70  -56.573 

LAGRANGIAN  FUNCTION  GRADIENT 

O.OOOOOE+OO  -674.48  580.53  O.OOOOOE+OO  686.51  211.71 


LENGTH  OF  SEARCH  DIRECTION  VECTOR-  1.48955 
SEARCH  DIRECTION 

-0.35937E-02  0.48331  -0.45321  0.84326E-02  0.48286  0.57251 


LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
LINE  SEARCH  STEP 
FITTING  PARABOLA 
LINE  SEARCH  STEP 


1  STEP  LENGTH-  2.0000E-02 

2  STEP  LENGTH-  4.0000E-02 

3  STEP  LENGTH-  8.0000E-02 

4  STEP  LENGTH-  1.6000E-01 

5  STEP  LENGTH-  7.8947E-02 
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DESIGN  NUMBER  6 

SCALED  STEP  LENGTH  *  0.7894733E-01 

SCALED  OBJECTIVE  (DF#  l)  =  10601.06 

SCALED  MAX  INEQUALITY  (DFi  4)  *  0 . 2239870E-03 

SCALED  MAX  EQUALITY  (DF#  S)  ■  0 . OOOOOOOE+OO 

ANALYSIS  CALLS  -  37  GRADIENT  CALLS  -  6 

UPDATING  HESSIAN 
OBJECTIVE  FUNCTION  GRADIENT 

0 . OOOOOE+OO  -601.03  233.64  O.OOOOOE+OO  717.64  87.084 

LAGRANGIAN  FUNCTION  GRADIENT 

O.OOOOOE+OO  -601.03  233.64  O.OOOOOE+OO  717.64  87.084 

LENGTH  OF  SEARCH  DIRECTION  VECTOR-  1.01017 
SEARCH  DIRECTION 

0. 22 125E-02-0. 17708  -0.63286  -0 .51917E-02-0. 17716  -0.73261 

LINE  SEARCH  STEP  1  STEP  LENGTH-  2.0000E-02 
LINE  SEARCH  STEP  2  STEP  LENGTH-  4.0000E-02 
LINE  SEARCH  STEP  3  STEP  LENGTH-  8.0000E-02 
FITTING  PARABOLA 

LINE  SEARCH  STEP  4  STEP  LENGTH-  5.8321E-02 
DESIGN  NUMBER  7 

SCALED  STEP  LENGTH  -  0.5832070E-01 

SCALED  OBJECTIVE  (DF#  1)  •  10S94.28 

SCALED  MAX  INEQUALITY  (DF#  4)  -  0 . 2479473E-03 

SCALED  MAX  EQUALITY  (DF#  5)  •  0 . OOOOOOOE+OO 

ANALYSIS  CALLS  -  41  GRADIENT  CALLS  •  7 

UPDATING  HESSIAN 
OBJECTIVE  FUNCTION  GRADIENT 

O.OOOOOE+OO  -662.03  17.795  O.OOOOOE+OO  691.63  -22.605 

LAGRANGIAN  FUNCTION  GRADIENT 

O.OOOOOE+OO  -662.03  17.795  O.OOOOOE+OO  691.63  -22.605 

LENGTH  OF  SEARCH  DIRECTION  VECTOR-  0.53193 
SEARCH  DIRECTION 

0.56900E-02-0. 46010  -0.18861  -0. 13352E-01-0. 45932  0.73590 

LINE  SEARCH  STEP  1  STEP  LENGTH-  2.0000E-02 
LINE  SEARCH  STEP  2  STEP  LENGTH-  9.0121E-03 

DESIGN  NUMBER  8  still  not  there,  but 

SCALED  STEP  LENGTH  -  0.9012106E-02  very,  very  close. 
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SCALED  OBJECTIVE  (DF#  1)  ■  10594.13  Ihe  step  length  is 

SCALED  MAI  INEQUALITY  (DF*  4)  *  0 .2119228E-03  small  compared  to  the 

SCALED  MAX  EQUALITY  (DF#  5)  •  0 . OOOOOOOE+OO  first  iteration. 

ANALYSIS  CALLS  «  43  GRADIENT  CALLS  -  8 


DESIGN  M£NU:PRO,DIS, SET, RUN, PAR, HIS, EXP.TOL, APP.COM, QUI, HELP 
<SEL£CT>  pro  pos  go  draw  another  display. 

we  see  that  we  are 
practically  at  the 
optimum.  This  objective 
pulled  the  objects  in 
toward  the  center. 


DESIGN  MENU:PRO,DIS. SET, RUN, PAR, HIS, EXP, TOL.APP.COM, QUI, HELP 
<SELECT>  qui 
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In  this  last  problem,  we  have  a  container  wail  and  a  cylinder.  We  will  use  rota¬ 
tional  degrees  of  freedom  in  this  problem.  The  model  is  created  in  NAVGRAPH  wiih  the 
following  commands. 


DEFI 

POIN 

ENTE 

1 

0  -4  -4 
2 

0  4-4 

3 

0  4  4 

4 

-4  4 
4 

0-4  4 
0 

DR 

SURF 

ENTE 

1 

12  3  4 

TRAN 

2 

1 

10  0 

SOLI 

PARA 

1 

1  2 


DISP  AUTO-C  ON 

DR 

POIN 

ENTE 

9 

4  1  0 
0 

ROTA 

10 
9 
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A. 


315 

7 

LINE 

ARC 

13 

9  10  11 
ARC 

14 

11  12  13 
ARC 

15 

15  14  13 
ARC 

16 

9  16  15 

SUR 

LINE 

7 

13  14  15  16 
TRAN 

8 
7 

3  0  0 

SOLI 

PARA 

2 

7  8 


gene 

crea 

ala 

sol  _ 

t  create  the  finite  elements  for  the 

mass  properties  calculation. 

Only  those  elements  that  are 
associated  with  the  solids  in  the 
j  j  j  group  are  included  in  the  mass 

properties  calculation  for  the  group. 
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si 


glo  grou 
new 
add 
sol 
1 

save 

CONT 

1 

NEW 

ADD 

SOL 

2 

SAVE 
CYL 
0  1 


QROTA  Rot at*  th*  cylinder  with 

tha  quaternion  rotate  command. 

5.5  0  0 
111 
20 
DR 


GENE  CREA  ELE  SOL  mesh  the  cylinder 


10  10  1 


S2 

We  calculate  the  mass  properties  with  MASSPROP.  The  cylinder's  mass  properties 
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output  file  is  listed  here  (the  container  is  the  same  as  in  example  one). 


HAVGRAPH  DATA  FILE  *  PCYL 
GROUP  NAME  -  CYL 

t  OF  HEXAHEDRON  ELEMENTS  *  100 

t  OF  HEDGE  ELEMENTS  -  0 

t  OF  PENTAHEDRON  ELEMENTS  -  0 

«  OF  TETRAHEDRON  ELEMENTS  •*  0 


SPECIFIC  WEIGHT 

490.0000 

GRAVITY  CONST. 

32.20000 

DENSITY 

15.21739 

< 

TOTAL  VOLUME 

9.386071 

< 

TOTAL  MASS 

142.8315 

GLOBAL  XBAR 

5.500001 

< 

GLOBAL  YBAR 

1.9488257E-07 

< 

GLOBAL  ZBAR 

-3.8444448E-08 

< 

GLOBAL  IXX 

76.76131 

GLOBAL  IYY 

4459.951 

GLOBAL  IZZ 

4461.088 

GLOBAL  IXY 

-14.51823 

GLOBAL  IXZ 

12.16859 

GLOBAL  IYZ 

2.371527 

CENTROID  IXX 

76.76131 

< 

CENTROID  IYY 

139.2964 

< 

CENTROID  IZZ 

140.4336 

< 

CENTROID  IXY 

-14.51807 

< 

CENTROID  IXZ 

12.16856 

< 

CENTROID  IYZ 

2.371527 

< 

The  mass  properties  input  file,  PCYL. MAS,  is  shown  here. 


64.0  15.21739  .5  0.  0. 

10388.41  5275.362  5275.362  0.  0.  0. 

9.386064  15.21739  5.5  .0  .0 

76.76131  139.2964  140.4336  -14.518  12.1686  2.37153 

Now  we  run  PREPAC.  We  give  the  NAVGRAPH  file  name  PCYL.DAT,  the  con¬ 
tainer  name,  and  PCYL. MAS.  We  call  the  problem  file  PCYL. PRO. 

Once  PREPAC  has  calculated  the  initial  separation  distance  and  the  composite  mass 
properties,  we  run  SETUP. 
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WELCOME  TO  THE  OPTDES.BYU  VERSION  4.0  SETUP  PROGRAM 
COPYRIGHT  (C)  1988  BY  BRIGHAM  YOUNG  UNIVERSITY 


PROMPT  TYPES: 

TYPE  RESPONSE 
PRESS  RETURN  KEY 

SPECIAL  CHARACTERS: 

ESCAPE  TO  MAIN  MENU 
"X"  EXECUTE  OPERATING  SYSTEM  COMMAND 
"#"  OPEN  AND  CLOSE  INTERACTION  FILES 

INTERACTION  FILE  NAMES: 

"LOG"  STARTING  LOG  FILE  NAME 
"RES"  STARTING  RESPONSE  FILE  NAME 
"NUL"  THERE  IS  NO  STARTING  BATCH  FILE 

<NAME  OF  PROBLEM  FILE>  PCYL.PRO 

SETUP  MENU :DIS, VAR, FUN, COM, QUI, HELP 
DIS  AV 


<SELECT>  AV 

AV#  VALUE  DESCRIPTION 


1 

0 . OOOOOE+OO 

1QUAT, 

GROUP 

CYL 

2 

0 . OOOOOE+OO 

2QUAT, 

GROUP 

CYL 

3 

0. OOOOOE+OO 

3QUAT , 

GROUP 

CYL 

4 

O 

o 

w 

8 

o 

o 

o 

o 

1TRAN , 

GROUP 

CYL 

5 

o 

o 

w 

8 

O 

o 

o 

o 

2TRAN, 

GROUP 

CYL 

6 

0. OOOOOE+OO 

3TRAN, 

GROUP 

CYL 

SETUP  MENU:DIS, VAR, FUN, COM, qUI, HELP 
<SELECT>  VAR 

EDIT  DESIGN  VARIABLES 

(1)  (2)  (3)  (4)  (5) 

DV#  #MP  AV#  COEFF  VALUE  MINIMUM 


1) 

<EDIT:  C,G,J,Q,HELP>  G  1  3  2  1  1  choss  just 

EDIT  DESIGN  VARIABLES 

(1)  (2)  (3)  (4)  (5) 

DV#  #MP  AV#  COEFF  VALUE  MINIMUM 


(6) 

MAXIMUM  DESCRIPTION 


tbs  rotations 


(6) 

MAXIMUM  DESCRIPTION 
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1  l  1  1.00  0.00E+00P  -1.00  1.00  1QUAT,  GROUP  CYL 

2)  2  1  2  1.00  O.OOE+OOP  -1.00  1.00  2QUAT ,  GROUP  CYL 

3)  3  i  3  1.00  O.OOE+OOP  -1.00  1.00  3QUAT,  GROUP  CYL 

4) 

<EDIT:  C,G,J ,Q,HELP>  Q 

SETUP  MENU:DIS, VAR, FUR, COM, QUI, HELP 

<SELECT>  DIS  AF  FUN 

AF#  VALUE  DESCRIPTION 


1  2.7796  ICON  SEP  A  CYL 

2  5.0000  ICON  CEN  A  CYL 

3  73.386  TOTAL  VOLUME 

4  1116.7  TOTAL  MASS 

5  1.1395  X  MASS  CENTROID 

6  0 . OOOOOE+OO  Y  MASS  CENTROID 

7  0. OOOOOE+OO  2  MASS  CENTROID 

8  10465.  XX  MOMENT  OF  INERTIA 

9  8528.7  YY  MOMENT  OF  INERTIA 

10  8529.9  ZZ  MOMENT  OF  INERTIA 

11  -14.518  XY  MOMENT  OF  INERTIA 

12  12.169  XZ  MOMENT  OF  INERTIA 

13  2.3715  YZ  MOMENT  OF  INERTIA 

14  0. OOOOOE+OO  X  PRESSURE  CENTER 

15  0. OOOOOE+OO  Y  PRESSURE  CENTER 

16  0. OOOOOE+OO  Z  PRESSURE  CENTER 

EDIT  DESIGN  FUNCTIONS 

(1)  (2)  (3)  (4)  (5)  (6) 

DFi  iMP  AFi  COEFF  VALUE  ALLOWABLE  RANGE  DESCRIPTION 


1) 

<EDIT:  C ,G, J ,Q ,HELP>  C128C14V 

just  an  objective  function  needed 
since  there  will  be  no 
interference. 

EDIT  DESIGN  FUNCTIONS 


(1) 

DFi  «MP 

(2) 

AFi 

(3) 

COEFF 

(4)  (5) 

VALUE  ALLOWABLE 

(6) 

RANGE 

DESCRIPTION 

_ -  ... 

1)  1  1 

8 

1.00 

0 . 10E+05vP  0 . OOOE+OO 

1.00 

XX  MOMENT  OF  INERTIA 

2) 
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<EDIT:  C,G,J,Q,HELP>  Q 


SETUP  MENU :DIS, VAR, FUN, COM, QUI, HELP 
< SELECT >  QUI 

<STORE  THIS  PROBLEM? >  Y 

<NAME  OF  PROBLEM  FILE>  PCYL.PRO 


We  now  run  PACKER. 


WELCOME  TO  THE  OPTDES.BYU  VERSION  4.0  DESIGN  PROGRAM 
COPYRIGHT  (C)  1988  BY  BRIGHAM  YOUNG  UNIVERSITY 

PROMPT  TYPES: 

TYPE  RESPONSE 
PRESS  RETURN  KEY 

SPECIAL  CHARACTERS: 

"*"  ESCAPE  TO  MAIN  MENU 
"X"  EXECUTE  OPERATING  SYSTEM  COMMAND 
OPEN  AND  CLOSE  INTERACTION  FILES 

INTERACTION  FILE  NAMES: 

"LOG"  STARTING  LOG  FILE  NAME 
"RES"  STARTING  RESPONSE  FILE  NAME 
"NUL"  THERE  IS  NO  STARTING  BATCH  FILE 

<NAME  OF  PROBLEM  FILE>  PCYL.PRO 


****The  PACKER  BANNER*****  (not  shown  here  to  save  space) 

Input  the  NAVGRAPH  file  (PCYL.DAT),  the  container  name,  and  the  mass  properties 
file  (PCYL.MAS).  The  initial  design  looked  like  this: 
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DESIGH  HUMBER  0 

SCALED  OBJECTIVE  (DF#  1)  *  10465.17 

AHALYSIS  CALLS  *  0  GRADIENT  CALLS  *  0 

DESIGH  MENU : PRO, DIS, SET .RUN, PAR, HIS, EXP, TOL.APP, COM, QUI, HELP 
<SELECT>  PAR  GRG  C  1  1  Q 

DESIGN  MENU:PRO, DIS. SET, RUN. PAR, HIS, EXP, TOL.APP, COM, QUI, HELP 
<SELECT>  PAR  GRA  C  1  2  Q  turn  on  explicit  gradients 

DESIGN  MENU:PRO, DIS, SET, RUN, PAR, HIS, EXP.TOL, APP.COM, QUI, HELP 
<SELECT>  DIS  DV  DIS  DF  display  design  variables 

DV*  SCALED  VALUE  UNSCALED  VALUE  SCALED  MIN  SCALED  MAX  DESCRIPTION 

1  0 . OOOOOOOE+OO  O.OOOOOOOE+OO  -1.00000  1.00000  1QUAT,  GROUP  CYL 

2  O.OOOOOOOE+OO  O.OOOOOOOE+OO  -1.00000  1.00000  2QUAT,  GROUP  CYL 

3  O.OOOOOOOE+OO  O.OOOOOOOE+OO  -1.00000  1.00000  3QUAT ,  GROUP  CYL 

DFi  SCALED  VALUE  UNSCALED  VALUE  DESCRIPTION 


1  10465.17  10465.17  XX  MOMENT  OF  INERTIA 

DESIGN  MENU : PRO , DIS . SET , RUN , PAR , HIS , EXP , TOL , APP , COM , QUI , HELP 
<SELECT>  RUN  GRG  5  run  5  GRG  iterations 

OBJECTIVE  FUNCTION  GRADIENT 
-0.35527E-12  7701.1  7705.8 

LAGRANGIAN  FUNCTION  GRADIENT 
0 . OOOOOE+OO  0 . OOOOOE+OO  O.OOOOOE+OO 
LENGTH  OF  SEARCH  DIRECTION  VECTOR*  1.41421 
SEARCH  DIRECTION 
0.48916E- 16-0. 70711  -0.70711 

LINE  SEARCH  STEP  1  STEP  LENGTH*  2.0000E-02 

LINE  SEARCH  STEP  2  STEP  LENGTH*  4.0000E-02  searching  for  the 

LINE  SEARCH  STEP  3  STEP  LENGTH*  8.0000E-02  optimum 

LINE  SEARCH  STEP  4  STEP  LENGTH*  1.6000E-01 

LINE  SEARCH  STEP  5  STEP  LENGTH*  3.2000E-01 

LINE  SEARCH  STEP  6  STEP  LENGTH*  5.7000E-01 
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FITTING  PARABOLA 

LINE  SEARCH  STEP  7  STEP  LENGTH®  2.7876E-01 
DESIGN  NUMBER  1 

SCALED  STEP  LENGTH  »  0.2787S71 

SCALED  OBJECTIVE  (DF#  1)  -  10459.79 

ANALYSIS  CALLS  -  7  GRADIENT  CALLS  *  1 

UPDATING  HESSIAN 
OBJECTIVE  FUNCTION  GRADIENT 
0.49921  7135.7  7140.6 

LAGRANGIAN  FUNCTION  GRADIENT 
0.49921  7135.7  7140.6 

LENGTH  OF  SEARCH  DIRECTION  VECTOR-  1.13546 
SEARCH  DIRECTION 

-0.43432E-03-0.70711  -0.70711 

LINE  SEARCH  STEP  1  STEP  LENGTH-  2.0000E-02 

LINE  SEARCH  STEP  2  STEP  LENGTH-  9.9985E-03 

DESIGN  NUMBER  2 

SCALED  STEP  LENGTH  -  0 . OOOOOOOE+OO 

SCALED  OBJECTIVE  (DF#  1)  -  104S9.79 

ANALYSIS  CALLS  -  10  GRADIENT  CALLS  -  2 

OPTIMUM  ASSUMED  (MINIMUM  STEPS IZE  REACHED)  optimum  reached. 

DESIGN  MENU: PRO , DIS , SET , RUN , PAR , HIS , EXP , TOL , APP , COM , QUI , HELP 
<SELECT>  PRO  POS  draw  the  optimal  solution. 

the  cylinder  has  indeed  rotated 
so  that  the  axis  of  symmetry  lines 
up  parallel  sith  the  global  x  axis 


Trouble  Shooting  Guide  For  NAVGRAPH  Model  Generation 


PROBLEM: 

Points  appear  to  be  coincident  and  cannot  be  selected 
accurately  with  thumb  wheels  or  mouse  select. 

POSSIBLE  SOLUTIONS: 

-rotate  and/or  magnify  until  unique  point  can  be  selected. 
PROBLEM: 

Point  Labels  cannot  be  seen  or  read. 

POSSIBLE  SOLUTIONS: 

-check  to  see  if  point  labels  axe  turned  on  (DISPLAY  OPTIONS) 
-turn  point  entities  on 

-rotate  and/or  magnify  (first  turn  AUTO-CENT  off) 

-if  the  NAVGRAPH  data  base  you  are  working  with  was  translated 
from  an  AFATL/ASE  format  data  base  turn  node  entities  on,  node 
labels  on,  and  HIDDEN  on  (under  DISPLAY  OPTIONS).  Node  numbers 
and  node  coordinates  are  the  same  as  the  point  numbers  and 
point  coordinates 

-plot  selected  points  by  using  the  PLOT  POINT  command 
PROBLEM: 

A  line  cannot  be  defined. 

POSSIBLE  SOLUTIONS: 

-check  to  see  if  each  point  (not  node)  defining  the  line  is 
defined  in  the  data  base  by  using  the  LIST  POINT  command 
-be  sure  the  correct  number  of  points  are  given  for  the  line 
-no  more  than  30  points  may  be  used  to  define  a  spline 

PROBLEM: 

A  line  or  spline  does  not  follow  the  shape  of  the  surface  or 
edge  closely  enough. 

POSSIBLE  SOLUTIONS: 

-redefine  the  line  by  overwriting 

-use  a  different  type  of  line 

-try  using  fewer,  more,  or  different  points 

-define  additional  points  on  the  surface  or  edge 

-consider  a  different  location  for  the  edge 

-consider  using  degenerate  edges 

PROBLEM: 

A  valid  solid  cannot  be  defined 
POSSIBLE  SOLUTIONS: 

-make  sure  that  the  lines  defining  the  solid  share  eight  (8) 
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7 


* 


y3 


Figure  18:  The  edge  direction  and  ordering  convention  for  NAVGRAPH  solids. 

common  end  points  (at  the  corners) 

-make  sure  that  the  solid  is  composed  of  twelve  (12)  lines 
-assure  that  line  directions  are  as  shown  in  Figure  17 
-define  surfaces  one  at  a  time  to  identify  problem  points  or 
lines 

PROBLEM: 

A  draw  command  does  not  draw  desired  entities 
POSSIBLE  SOLUTIONS: 

-list  entities  to  make  sure  they  axe  in  the  data  base 
-turn  appropriate  entities  on 
-try  drawing  points  only  first 

PROBLEM: 

A  plot  command  does  not  respond. 

POSSIBLE  SOLUTIONS: 

-check  to  make  sure  desired  entities  are  present  in  data  base 
-enter  a  draw  command  first  followed  by  a  plot  command 
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Suggestions  and  Tips 


1.  Before  creating  the  model,  make  a  rough  sketch  on  paper  of 
each  line  including  direction,  line,  and  point  labels. 

2.  To  save  time,  take  advantage  of  EXTRUDE  and  SWEEP  commands. 

3.  Creating  half  of  a  model  that  is  symmetric  around  at  least  one 
axis  is  advisable.  The  model  can  then  be  reflected  using  the 
MIRROR  command.  This  saves  time  and  is  usually  more  accurate. 

4.  Merging  two  or  more  lines  can  be  helpful  in  curve  fitting. 

5.  Use  the  surface  fill  option  (under  MORE  OPTIONS)  to  see  if  a 
solid  follows  the  shape  of  the  desired  object. 

6.  Become  as  familiar  an  possible  with  the  shape  of  the  object  to 
be  modelled  in  order  to  choose  the  best  points  for  line  and 
spline  definition. 

7.  Try  to  define  lines  using  points  that  are  equally  spaced. 

8.  3-D  graphics  can  be  deceiving  at  time  --  rotate  and/or  magnify 
to  verify  the  correctness  of  the  geometry.  (Listing  the  data 
base  may  not  always  be  helpful.) 

9.  To  get  a  better  idea  of  the  shape  of  the  object  to  be 
modelled,  turn  HIDDEN  on  and  node  and  element  labels  off. 

View  from  several  different  angles. 

10.  View  the  shape  as  a  deformed  cube  of  several  deformed  cubes  and 
try  to  place  twelve  (12)  lines  in  a  natural  way. 

11.  Time  can  be  saved  by  turning  HIDDEN  off  when  executing  a  draw 
command  (the  hidden  line  algorithm  is  by-passed) . 
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Curve  Fitting  Technique 

Figure  18  shows  eight  (8)  points  around  the  nose  section  of  the  F-15.  Points  1, 
3,  4,  5,  6,  9,  8,  and  7  axe  at  the  surface  and  when  joined  by  straight  lines  show  the 
general  shape.  The  other  number  labels  on  the  figure  label  the  finite  elements.  A  spline 
is  best  suited  to  fit  the  shape.  At  the  appropriate  prompt  11  point  numbers  are  given 
in  this  order:  7,  1,  3,  4,  5,  6,  9,  8,  7,  1,  and  3.  This  method  of  overlapping  beginning 
and  ending  points  is  necessary  to  give  a  smooth  curve  through  point  1  (see  Figure  19). 
This  procedure  creates  11  lines  as  follows:  line  1  is  the  entire  spline  (it  is  defined  as  a 
merged  line  in  the  data  base),  lines  2-11  are  the  individual  parametric  cubic  definition 
lines  created  between  each  pair  of  points.  Lines  1 ,  2,  and  1 1  can  then  be  deleted  to  leave 
only  one  line  between  each  pair  of  points  (see  Figure  20).  These  lines  are  now  available 
for  merging  and  creating  solids. 
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